MBNL1 alternative splicing isoforms play opposing roles in cancer

RBFOX2 deregulation promotes pancreatic cancer progression and metastasis through alternative splicing

RNA splicing is an important biological process associated with cancer initiation and progression. However, the contribution of alternative splicing to pancreatic cancer (PDAC) development is not well understood. Here, we identify an enrichment of RNA binding proteins (RBPs) involved in splicing regulation linked to PDAC progression from a forward genetic screen using Sleeping Beauty insertional mutagenesis in a mouse model of pancreatic cancer. We demonstrate downregulation of RBFOX2, an RBP of the FOX family, promotes pancreatic cancer progression and liver metastasis. Specifically, we show RBFOX2 regulates exon splicing events in transcripts encoding proteins involved in cytoskeletal remodeling programs. These exons are differentially spliced in PDAC patients, with enhanced exon skipping in the classical subtype for several RBFOX2 targets. RBFOX2 mediated splicing of ABI1, encoding the Abelson-interactor 1 adapter protein, controls the abundance and localization of ABI1 protein isoforms in pancreatic cancer cells and promotes the relocalization of ABI1 from the cytoplasm to the periphery of migrating cells. Using splice-switching antisense oligonucleotides (AONs) we demonstrate the ABI1 ∆Ex9 isoform enhances cell migration. Together, our data identify a role for RBFOX2 in promoting PDAC progression through alternative splicing regulation.

FMRP deficiency leads to multifactorial dysregulation of splicing and mislocalization of MBNL1 to the cytoplasm

Fragile X syndrome (FXS) is a neurodevelopmental disorder that is often modeled in Fmr1 knockout mice where the RNA-binding protein FMRP is absent. Here, we show that in Fmr1-deficient mice, RNA mis-splicing occurs in several brain regions and peripheral tissues. To assess molecular mechanisms of splicing mis-regulation, we employed N2A cells depleted of Fmr1. In the absence of FMRP, RNA-specific exon skipping events are linked to the splicing factors hnRNPF, PTBP1, and MBNL1. FMRP regulates the translation of Mbnl1 mRNA as well as Mbnl1 RNA auto-splicing. Elevated Mbnl1 auto-splicing in FMRP-deficient cells results in the loss of a nuclear localization signal (NLS)-containing exon. This in turn alters the nucleus-to-cytoplasm ratio of MBNL1. This redistribution of MBNL1 isoforms in Fmr1-deficient cells could result in downstream splicing changes in other RNAs. Indeed, further investigation revealed that splicing disruptions resulting from Fmr1 depletion could be rescued by overexpression of nuclear MBNL1. Altered Mbnl1 auto-splicing also occurs in human FXS postmortem brain. These data suggest that FMRP-controlled translation and RNA processing may cascade into a general dys-regulation of splicing in Fmr1-deficient cells.

Modulatory role of RNA helicases in MBNL-dependent alternative splicing regulation

Muscleblind-like splicing regulators (MBNLs) activate or repress the inclusion of alternative splicing (AS) events, enabling the developmental transition of fetal mRNA splicing isoforms to their adult forms. Herein, we sought to elaborate the mechanism by which MBNLs mediate AS related to biological processes. We evaluated the functional role of DEAD-box (DDX) RNA helicases, DDX5 and DDX17 in MBNL-dependent AS regulation. Whole-transcriptome analysis and validation approaches revealed a handful of MBNLs-dependent AS events to be affected by DDX5 and DDX17 in mostly an opposite manner. The opposite expression patterns of these two groups of factors during muscle development and coordination of fetal-to-adult splicing transition indicate the importance of these proteins at early stages of development. The identified pathways of how the helicases modulate MBNL splicing activity include DDX5 and DDX17-dependent changes in the ratio of MBNL splicing isoforms and most likely changes in accessibility of MBNL-binding sites. Another pathway involves the mode of action of the helicases independent of MBNL activity. These findings lead to a deeper understanding of the network of interdependencies between RNA-binding proteins and constitute a valuable element in the discussion on developmental homeostasis and pathological states in which the studied protein factors play a significant role.

Regulation of Epithelial-Mesenchymal Transitions by Alternative Splicing: Potential New Area for Cancer Therapeutics

The epithelial-mesenchymal transition (EMT) is a complicated biological process in which cells with epithelial phenotype are transformed into mesenchymal cells with loss of cell polarity and cell-cell adhesion and gain of the ability to migrate. EMT and the reverse mesenchymal-epithelial transitions (METs) are present during cancer progression and metastasis. Using the dynamic switch between EMT and MET, tumour cells can migrate to neighbouring organs or metastasize in the distance and develop resistance to traditional chemotherapy and targeted drug treatments. Growing evidence shows that reversing or inhibiting EMT may be an advantageous approach for suppressing the migration of tumour cells or distant metastasis. Among different levels of modulation of EMT, alternative splicing (AS) plays an important role. An in-depth understanding of the role of AS and EMT in cancer is not only helpful to better understand the occurrence and regulation of EMT in cancer progression, but also may provide new therapeutic strategies. This review will present and discuss various splice variants and splicing factors that have been shown to play a crucial role in EMT.

Muscleblind-like proteins use modular domains to localize RNAs by riding kinesins and docking to membranes

RNA binding proteins (RBPs) act as critical facilitators of spatially regulated gene expression. Muscleblind-like (MBNL) proteins, implicated in myotonic dystrophy and cancer, localize RNAs to myoblast membranes and neurites through unknown mechanisms. We find that MBNL forms motile and anchored granules in neurons and myoblasts, and selectively associates with kinesins Kif1bα and Kif1c through its zinc finger (ZnF) domains. Other RBPs with similar ZnFs associate with these kinesins, implicating a motor-RBP specificity code. MBNL and kinesin perturbation leads to widespread mRNA mis-localization, including depletion of Nucleolin transcripts from neurites. Live cell imaging and fractionation reveal that the unstructured carboxy-terminal tail of MBNL1 allows for anchoring at membranes. An approach, termed RBP Module Recruitment and Imaging (RBP-MRI), reconstitutes kinesin- and membrane-recruitment functions using MBNL-MS2 coat protein fusions. Our findings decouple kinesin association, RNA binding, and membrane anchoring functions of MBNL while establishing general strategies for studying multi-functional, modular domains of RBPs.

Genome-wide landscape of RNA-binding protein (RBP) networks as potential molecular regulators of psychiatric co-morbidities: a computational analysis

Background

Psychiatric disorders are a major burden on global health. These illnesses manifest as co-morbid conditions, further complicating the treatment. There is a limited understanding of the molecular and regulatory basis of psychiatric co-morbidities. The existing research in this regard has largely focused on epigenetic modulators, non-coding RNAs, and transcription factors. RNA-binding proteins (RBPs) functioning as multi-protein complexes are now known to be predominant controllers of multiple gene regulatory processes. However, their involvement in gene expression dysregulation in psychiatric co-morbidities is yet to be understood.

Results

Ten RBPs (QKI, ELAVL2, EIF2S1, SRSF3, IGF2BP2, EIF4B, SNRNP70, FMR1, DAZAP1, and MBNL1) were identified to be associated with psychiatric disorders such as schizophrenia, major depression, and bipolar disorders. Analysis of transcriptomic changes in response to individual depletion of these RBPs showed the potential influence of a large number of RBPs driving differential gene expression, suggesting functional cross-talk giving rise to multi-protein networks. Subsequent transcriptome analysis of post-mortem human brain samples from diseased and control individuals also suggested the involvement of ~ 100 RBPs influencing gene expression changes. These RBPs were found to regulate various processes including transcript splicing, mRNA transport, localization, stability, and translation. They were also found to form an extensive interactive network. Further, hnRNP, SRSF, and PCBP family RBPs, Matrin3, U2AF2, KHDRBS1, PTBP1, and also PABPN1 were found to be the hub proteins of the RBP network.

Conclusions

Extensive RBP networks involving a few hub proteins could result in transcriptome-wide dysregulation of post-transcriptional modifications, potentially driving multiple psychiatric disorders. Understanding the functional involvement of RBP networks in psychiatric disorders would provide insights into the molecular basis of psychiatric co-morbidities.

RNA Foci Formation in a Retinal Glial Model for Spinocerebellar Ataxia Type 7

Spinocerebellar ataxia type 7 (SCA7) is a neurodegenerative disorder characterized by cerebellar ataxia and retinopathy. SCA7 is caused by a CAG expansion in the ATXN7 gene, which results in an extended polyglutamine (polyQ) tract in the encoded protein, the ataxin-7. PolyQ expanded ataxin-7 elicits neurodegeneration in cerebellar Purkinje cells, however, its impact on the SCA7-associated retinopathy remains to be addressed. Since Müller glial cells play an essential role in retinal homeostasis, we generate an inducible model for SCA7, based on the glial Müller MIO-M1 cell line. The SCA7 pathogenesis has been explained by a protein gain-of-function mechanism, however, the contribution of the mutant RNA to the disease cannot be excluded. In this direction, we found nuclear and cytoplasmic foci containing mutant RNA accompanied by subtle alternative splicing defects in MIO-M1 cells. RNA foci were also observed in cells from different lineages, including peripheral mononuclear leukocytes derived from SCA7 patient, suggesting that this molecular mark could be used as a blood biomarker for SCA7. Collectively, our data showed that our glial cell model exhibits the molecular features of SCA7, which makes it a suitable model to study the RNA toxicity mechanisms, as well as to explore therapeutic strategies aiming to alleviate glial dysfunction.

Functional analysis of the short splicing variant encoded by CHI3L1/YKL-40 in glioblastoma

The glycoprotein YKL-40 has been well studied as a serum biomarker of prognosis and disease status in glioblastoma. YKL-40 is a chitinase-like protein with defective chitinase activity that plays an important role in promoting cell proliferation, migration, and metastasis in glioblastoma multiforme (GBM). The short variant (SV) of YKL-40, generated by an alternative splicing event that splices out exon 8, was reported in the early developing human musculoskeletal system, although its role in GBM is still unknown. Our results showed that individual glioblastoma cell lines displayed increased expression of the short variant of YKL-40 after low serum treatment. In addition, unlike the full-length (FL) version, which was localized to all cell compartments, the short isoform could not be secreted and was localized only to the cytoplasm. Functionally, FL YKL-40 promoted cell proliferation and migration, whereas SV YKL-40 suppressed them. Transcriptome analysis revealed that these opposing roles of the two isoforms may be modulated by differentially regulating several oncogenic-related pathways, including p53, the G2/M checkpoint, and MYC-related signaling. This study may provide new ideas for the development of targeted anti-YKL-40 therapy in GBM treatment.

An interplay between BRD4 and G9a regulates skeletal myogenesis

Histone acetylation and methylation are epigenetic modifications that are dynamically regulated by chromatin modifiers to precisely regulate gene expression. However, the interplay by which histone modifications are synchronized to coordinate cellular differentiation is not fully understood. In this study, we demonstrate a relationship between BRD4, a reader of acetylation marks, and G9a, a writer of methylation marks in the regulation of myogenic differentiation. Using loss- and gain-of-function studies, as well as a pharmacological inhibition of its activity, we examined the mechanism by which BRD4 regulates myogenesis. Transcriptomic analysis using RNA sequencing revealed that a number of myogenic differentiation genes are downregulated in Brd4-depleted cells. Interestingly, some of these genes were upregulated upon G9a knockdown, indicating that BRD4 and G9a play opposing roles in the control of myogenic gene expression. Remarkably, the differentiation defect caused by Brd4 knockdown was rescued by inhibition of G9a methyltransferase activity. These findings demonstrate that the absence of BRD4 results in the upregulation of G9a activity and consequently impaired myogenic differentiation. Collectively, our study identifies an interdependence between BRD4 and G9a for the precise control of transcriptional outputs to regulate myogenesis.

The MBNL1/circNTRK2/PAX5 pathway regulates aerobic glycolysis in glioblastoma cells by encoding a novel protein NTRK2-243aa

Glioblastoma multiforme (GBM) is the most common tumor of the human central nervous system. Aerobic glycolysis has been strongly related to tumor development and malignant behavior. In this study, we found that MBNL1, circNTRK2, and NTRK2-243aa were markedly downregulated and inhibited glycolysis in GBM, whereas PAX5 was upregulated and promoted glycolysis. Functionally, MBNL1 promoted the expression of circNTRK2 by binding to NTRK2 pre-mRNA, as validated using RNA pull-down and nascent RNA immunoprecipitation assays. Mass spectrometry, western blotting, and immunofluorescence staining methods were used to detect the expression of NTRK2-243aa. NTRK2-243aa-encoded by circNTRK2-phosphorylated PAX5 at Y102, leading to the attenuation of the half-life of PAX5, as validated by in vitro kinase and MG132 rescue assays. Besides, PAX5 transcriptionally facilitated the expression of PKM2 and HK2 by binding to their promoter regions, as verified by luciferase reporter and chromatin immunoprecipitation assays. Finally, overexpression of MBNL1 and circNTRK2 combined with PAX5 knockdown effectively inhibited the formation of GBM xenograft tumors and significantly prolonged the survival of orthotopic nude mice. We have delineated that the MBNL1/circNTRK2/PAX5 pathway plays a crucial role in regulating GBM glycolysis and could provide potential targets and alternative strategies for the treatment of GBM.

DEAD-Box RNA Helicases DDX3X and DDX5 as Oncogenes or Oncosuppressors: A Network Perspective

Simple Summary

The transformation of a normal cell into a cancerous one is caused by the deregulation of different metabolic pathways, involving a complex network of protein–protein interactions. The cellular enzymes DDX3X and DDX5 play important roles in the maintenance of normal cell metabolism, but their deregulation can accelerate tumor transformation. Both DDX3X and DDX5 interact with hundreds of different cellular proteins, and depending on the specific pathways in which they are involved, both proteins can either act as suppressors of cancer or as oncogenes. In this review, we summarize the current knowledge about the roles of DDX3X and DDX5 in different tumors. In addition, we present a list of interacting proteins and discuss the possible contribution of some of these protein–protein interactions in determining the roles of DDX3X and DDX5 in the process of cancer proliferation, also suggesting novel hypotheses for future studies.

Abstract

RNA helicases of the DEAD-box family are involved in several metabolic pathways, from transcription and translation to cell proliferation, innate immunity and stress response. Given their multiple roles, it is not surprising that their deregulation or mutation is linked to different pathological conditions, including cancer. However, while in some cases the loss of function of a given DEAD-box helicase promotes tumor transformation, indicating an oncosuppressive role, in other contexts the overexpression of the same enzyme favors cancer progression, thus acting as a typical oncogene. The roles of two well-characterized members of this family, DDX3X and DDX5, as both oncogenes and oncosuppressors have been documented in several cancer types. Understanding the interplay of the different cellular contexts, as defined by the molecular interaction networks of DDX3X and DDX5 in different tumors, with the cancer-specific roles played by these proteins could help to explain their apparently conflicting roles as cancer drivers or suppressors.

High miR203a-3p and miR-375 expression in the airways of smokers with and without COPD

Smoking is a leading cause of chronic obstructive pulmonary disease (COPD). It is known to have a significant impact on gene expression and (inflammatory) cell populations in the airways involved in COPD pathogenesis. In this study, we investigated the impact of smoking on the expression of miRNAs in healthy and COPD individuals. We aimed to elucidate the overall smoking-induced miRNA changes and those specific to COPD. In addition, we investigated the downstream effects on regulatory gene expression and the correlation to cellular composition. We performed a genome-wide miRNA expression analysis on a dataset of 40 current- and 22 ex-smoking COPD patients and a dataset of 35 current- and 38 non-smoking respiratory healthy controls and validated the results in an independent dataset. miRNA expression was then correlated with mRNA expression in the same patients to assess potential regulatory effects of the miRNAs. Finally, cellular deconvolution analysis was used to relate miRNAs changes to specific cell populations. Current smoking was associated with increased expression of three miRNAs in the COPD patients and 18 miRNAs in the asymptomatic smokers compared to respiratory healthy controls. In comparison, four miRNAs were lower expressed with current smoking in asymptomatic controls. Two of the three smoking-related miRNAs in COPD, miR-203a-3p and miR-375, were also higher expressed with current smoking in COPD patients and the asymptomatic controls. The other smoking-related miRNA in COPD patients, i.e. miR-31-3p, was not present in the respiratory healthy control dataset. miRNA-mRNA correlations demonstrated that miR-203a-3p, miR-375 and also miR-31-3p expression were negatively associated with genes involved in pro-inflammatory pathways and positively associated with genes involved in the xenobiotic pathway. Cellular deconvolution showed that higher levels of miR-203a-3p were associated with higher proportions of proliferating-basal cells and secretory (club and goblet) cells and lower levels of fibroblasts, luminal macrophages, endothelial cells, B-cells, amongst other cell types. MiR-375 expression was associated with lower levels of secretory cells, ionocytes and submucosal cells, but higher levels of endothelial cells, smooth muscle cells, and mast cells, amongst other cell types. In conclusion, we identified two smoking-induced miRNAs (miR-375 and miR-203a-3p) that play a role in regulating inflammation and detoxification pathways, regardless of the presence or absence of COPD. Additionally, in patients with COPD, we identified miR-31-3p as a miRNA induced by smoking. Our identified miRNAs should be studied further to unravel which smoking-induced inflammatory mechanisms are reactive and which are involved in COPD pathogenesis.

Disrupting the Molecular Pathway in Myotonic Dystrophy

Myotonic dystrophy is the most common muscular dystrophy in adults. It consists of two forms: type 1 (DM1) and type 2 (DM2). DM1 is associated with a trinucleotide repeat expansion mutation, which is transcribed but not translated into protein. The mutant RNA remains in the nucleus, which leads to a series of downstream abnormalities. DM1 is widely considered to be an RNA-based disorder. Thus, we consider three areas of the RNA pathway that may offer targeting opportunities to disrupt the production, stability, and degradation of the mutant RNA.

Targeted SMN Exon Skipping: A Useful Control to Assess In Vitro and In Vivo Splice-Switching Studies

The literature surrounding the use of antisense oligonucleotides continues to grow, with new disease and mechanistic applications constantly evolving. Furthermore, the discovery and advancement of novel chemistries continues to improve antisense delivery, stability and effectiveness. For each new application, a rational sequence design is recommended for each oligomer, as is chemistry and delivery optimization. To confirm oligomer delivery and antisense activity, a positive control AO sequence with well characterized target-specific effects is recommended. Here, we describe splice-switching antisense oligomer sequences targeting the ubiquitously expressed human and mouse SMN and Smn genes for use as control AOs for this purpose. We report two AO sequences that induce targeted skipping of SMN1/SMN2 exon 7 and two sequences targeting the Smn gene, that induce skipping of exon 5 and exon 7. These antisense sequences proved effective in inducing alternative splicing in both in vitro and in vivo models and are therefore broadly applicable as controls. Not surprisingly, we discovered a number of differences in efficiency of exon removal between the two species, further highlighting the differences in splice regulation between species.

Altered Expression of MBNL Family of Alternative Splicing Factors in Colorectal Cancer

BACKGROUND/AIM

Colorectal cancer is currently the third leading cause of cancer-related deaths and recently, alternative splicing has risen as its important regulator and potential treatment target. In the present study, we analyzed gene expression of the MBNL family of regulators of alternative splicing in various stages of colorectal cancer development, together with the MBNL-target splicing events in FOXP1 and EPB41L3 genes and tumor-related CD44 variants.

MATERIALS AND METHODS

Samples of tumor tissue and non-malignant mucosa from 108 patients were collected. After RNA isolation and reverse transcription, the relative gene expression of a selected gene panel was tested by quantitative real-time PCR, followed by statistical analysis.

RESULTS

MBNL expression was decreased in tumor tissue compared to non-tumor mucosa. In addition, lower expression was observed for the variants of FOXP1 and EPB41L3, while higher expression in tumor tissue was detected both for total CD44 and its cancer-related variants 3 and 6. Transcript levels of the MBNL genes were not found to be related to any of the studied clinicopathological characteristics. Multiple significant associations were identified in the target gene panel, including higher transcript levels of FOXP1 and CD44v3 in patients with distant metastases and connections between recurrence-free survival and altered levels of FOXP1 and CD44v3.

CONCLUSION

Our results identified for the first-time deregulation of MBNL genes in colorectal cancer. Down-regulation of their transcripts in tumor tissue compared to matched non-tumor mucosa can lead to transition of alternative splicing patterns towards a less differentiated phenotype, which highlights the importance of alternative splicing regulation for tumor growth and propagation.

A cell-to-patient machine learning transfer approach uncovers novel basal-like breast cancer prognostic markers amongst alternative splice variants

BACKGROUND

Breast cancer is amongst the 10 first causes of death in women worldwide. Around 20% of patients are misdiagnosed leading to early metastasis, resistance to treatment and relapse. Many clinical and gene expression profiles have been successfully used to classify breast tumours into 5 major types with different prognosis and sensitivity to specific treatments. Unfortunately, these profiles have failed to subclassify breast tumours into more subtypes to improve diagnostics and survival rate. Alternative splicing is emerging as a new source of highly specific biomarkers to classify tumours in different grades. Taking advantage of extensive public transcriptomics datasets in breast cancer cell lines (CCLE) and breast cancer tumours (TCGA), we have addressed the capacity of alternative splice variants to subclassify highly aggressive breast cancers.

RESULTS

Transcriptomics analysis of alternative splicing events between luminal, basal A and basal B breast cancer cell lines identified a unique splicing signature for a subtype of tumours, the basal B, whose classification is not in use in the clinic yet. Basal B cell lines, in contrast with luminal and basal A, are highly metastatic and express epithelial-to-mesenchymal (EMT) markers, which are hallmarks of cell invasion and resistance to drugs. By developing a semi-supervised machine learning approach, we transferred the molecular knowledge gained from these cell lines into patients to subclassify basal-like triple negative tumours into basal A- and basal B-like categories. Changes in splicing of 25 alternative exons, intimately related to EMT and cell invasion such as ENAH, CD44 and CTNND1, were sufficient to identify the basal-like patients with the worst prognosis. Moreover, patients expressing this basal B-specific splicing signature also expressed newly identified biomarkers of metastasis-initiating cells, like CD36, supporting a more invasive phenotype for this basal B-like breast cancer subtype.

CONCLUSIONS

Using a novel machine learning approach, we have identified an EMT-related splicing signature capable of subclassifying the most aggressive type of breast cancer, which are basal-like triple negative tumours. This proof-of-concept demonstrates that the biological knowledge acquired from cell lines can be transferred to patients data for further clinical investigation. More studies, particularly in 3D culture and organoids, will increase the accuracy of this transfer of knowledge, which will open new perspectives into the development of novel therapeutic strategies and the further identification of specific biomarkers for drug resistance and cancer relapse.

Splice-Switching Antisense Oligonucleotides as a Targeted Intrinsic Engineering Tool for Generating Armored Redirected T Cells

Modification of specificity of T cells for the use in adoptive transfer (CAR- or TCR-redirected T cells) has revolutionized the therapy of liquid tumors and some infectious diseases. However, several obstacles are still hampering the efficacy of such potent therapy, hence concurrent modification of the function is also required to obtain successful results. Here we show the use of splice-switching antisense oligonucleotides (SSOs) as a tool to transiently modify T cell function. We demonstrate the possibility to transfect SSOs and an exogenous TCR into primary human T cells in the same electroporation reaction, without affecting viability and function of the transfected T lymphocytes. Moreover, we show that SSOs targeting T cell-specific mRNAs induce the skipping of the targeted exons, and the reduction of the protein and consequent modification of T cell function. This technical work paves the way to the use of SSOs in immune cells, not only for the knockdown of the functional isoform of the targeted proteins, but also for the protein manipulation by elimination of specific domains encoded by targeted exons.

Identification of RNA-binding protein SNRPA1 for prognosis in prostate cancer

Prostate cancer is one of the deadliest cancers in men. RNA-binding proteins play a critical role in human cancers; however, whether they have a significant effect on the prognosis of prostate cancer has yet to be elucidated. In the present study, we performed a comprehensive analysis of RNA sequencing and clinical data from the Cancer Genome Atlas dataset and obtained differentially expressed RNA-binding proteins between prostate cancer and benign tissues. We constructed a protein-protein interaction network and Cox regression analyses were conducted to identify prognostic hub RNA-binding proteins. SNRPA1 was associated with the highest risk of poor prognosis and was therefore selected for further analysis. SNRPA1 expression was positively correlated with Gleason score and pathological TNM stage in prostate cancer patients. Furthermore, the expression profile of SNRPA1 was validated using the Oncomine, Human Protein Atlas, and Cancer Cell Line Encyclopedia databases. Meanwhile, the prognostic profile of SNRPA1 was successfully verified in GSE70769. Additionally, the results of molecular experiments revealed the proliferative role of SNRPA1 in prostate cancer cells. In summary, our findings evidenced a relationship between RNA-binding proteins and prostate cancer and indicated the prognostic significance of SNRPA1 in prostate cancer.

Splicing factors: Insights into their regulatory network in alternative splicing in cancer

More than 95% of all human genes are alternatively spliced after transcription, which enriches the diversity of proteins and regulates transcript and/or protein levels. The splicing isoforms produced from the same gene can manifest distinctly, even exerting opposite effects. Mounting evidence indicates that the alternative splicing (AS) mechanism is ubiquitous in various cancers and drives the generation and maintenance of various hallmarks of cancer, such as enhanced proliferation, inhibited apoptosis, invasion and metastasis, and angiogenesis. Splicing factors (SFs) play pivotal roles in the recognition of splice sites and the assembly of spliceosomes during AS. In this review, we mainly discuss the similarities and differences of SF domains, the details of SF function in AS, the effect of SF-driven pathological AS on different hallmarks of cancer, and the main drivers of SF expression level and subcellular localization. In addition, we briefly introduce the application prospects of targeted therapeutic strategies, including small-molecule inhibitors, siRNAs and splice-switching oligonucleotides (SSOs), from three perspectives (drivers, SFs and pathological AS). Finally, we share our insights into the potential direction of research on SF-centric AS-related regulatory networks.

Poison Exon Splicing Regulates a Coordinated Network of SR Protein Expression during Differentiation and Tumorigenesis

The RNA isoform repertoire is regulated by splicing factor (SF) expression, and alterations in SF levels are associated with disease. SFs contain ultraconserved poison exon (PE) sequences that exhibit greater identity across species than nearby coding exons, but their physiological role and molecular regulation is incompletely understood. We show that PEs in serine-arginine-rich (SR) proteins, a family of 14 essential SFs, are differentially spliced during induced pluripotent stem cell (iPSC) differentiation and in tumors versus normal tissues. We uncover an extensive cross-regulatory network of SR proteins controlling their expression via alternative splicing coupled to nonsense-mediated decay. We define sequences that regulate PE inclusion and protein expression of the oncogenic SF TRA2β using an RNA-targeting CRISPR screen. We demonstrate location dependency of RS domain activity on regulation of TRA2β-PE using CRISPR artificial SFs. Finally, we develop splice-switching antisense oligonucleotides to reverse the increased skipping of TRA2β-PE detected in breast tumors, altering breast cancer cell viability, proliferation, and migration.

The SRSF3-MBNL1-Acin1 circuit constitutes an emerging axis to lessen DNA fragmentation in colorectal cancer via an alternative splicing mechanism

Altered alternative splicing (AS) events are considered pervasive causes that result in the development of carcinogenesis. Herein, we identified reprogrammed expression and splicing profiles of Muscle blind-like protein 1 (MBNL1) transcripts in tumorous tissues compared to those of adjacent normal tissues dissected from individual colorectal cancer (CRC) patients using whole-transcriptome analyses. MBNL1 transcript 8 (MBNL1₈) containing exons 5 and 7 was majorly generated by cancerous tissues and CRC-derived cell lines compared with those of the normal counterparts. Interplay between the exonic CA-rich element and upregulated SRSF3 facilitated the inclusion of MBNL1 exons 5 and 7, which encode a bipartite nuclear localization signal (NLS) and conformational NLS. Moreover, abundant SRSF3 interfered with the autoregulatory mechanism involved in utilization of MBNL1 exons 5 and 7, resulting in enrichment of the MBNL1₈ isoform in cultured CRC cell lines. Subsequently, an increase in the MBNL1₈ isoform drove a shift in the apoptotic chromatin condensation inducer in nucleus 1-S (Acin1-S) isoform to the Acin1-L isoform, leading to diminished DNA fragmentation in cultured CRC cells under oxidative stress. Taken together, SRSF3-MBNL1-Acin1 was demonstrated to constitute an emerging axis which is relevant to proapoptotic signatures and post-transcriptional events of CRC cells.

Tumorigenic de-differentiation: the alternative splicing way

The mechanism of acquisition of tumorigenic properties by somatic cells at the onset of cancer and later during relapse is a question of paramount importance in cancer biology. We have recently discovered a Muscleblind like-1 (MBNL1)-driven alternative-splicing mediated mechanism of tumorigenic de-differentiation that is associated with poor prognosis, relapse and metastasis in common cancer types.

The transcriptional repressor BCL11A promotes breast cancer metastasis

The phenotypes of each breast cancer subtype are defined by their transcriptomes. However, the transcription factors that regulate differential patterns of gene expression that contribute to specific disease outcomes are not well understood. Here, using gene silencing and overexpression approaches, RNA-Seq, and splicing analysis, we report that the transcription factor B-cell leukemia/lymphoma 11A (BCL11A) is highly expressed in triple-negative breast cancer (TNBC) and drives metastatic disease. Moreover, BCL11A promotes cancer cell invasion by suppressing the expression of muscleblind-like splicing regulator 1 (MBNL1), a splicing regulator that suppresses metastasis. This ultimately increases the levels of an alternatively spliced isoform of integrin-α6 (ITGA6), which is associated with worse patient outcomes. These results suggest that BCL11A sustains TNBC cell invasion and metastatic growth by repressing MBNL1-directed splicing of ITGA6 Our findings also indicate that BCL11A lies at the interface of transcription and splicing and promotes aggressive TNBC phenotypes.

Modelling the developmental spliceosomal craniofacial disorder Burn-McKeown syndrome using induced pluripotent stem cells

The craniofacial developmental disorder Burn-McKeown Syndrome (BMKS) is caused by biallelic variants in the pre-messenger RNA splicing factor gene TXNL4A/DIB1. The majority of affected individuals with BMKS have a 34 base pair deletion in the promoter region of one allele of TXNL4A combined with a loss-of-function variant on the other allele, resulting in reduced TXNL4A expression. However, it is unclear how reduced expression of this ubiquitously expressed spliceosome protein results in craniofacial defects during development. Here we reprogrammed peripheral mononuclear blood cells from a BMKS patient and her unaffected mother into induced pluripotent stem cells (iPSCs) and differentiated the iPSCs into induced neural crest cells (iNCCs), the key cell type required for correct craniofacial development. BMKS patient-derived iPSCs proliferated more slowly than both mother- and unrelated control-derived iPSCs, and RNA-Seq analysis revealed significant differences in gene expression and alternative splicing. Patient iPSCs displayed defective differentiation into iNCCs compared to maternal and unrelated control iPSCs, in particular a delay in undergoing an epithelial-to-mesenchymal transition (EMT). RNA-Seq analysis of differentiated iNCCs revealed widespread gene expression changes and mis-splicing in genes relevant to craniofacial and embryonic development that highlight a dampened response to WNT signalling, the key pathway activated during iNCC differentiation. Furthermore, we identified the mis-splicing of TCF7L2 exon 4, a key gene in the WNT pathway, as a potential cause of the downregulated WNT response in patient cells. Additionally, mis-spliced genes shared common sequence properties such as length, branch point to 3’ splice site (BPS-3’SS) distance and splice site strengths, suggesting that splicing of particular subsets of genes is particularly sensitive to changes in TXNL4A expression. Together, these data provide the first insight into how reduced TXNL4A expression in BMKS patients might compromise splicing and NCC function, resulting in defective craniofacial development in the embryo.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

A tumor-associated splice-isoform of MAP2K7 drives dedifferentiation in MBNL1-low cancers via JNK activation

Master splicing regulator MBNL1 shapes large transcriptomic changes that drive cellular differentiation during development. Here we demonstrate that MBNL1 is a suppressor of tumor dedifferentiation. We surveyed MBNL1 expression in matched tumor/normal pairs across The Cancer Genome Atlas and found that MBNL1 was down-regulated in several common cancers. Down-regulation of MBNL1 predicted poor overall survival in breast, lung, and stomach adenocarcinomas and increased relapse and distant metastasis in triple-negative breast cancer. Down-regulation of MBNL1 led to increased tumorigenic and stem/progenitor-like properties in vitro and in vivo. A discrete set of alternative splicing events (ASEs) are shared between MBNL1-low cancers and embryonic stem cells including a MAP2K7∆exon2 splice variant that leads to increased stem/progenitor-like properties via JNK activation. Accordingly, JNK inhibition is capable of reversing MAP2K7∆exon2-driven tumor dedifferentiation in MBNL1-low cancer cells. Our work elucidates an alternative-splicing mechanism that drives tumor dedifferentiation and identifies biomarkers that predict enhanced susceptibility to JNK inhibition.

Mammalian CBX7 isoforms p36 and p22 exhibit differential responses to serum, varying functions for proliferation, and distinct subcellular localization

CBX7 is a polycomb group protein, and despite being implicated in many diseases, its role in cell proliferation has been controversial: some groups described its pro-proliferative properties, but others illustrated its inhibitory effects on cell growth. To date, the reason for the divergent observations remains unknown. While several isoforms for CBX7 were reported, no studies investigated whether the divergent roles of CBX7 could be due to distinct functions of CBX7 isoforms. In this study, we newly identified mouse CBX7 transcript variant 1 (mCbx7v1), which is homologous to the human CBX7 gene (hCBX7v1) and is expressed in various mouse organs. We revealed that mCbx7v1 and hCBX7v1 encode a 36 kDa protein (p36^CBX7) whereas mCbx7 and hCBX7v3 encode a 22 kDa protein (p22^CBX7). This study further demonstrated that p36^CBX7 was localized to the nucleus and endogenously expressed in proliferating cells whereas p22^CBX7 was localized to the cytoplasm, induced by serum starvation in both human and mouse cells, and inhibited cell proliferation. Together, these data indicate that CBX7 isoforms are localized in different locations in a cell and play differing roles in cell proliferation. This varying function of CBX7 isoforms may help us understand the distinct function of CBX7 in various studies.

Muscleblind-like 2 controls the hypoxia response of cancer cells

Hypoxia is a hallmark of solid cancers, supporting proliferation, angiogenesis, and escape from apoptosis. There is still limited understanding of how cancer cells adapt to hypoxic conditions and survive. We analyzed transcriptome changes of human lung and breast cancer cells under chronic hypoxia. Hypoxia induced highly concordant changes in transcript abundance, but divergent splicing responses, underlining the cell type-specificity of alternative splicing programs. While RNA-binding proteins were predominantly reduced, hypoxia specifically induced muscleblind-like protein 2 (MBNL2). Strikingly, MBNL2 induction was critical for hypoxia adaptation by controlling the transcript abundance of hypoxia response genes, such as vascular endothelial growth factor A (VEGFA) MBNL2 depletion reduced the proliferation and migration of cancer cells, demonstrating an important role of MBNL2 as cancer driver. Hypoxia control is specific for MBNL2 and not shared by its paralog MBNL1. Thus, our study revealed MBNL2 as central mediator of cancer cell responses to hypoxia, regulating the expression and alternative splicing of hypoxia-induced genes.

A novel alternative splicing-based prediction model for uteri corpus endometrial carcinoma

Alternative splicing (AS) is crucial a mechanism by which the complexity of mammalian and viral proteom increased overwhelmingly. There lacks systematic and comprehensive analysis of the prognostic significance of AS profiling landscape for uteri corpus endometrial carcinoma (UCEC). In this study, univariate and multivariate Cox regression analyses were conducted to identify candidate survival-associated AS events curated from SpliceSeq for the construction of prognostic index (PI) models. A correlation network between splicing factor-related AS events and significant survival-associated AS events were constructed using Cytoscape 3.5. As consequence, 28281 AS events from 8137 genes were detected from 506 UCEC patients, including 2630 survival-associated AS events. Kaplan Meier survival analysis revealed that six of the seven PI models (AD, AP, AT, ME, RI and ALL) exhibited good performance in stratifying the prognosis of low risk and high risk group (P<0.05). Among the six PI models, PI-AT performed best with an area under curves (AUC) value of 0.758 from time-dependent receiver operating characteristic. Correlation network implicated potential regulatory mechanism of AS events in UCEC. PI models based on survival-associated AS events for UCEC in this study showed preferable prognosis-predicting ability and may be promising prognostic indicators for UCEC patients.

Summary: This is the first study to systematically investigate the prognostic value of AS in UCEC. Findings in the presents study supported the clinical potential of AS for UCEC and shed light on the potential AS-associated molecular basis of UCEC.

Integration of multi-omics data to mine cancer-related gene modules

The identification of cancer-related genes is a major research goal, with implications for determining the pathogenesis of cancer and identifying biomarkers for early diagnosis and treatment. In this study, by integrating multi-omics data, including gene expression, DNA copy number variation, DNA methylation, transcription factors, miRNA, and lncRNA data, we propose a method for mining cancer-related genes based on network models. First, using random forest-based feature selection method multi-omics data are integrated to identify key regulatory factors that affect gene expression, and then genome-wide regulatory networks are constructed. Next, by comparing the regulatory networks of key candidate genes in variant samples and non-variant samples, a differential expression regulatory network is generated. The differential network contains a collection of abnormal regulatory genes of key candidate genes. Then, by introducing the functional similarity as a distance metric for gene sets, a density-based clustering method is used to mine gene modules related to cancer. We applied this method to LUSC (lung squamous cell carcinoma) and mined cancer-related gene modules composed of 20 genes. GO function and KEGG pathway analyses indicated that the modules were closely related to cancer. A survival analysis was used to verify that the excavated gene modules can effectively distinguish between high- and low-risk groups. Overall, these results suggest that the proposed method can be used to identify cancer-related gene modules, providing a basis for the development of biomarkers for diagnosis and treatment.

Alternatively spliced MBNL1 isoforms exhibit differential influence on enhancing brown adipogenesis

Browning of white adipocytes (WAs) (also referred as beige cells) was demonstrated to execute thermogenesis by consuming stored lipids as do brown adipocytes (BAs), and this is highly related to metabolic homeostasis. Alternative splicing (AS) constitutes a pivotal mechanism for defining cellular fates and functional specifications. Nevertheless, the impacts of AS regulation on the browning of WAs have not been comprehensively investigated. In this study, we first identified the discriminative expression and splicing profiles of the muscleblind-like 1 (MBNL1) gene in postnatal brown adipose tissues (BATs) compared to those of embryonic BATs. A shift in the MBNL1^{+ex 5} isoform 7 (MBNL1₇) to MBNL1^{-ex 5} isoform 1 (MBNL1₁) was characterized throughout BAT development or during the in vitro browning of pre-WAs, 3T3-L1 cells. The interplay between MBNL1 and the exonic CCUG motif constitutes an autoregulatory mechanism for excluding MBNL1 exon 5. The simultaneous association of RNA-binding motif protein 4a (RBM4a) with exonic and intronic CU elements collaboratively mediates the skipping of MBNL1 exon 5. Overexpressing the MBNL1₁ isoform exhibited a more-prominent effect than that of the MBNL1₇ isoform on programming its own transcripts and beige cell-related splicing events in a CCUG motif-mediated manner. In addition to splicing regulation, overexpression of the MBNL1₁ and MBNL1₇ isoforms differentially enhanced beige adipogenic signatures of 3T3-L1 cells. Our findings demonstrated that MBNL1 constitutes an emerging and autoregulatory mechanism involved in development of beige cells.

Quantitative screening of the effects of hyper-osmotic stress on cancer cells cultured in 2- or 3-dimensional settings

The maintenance of precise cell volume is critical for cell survival. Changes in extracellular osmolarity affect cell volume and may impact various cellular processes such as mitosis, mitochondrial functions, DNA repair as well as cell migration and proliferation. Much of what we know about the mechanisms of cell osmoregulation comes from in vitro two-dimensional (2D) assays that are less physiologically relevant than three-dimensional (3D) in vitro or in vivo settings. Here, we developed a microfluidic model to study the impact of hyper-osmotic stress on the migration, proliferation and ion channel/transporter expression changes of three metastatic cell lines (MDA-MB-231, A549, T24) in 2D versus 3D environments. We observed a global decrease in cell migration and proliferation upon hyper-osmotic stress treatment, with similar responses between 2D and 3D conditions. Specific ion channels/aquaporins are over-expressed in metastatic cells and play a central role during osmo-regulation. Therefore, the effects of hyper-osmotic stress on two transporters, aquaporin 5 (AQP5) and the transient receptor potential cation channel (TRPV4), was investigated. While hyper-osmotic stress had no major impact on the transporters of cells cultured in 2D, cells embedded in collagen gel (3D) decreased their AQP5 expression and exhibited a reduction in intra-cellular translocation of TRPV4. Furthermore, cell dispersion from T24 aggregates embedded in 3D collagen gel decreased with higher levels of hyper-osmotic stress. In conclusion, this study provides evidence on the impact of hyper-osmotic stress on various aspects of metastatic cell progression and highlights the importance of having a 3D cell culture platform in investigating molecular players involved in cancer cell migration.

A Computational Analysis of Alternative Splicing across Mammalian Tissues Reveals Circadian and Ultradian Rhythms in Splicing Events

Mounting evidence points to a role of the circadian clock in the temporal regulation of post-transcriptional processes in mammals, including alternative splicing (AS). In this study, we carried out a computational analysis of circadian and ultradian rhythms on the transcriptome level to characterise the landscape of rhythmic AS events in published datasets covering 76 tissues from mouse and olive baboon. Splicing-related genes with 24-h rhythmic expression patterns showed a bimodal distribution of peak phases across tissues and species, indicating that they might be controlled by the circadian clock. On the output level, we identified putative oscillating AS events in murine microarray data and pairs of differentially rhythmic splice isoforms of the same gene in baboon RNA-seq data that peaked at opposing times of the day and included oncogenes and tumour suppressors. We further explored these findings using a new circadian RNA-seq dataset of human colorectal cancer cell lines. Rhythmic isoform expression patterns differed between the primary tumour and the metastatic cell line and were associated with cancer-related biological processes, indicating a functional role of rhythmic AS that might be implicated in tumour progression. Our data shows that rhythmic AS events are widespread across mammalian tissues and might contribute to a temporal diversification of the proteome.

Short Tandem Repeat Expansions and RNA-Mediated Pathogenesis in Myotonic Dystrophy

Short tandem repeat (STR) or microsatellite, expansions underlie more than 50 hereditary neurological, neuromuscular and other diseases, including myotonic dystrophy types 1 (DM1) and 2 (DM2). Current disease models for DM1 and DM2 propose a common pathomechanism, whereby the transcription of mutant DMPK (DM1) and CNBP (DM2) genes results in the synthesis of CUG and CCUG repeat expansion (CUG^exp, CCUG^exp) RNAs, respectively. These CUG^exp and CCUG^exp RNAs are toxic since they promote the assembly of ribonucleoprotein (RNP) complexes or RNA foci, leading to sequestration of Muscleblind-like (MBNL) proteins in the nucleus and global dysregulation of the processing, localization and stability of MBNL target RNAs. STR expansion RNAs also form phase-separated gel-like droplets both in vitro and in transiently transfected cells, implicating RNA-RNA multivalent interactions as drivers of RNA foci formation. Importantly, the nucleation and growth of these nuclear foci and transcript misprocessing are reversible processes and thus amenable to therapeutic intervention. In this review, we provide an overview of potential DM1 and DM2 pathomechanisms, followed by a discussion of MBNL functions in RNA processing and how multivalent interactions between expanded STR RNAs and RNA-binding proteins (RBPs) promote RNA foci assembly.

Identification of RBPMS as a mammalian smooth muscle master splicing regulator via proximity of its gene with super-enhancers

Alternative splicing (AS) programs are primarily controlled by regulatory RNA-binding proteins (RBPs). It has been proposed that a small number of master splicing regulators might control cell-specific splicing networks and that these RBPs could be identified by proximity of their genes to transcriptional super-enhancers. Using this approach we identified RBPMS as a critical splicing regulator in differentiated vascular smooth muscle cells (SMCs). RBPMS is highly down-regulated during phenotypic switching of SMCs from a contractile to a motile and proliferative phenotype and is responsible for 20% of the AS changes during this transition. RBPMS directly regulates AS of numerous components of the actin cytoskeleton and focal adhesion machineries whose activity is critical for SMC function in both phenotypes. RBPMS also regulates splicing of other splicing, post-transcriptional and transcription regulators including the key SMC transcription factor Myocardin, thereby matching many of the criteria of a master regulator of AS in SMCs.