Human Tra2 proteins jointly control a CHEK1 splicing switch among alternative and constitutive target exons

Incomplete paralog compensation generates selective dependency on TRA2A in cancer

Paralogs often exhibit functional redundancy, allowing them to effectively compensate for each other's loss. However, this buffering mechanism is frequently disrupted in cancer, exposing unique paralog-specific vulnerabilities. Here, we identify a selective dependency on the splicing factor TRA2A. We find that TRA2A and its paralog TRA2B are synthetic lethal partners that function as widespread and largely redundant activators of both alternative and constitutive splicing. While loss of TRA2A alone is typically neutral due to compensation by TRA2B, we discover that a subset of cancer cell lines are highly TRA2A-dependent. Upon TRA2A depletion, these cell lines exhibit a lack of paralog buffering specifically on shared splicing targets, leading to defects in mitosis and cell death. Notably, TRA2B overexpression rescues both the aberrant splicing and lethality associated with TRA2A loss, indicating that paralog compensation is dosage-sensitive. Together, these findings reveal a complex dosage-dependent relationship between paralogous splicing factors, and highlight how dysfunctional paralog buffering can create a selective dependency in cancer.

Antisense oligonucleotide-mediated TRA2β poison exon inclusion induces the expression of a lncRNA with anti-tumor effects

Upregulated expression of the oncogenic splicing factor TRA2β occurs in human tumors partly through decreased inclusion of its autoregulatory non-coding poison exon (PE). Here, we reveal that low TRA2β-PE inclusion negatively impacts patient survival across several tumor types. We demonstrate the ability of splice-switching antisense oligonucleotides (ASOs) to promote TRA2β-PE inclusion and lower TRA2β protein levels in pre-clinical cancer models. TRA2β-PE-targeting ASOs induce anti-cancer phenotypes and widespread transcriptomic alterations with functional impact on RNA processing, mTOR, and p53 signaling pathways. Surprisingly, the effect of TRA2β-PE-targeting ASOs on cell viability are not phenocopied by TRA2β knockdown. Mechanistically, we find that the ASO functions by both decreasing TRA2β protein and inducing the expression of TRA2β-PE-containing transcripts that act as long non-coding RNAs to sequester nuclear proteins. Finally, TRA2β-PE-targeting ASOs are toxic to preclinical 3D organoid and in vivo patient-derived xenograft models. Together, we demonstrate that TRA2β-PE acts both as a regulator of protein expression and a long-noncoding RNA to control cancer cell growth. Drugging oncogenic splicing factors using PE-targeting ASOs is a promising therapeutic strategy.

An ultra-conserved poison exon in the Tra2b gene encoding a splicing activator is essential for male fertility and meiotic cell division

The cellular concentrations of splicing factors (SFs) are critical for controlling alternative splicing. Most serine and arginine-enriched (SR) protein SFs regulate their own concentration via a homeostatic feedback mechanism that involves regulation of inclusion of non-coding 'poison exons' (PEs) that target transcripts for nonsense-mediated decay. The importance of SR protein PE splicing during animal development is largely unknown despite PE ultra-conservation across animal genomes. To address this, we used mouse genetics to disrupt an ultra-conserved PE in the Tra2b gene encoding the SR protein Tra2β. Focussing on germ cell development, we found that Tra2b PE deletion causes azoospermia due to catastrophic cell death during meiotic prophase. Failure to proceed through meiosis was associated with increased Tra2b expression sufficient to drive aberrant Tra2β protein hyper-responsive splice patterns. Although critical for meiotic prophase, Tra2b PE deletion spared earlier mitotically active germ cells, even though these still required Tra2b gene function. Our data indicate that PE splicing control prevents the accumulation of toxic levels of Tra2β protein that are incompatible with meiotic prophase. This unexpected connection with male fertility helps explain Tra2b PE ultra-conservation and indicates the importance of evaluating PE function in animal models.

Differential quantification of alternative splicing events on spliced pangenome graphs

Pangenomes are becoming a powerful framework to perform many bioinformatics analyses taking into account the genetic variability of a population, thus reducing the bias introduced by a single reference genome. With the wider diffusion of pangenomes, integrating genetic variability with transcriptome diversity is becoming a natural extension that demands specific methods for its exploration. In this work, we extend the notion of spliced pangenomes to that of annotated spliced pangenomes; this allows us to introduce a formal definition of Alternative Splicing (AS) events on a graph structure. To investigate the usage of graph pangenomes for the quantification of AS events across conditions, we developed pantas, the first pangenomic method for the detection and differential analysis of AS events from short RNA-Seq reads. A comparison with state-of-the-art linear reference-based approaches proves that pantas achieves competitive accuracy, making spliced pangenomes effective for conducting AS events quantification and opening future directions for the analysis of population-based transcriptomes.

Autoregulated splicing of TRA2β programs T cell fate in response to antigen-receptor stimulation

T cell receptor (TCR) sensitivity to peptide-major histocompatibility complex (MHC) dictates T cell fate. Canonical models of TCR sensitivity cannot be fully explained by transcriptional regulation. In this work, we identify a posttranscriptional regulatory mechanism of TCR sensitivity that guides alternative splicing of TCR signaling transcripts through an evolutionarily ultraconserved poison exon (PE) in the RNA-binding protein (RBP) TRA2β in mouse and human. TRA2β-PE splicing, seen during cancer and infection, was required for TCR-induced effector T cell expansion and function. Tra2β-PE skipping enhanced T cell response to antigen by increasing TCR sensitivity. As antigen levels decreased, Tra2β-PE reinclusion allowed T cell survival. Finally, we found that TRA2β-PE was first included in the genome of jawed vertebrates that were capable of TCR gene rearrangements. We propose that TRA2β-PE splicing acts as a gatekeeper of TCR sensitivity to shape T cell fate.

SpliceWiz: interactive analysis and visualization of alternative splicing in R

Alternative splicing (AS) is a crucial mechanism for regulating gene expression and isoform diversity in eukaryotes. However, the analysis and visualization of AS events from RNA sequencing data remains challenging. Most tools require a certain level of computer literacy and the available means of visualizing AS events, such as coverage and sashimi plots, have limitations and can be misleading. To address these issues, we present SpliceWiz, an R package with an interactive Shiny interface that allows easy and efficient AS analysis and visualization at scale. A novel normalization algorithm is implemented to aggregate splicing levels within sample groups, thereby allowing group differences in splicing levels to be accurately visualized. The tool also offers downstream gene ontology enrichment analysis, highlighting ASEs belonging to functional pathways of interest. SpliceWiz is optimized for speed and efficiency and introduces a new file format for coverage data storage that is more efficient than BigWig. Alignment files are processed orders of magnitude faster than other R-based AS analysis tools and on par with command-line tools. Overall, SpliceWiz streamlines AS analysis, enabling reliable identification of functionally relevant AS events for further characterization. SpliceWiz is a Bioconductor package and is also available on GitHub (https://github.com/alexchwong/SpliceWiz).

Lnc-ZEB2-19 Inhibits the Progression and Lenvatinib Resistance of Hepatocellular Carcinoma by Attenuating the NF-κB Signaling Pathway through the TRA2A/RSPH14 Axis

Long non-coding RNAs have been reported to play a crucial role in tumor progression in hepatocellular carcinoma (HCC). Lnc-ZEB2-19 has been validated to be deficiently expressed in HCC. However, the capabilities and underlying mechanisms of lnc-ZEB2-19 remain uncertain. In this study, we verified that the downregulation of lnc-ZEB2-19 was prevalent in HCC and significantly correlated with the unfavorable prognosis. Further in vitro and in vivo verified that lnc-ZEB2-19 notably inhibited the proliferation, metastasis, stemness, and lenvatinib resistance (LR) of HCC cells. Mechanistically, lnc-ZEB2-19 inhibited HCC progression and LR by specifically binding to transformer 2α (TRA2A) and promoting its degradation, which resulted in the instability of RSPH14 mRNA, leading to the downregulation of Rela(p65) and p-Rela(p-p65). Furthermore, rescue assays showed that silencing RSPH14 partially restrained the effect of knockdown expression of lnc-ZEB2-19 on HCC cell metastatic ability and stemness. The findings describe a novel regulatory axis, lnc-ZEB2-19/TRA2A/RSPH14, downregulating the nuclear factor kappa B to inhibit HCC progression and LR.

TRA2: The dominant power of alternative splicing in tumors

The dysregulation of alternative splicing (AS) is frequently found in cancer and considered as key markers for cancer progression and therapy. Transformer 2 (TRA2), a nuclear RNA binding protein, consists of transformer 2 alpha homolog (TRA2A) and transformer 2 beta homolog (TRA2B), and plays a role in the regulation of pre-mRNA splicing. Growing evidence has been provided that TRA2A and TRA2B are dysregulated in several types of tumors, and participate in the regulation of proliferation, migration, invasion, and chemotherapy resistance in cancer cells through alteration of AS of cancer-related genes. In this review, we highlight the role of TRA2 in tumorigenesis and metastasis, and discuss potential molecular mechanisms how TRA2 influences tumorigenesis and metastasis via controlling AS of pre-mRNA. We propose that TRA2Ais a novel biomarker and therapeutic target for cancer progression and therapy.

DNA Repair and Therapeutic Strategies in Cancer Stem Cells

First-line cancer treatments successfully eradicate the differentiated tumour mass but are comparatively ineffective against cancer stem cells (CSCs), a self-renewing subpopulation thought to be responsible for tumour initiation, metastasis, heterogeneity, and recurrence. CSCs are thus presented as the principal target for elimination during cancer treatment. However, CSCs are challenging to drug target because of numerous intrinsic and extrinsic mechanisms of drug resistance. One such mechanism that remains relatively understudied is the DNA damage response (DDR). CSCs are presumed to possess properties that enable enhanced DNA repair efficiency relative to their highly proliferative bulk progeny, facilitating improved repair of double-strand breaks induced by radiotherapy and most chemotherapeutics. This can occur through multiple mechanisms, including increased expression and splicing fidelity of DNA repair genes, robust activation of cell cycle checkpoints, and elevated homologous recombination-mediated DNA repair. Herein, we summarise the current knowledge concerning improved genome integrity in non-transformed stem cells and CSCs, discuss therapeutic opportunities within the DDR for re-sensitising CSCs to genotoxic stressors, and consider the challenges posed regarding unbiased identification of novel DDR-directed strategies in CSCs. A better understanding of the DDR mediating chemo/radioresistance mechanisms in CSCs could lead to novel therapeutic approaches, thereby enhancing treatment efficacy in cancer patients.

Clustered variants in the 5' coding region of TRA2B cause a distinctive neurodevelopmental syndrome

PURPOSE

Transformer2 proteins (Tra2α and Tra2β) control splicing patterns in human cells, and no human phenotypes have been associated with germline variants in these genes. The aim of this work was to associate germline variants in the TRA2B gene to a novel neurodevelopmental disorder.

METHODS

A total of 12 individuals from 11 unrelated families who harbored predicted loss-of-function monoallelic variants, mostly de novo, were recruited. RNA sequencing and western blot analyses of Tra2β-1 and Tra2β-3 isoforms from patient-derived cells were performed. Tra2β1-GFP, Tra2β3-GFP and CHEK1 exon 3 plasmids were transfected into HEK-293 cells.

RESULTS

All variants clustered in the 5' part of TRA2B, upstream of an alternative translation start site responsible for the expression of the noncanonical Tra2β-3 isoform. All affected individuals presented intellectual disability and/or developmental delay, frequently associated with infantile spasms, microcephaly, brain anomalies, autism spectrum disorder, feeding difficulties, and short stature. Experimental studies showed that these variants decreased the expression of the canonical Tra2β-1 isoform, whereas they increased the expression of the Tra2β-3 isoform, which is shorter and lacks the N-terminal RS1 domain. Increased expression of Tra2β-3-GFP were shown to interfere with the incorporation of CHEK1 exon 3 into its mature transcript, normally incorporated by Tra2β-1.

CONCLUSION

Predicted loss-of-function variants clustered in the 5' portion of TRA2B cause a new neurodevelopmental syndrome through an apparently dominant negative disease mechanism involving the use of an alternative translation start site and the overexpression of a shorter, repressive Tra2β protein.

Differential Degradation of TRA2A and PYCR2 Mediated by Ubiquitin E3 Ligase E4B

E4B belongs to the U-box E3 ligase family and functions as either an E3 or an E4 enzyme in protein ubiquitination. Transformer2A (TRA2A) and Pyrroline-5-carboxylate reductase 2 (PYCR2) are related to cancer development and are overexpressed in many cancer cells. The degradation of TRA2A and PYCR2 mediated by the ubiquitin-proteasome system (UPS) has not been reported. This study validated that E4B could ubiquitinate TRA2A and PYCR2 as an E3 ligase both in vitro and in the HEK293 cells. E4B mediated the degradation by forming K11- and K48- linked polyubiquitin chains on TRA2A and PYCR2, respectively. E4B regulated the alternative splicing function of TRA2A and affected RSRC2 transcription in the HEK293 cells. Although E4B is highly expressed, it hardly degrades TRA2A and PYCR2 in hepatocellular carcinoma (HCC) cells, suggesting other mechanisms exist for degradation of TRA2A and PYCR2 in the HCC cells. We finally reported that E4B interacted with substrates via its variable region.

Tra2beta-Dependent Regulation of RIO Kinase 3 Splicing During Rift Valley Fever Virus Infection Underscores the Links Between Alternative Splicing and Innate Antiviral Immunity

Rift Valley fever virus (RVFV) is an emerging pathogen that has potential to cause severe disease in humans and domestic livestock. Propagation of RVFV strain MP-12 is negatively impacted by the actions of RIOK3, a protein involved in the cellular immune response to viral infection. During RVFV infection, RIOK3 mRNA is alternatively spliced to produce an isoform that correlates with the inhibition of interferon β signaling. Here, we identify splicing factor TRA2-β (also known as TRA2beta and hTRA2-β) as a key regulator governing the relative abundance of RIOK3 splicing isoforms. Using RT-PCR and minigenes, we determined that TRA2-β interaction with RIOK3 pre-mRNA was necessary for constitutive splicing of RIOK3 mRNA, and conversely, lack of TRA2-β engagement led to increased alternative splicing. Expression of TRA2-β was found to be necessary for RIOK3's antiviral effect against RVFV. Intriguingly, TRA2-β mRNA is also alternatively spliced during RVFV infection, leading to a decrease in cellular TRA2-β protein levels. These results suggest that splicing modulation serves as an immune evasion strategy by RVFV and/or is a cellular mechanism to prevent excessive immune response. Furthermore, the results suggest that TRA2-β can act as a key regulator of additional steps of the innate immune response to viral infection.

Shark: fishing relevant reads in an RNA-Seq sample

Motivation

Recent advances in high-throughput RNA-Seq technologies allow to produce massive datasets. When a study focuses only on a handful of genes, most reads are not relevant and degrade the performance of the tools used to analyze the data. Removing irrelevant reads from the input dataset leads to improved efficiency without compromising the results of the study.

Results

We introduce a novel computational problem, called gene assignment and we propose an efficient alignment-free approach to solve it. Given an RNA-Seq sample and a panel of genes, a gene assignment consists in extracting from the sample, the reads that most probably were sequenced from those genes. The problem becomes more complicated when the sample exhibits evidence of novel alternative splicing events. We implemented our approach in a tool called Shark and assessed its effectiveness in speeding up differential splicing analysis pipelines. This evaluation shows that Shark is able to significantly improve the performance of RNA-Seq analysis tools without having any impact on the final results.

Availability and implementation

The tool is distributed as a stand-alone module and the software is freely available at https://github.com/AlgoLab/shark.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Identification of protein-protected mRNA fragments and structured excised intron RNAs in human plasma by TGIRT-seq peak calling

Human plasma contains > 40,000 different coding and non-coding RNAs that are potential biomarkers for human diseases. Here, we used thermostable group II intron reverse transcriptase sequencing (TGIRT-seq) combined with peak calling to simultaneously profile all RNA biotypes in apheresis-prepared human plasma pooled from healthy individuals. Extending previous TGIRT-seq analysis, we found that human plasma contains largely fragmented mRNAs from > 19,000 protein-coding genes, abundant full-length, mature tRNAs and other structured small non-coding RNAs, and less abundant tRNA fragments and mature and pre-miRNAs. Many of the mRNA fragments identified by peak calling correspond to annotated protein-binding sites and/or have stable predicted secondary structures that could afford protection from plasma nucleases. Peak calling also identified novel repeat RNAs, miRNA-sized RNAs, and putatively structured intron RNAs of potential biological, evolutionary, and biomarker significance, including a family of full-length excised intron RNAs, subsets of which correspond to mirtron pre-miRNAs or agotrons.

Arsenite Exposure Displaces Zinc from ZRANB2 Leading to Altered Splicing

Exposure to arsenic, a class I carcinogen, affects 200 million people globally. Skin is the major target organ, but the molecular etiology of arsenic-induced skin carcinogenesis remains unclear. Arsenite (As³⁺)-induced disruption of alternative splicing could be involved, but the mechanism is unknown. Zinc finger proteins play key roles in alternative splicing. As³⁺ can displace zinc (Zn²⁺) from C3H1 and C4 zinc finger motifs (zfm's), affecting protein function. ZRANB2, an alternative splicing regulator with two C4 zfm's integral to its structure and splicing function, was chosen as a candidate for this study. We hypothesized that As³⁺ could displace Zn²⁺ from ZRANB2, altering its structure, expression, and splicing function. As³⁺/Zn²⁺ binding and mutual displacement experiments were performed with synthetic apo-peptides corresponding to each ZRANB2 zfm, employing a combination of intrinsic fluorescence, ultraviolet spectrophotometry, zinc colorimetric assay, and liquid chromatography-tandem mass spectrometry. ZRANB2 expression in HaCaT cells acutely exposed to As³⁺ (0 or 5 μM, 0-72 h; or 0-5 μM, 6 h) was examined by RT-qPCR and immunoblotting. ZRANB2-dependent splicing of TRA2B mRNA, a known ZRANB2 target, was monitored by reverse transcription-polymerase chain reaction. As³⁺ bound to, as well as displaced Zn²⁺ from, each zfm. Also, Zn²⁺ displaced As³⁺ from As³⁺-bound zfm's acutely, albeit transiently. As³⁺ exposure induced ZRANB2 protein expression between 3 and 24 h and at all exposures tested but not ZRANB2 mRNA expression. ZRANB2-directed TRA2B splicing was impaired between 3 and 24 h post-exposure. Furthermore, ZRANB2 splicing function was also compromised at all As³⁺ exposures, starting at 100 nm. We conclude that As³⁺ exposure displaces Zn²⁺ from ZRANB2 zfm's, changing its structure and compromising splicing of its targets, and increases ZRANB2 protein expression as a homeostatic response both at environmental/toxicological exposures and therapeutically relevant doses.

Clinical and functional significance of CHK1-S, an alternatively spliced isoform of the CHK1 gene, in hepatocellular carcinoma

Alternative splicing plays critical roles in many disease processes and splicing dysregulation is a hallmark of cancer. The different splicing isoforms may have significantly different effects on the malignant progression of cancer. Checkpoint kinase 1 (CHK1) is a serine/threonine kinase and regulates DNA damage response. In this study, we measured the expression of an alternative CHK1 transcript (CHK1-S, excluded exon 3) in hepatocellular carcinoma (HCC) tissues. Our results showed that CHK1-S was significantly upregulated in HCC tissues compared with paired adjacent noncancerous hepatic tissues. The levels of full-length CHK1(CHK1-L), CHK1-S and the ratio of CHK1-S/L in tumor tissue were associated with relapse free survival (RFS) of postoperative HCC patients, respectively, but not the levels of CHK1-L, CHK1-S and the ratio of CHK1-S/L in adjacent normal tissue. To further demonstrate the role of CHK1-S in HCC, CCK-8 assays, EdU incorporation assays and colony formation assays were used. The results showed that overexpression of CHK1-S significantly accelerated HCC cell proliferation, compared with CHK1-L. In addition, we found that serine-arginine protein kinase 1 (SRPK1), as an upstream regulator kinase of splicing factor, could upregulate the expression of CHK1-S and its expression level was significantly higher in HCC tumors than the paired normal tissues and was associated with the levels of CHK1-S (P=0.016). In conclusion, our study demonstrated that CHK1-S, acts as an oncogene, which was upregulated and associated with RFS in HCC patients. SRPK1 may mediate its mRNA splicing in HCC. All these data indicated that the expression of CHK1-S would have potential prognostic values and splicing kinase SRPK1 might be developed as therapeutic target in HCC.

BANDITS: Bayesian differential splicing accounting for sample-to-sample variability and mapping uncertainty

Alternative splicing is a biological process during gene expression that allows a single gene to code for multiple proteins. However, splicing patterns can be altered in some conditions or diseases. Here, we present BANDITS, a R/Bioconductor package to perform differential splicing, at both gene and transcript level, based on RNA-seq data. BANDITS uses a Bayesian hierarchical structure to explicitly model the variability between samples and treats the transcript allocation of reads as latent variables. We perform an extensive benchmark across both simulated and experimental RNA-seq datasets, where BANDITS has extremely favourable performance with respect to the competitors considered.

Characterizing the interplay between gene nucleotide composition bias and splicing

BACKGROUND

Nucleotide composition bias plays an important role in the 1D and 3D organization of the human genome. Here, we investigate the potential interplay between nucleotide composition bias and the regulation of exon recognition during splicing.

RESULTS

By analyzing dozens of RNA-seq datasets, we identify two groups of splicing factors that activate either about 3200 GC-rich exons or about 4000 AT-rich exons. We show that splicing factor-dependent GC-rich exons have predicted RNA secondary structures at 5' ss and are dependent on U1 snRNP-associated proteins. In contrast, splicing factor-dependent AT-rich exons have a large number of decoy branch points, SF1- or U2AF2-binding sites and are dependent on U2 snRNP-associated proteins. Nucleotide composition bias also influences local chromatin organization, with consequences for exon recognition during splicing. Interestingly, the GC content of exons correlates with that of their hosting genes, isochores, and topologically associated domains.

CONCLUSIONS

We propose that regional nucleotide composition bias over several dozens of kilobase pairs leaves a local footprint at the exon level and induces constraints during splicing that can be alleviated by local chromatin organization at the DNA level and recruitment of specific splicing factors at the RNA level. Therefore, nucleotide composition bias establishes a direct link between genome organization and local regulatory processes, like alternative splicing.

Alternative splicing in lung cancer

Lung cancer has the highest mortality rate of all cancers worldwide. Lung cancer is a very heterogeneous disease that is often diagnosed at later stages which have a poor prognosis. Aberrant alternative splicing patterns found in lung cancer contribute to important cell functions. These include changes in splicing for the BCL2L1, MDM2, MDM4, NUMB and MET genes during lung tumourigenesis, to affect pathways involved in apoptosis, cell proliferation and cellular cohesion. Global analyses of RNASeq datasets suggest there may be many more potentially influential aberrant splicing events that need to be investigated in lung cancer. Changes in expression of the splicing factors that regulate alternative splicing events have also been identified in lung cancer. Of these, changes in expression of QKI, RBM4, RBM5, RBM6, RBM10 and SRSF1 proteins regulate many of the most frequently referenced aberrant splicing events in lung cancer. The expanding list of genes known to be aberrantly spliced in lung cancer along with the altered expression of splicing factors that regulate them are providing new clues as to how lung cancer develops, and how these events can be exploited for better treatment. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

RBMX family proteins connect the fields of nuclear RNA processing, disease and sex chromosome biology

RBMX is a ubiquitously expressed nuclear RNA binding protein that is encoded by a gene on the X chromosome. RBMX belongs to a small protein family with additional members encoded by paralogs on the mammalian Y chromosome and other chromosomes. These RNA binding proteins are important for normal development, and also implicated in cancer and viral infection. At the molecular level RBMX family proteins contribute to splicing control, transcription and genome integrity. Establishing what endogenous genes and pathways are controlled by RBMX and its paralogs will have important implications for understanding chromosome biology, DNA repair and mammalian development. Here we review what is known about this family of RNA binding proteins, and identify important current questions about their functions.

ASGAL: aligning RNA-Seq data to a splicing graph to detect novel alternative splicing events

Background

While the reconstruction of transcripts from a sample of RNA-Seq data is a computationally expensive and complicated task, the detection of splicing events from RNA-Seq data and a gene annotation is computationally feasible. This latter task, which is adequate for many transcriptome analyses, is usually achieved by aligning the reads to a reference genome, followed by comparing the alignments with a gene annotation, often implicitly represented by a graph: the splicing graph.

Results

We present ASGAL (Alternative Splicing Graph ALigner): a tool for mapping RNA-Seq data to the splicing graph, with the specific goal of detecting novel splicing events, involving either annotated or unannotated splice sites. ASGAL takes as input the annotated transcripts of a gene and a RNA-Seq sample, and computes (1) the spliced alignments of each read in input, and (2) a list of novel events with respect to the gene annotation.

Conclusions

An experimental analysis shows that ASGAL allows to enrich the annotation with novel alternative splicing events even when genes in an experiment express at most one isoform. Compared with other tools which use the spliced alignment of reads against a reference genome for differential analysis, ASGAL better predicts events that use splice sites which are novel with respect to a splicing graph, showing a higher accuracy. To the best of our knowledge, ASGAL is the first tool that detects novel alternative splicing events by directly aligning reads to a splicing graph.

Availability

Source code, documentation, and data are available for download at http://asgal.algolab.eu.

TRA2A promotes proliferation, migration, invasion and epithelial mesenchymal transition of glioma cells

TRA2A, Transformer2A proteins, plays important roles in paclitaxel resistance and progression of breast cancer. However, whether TRA2A was involved in the progression of glioma remains to be elucidated. In this study, our results demonstrated that the expression of TRA2A was higher in the glioma tissue than that of normal tissue. Overexpression of TRA2A in glioma SHG44 cell lines promoted the tumor cells proliferation, migration, invasion and epithelial mesenchymal transition (EMT), while, knockdown of TRA2A showed the opposite effect. Thus, our findings provide new insights into the role of TRA2A in the progression of glioma, and implicate the potential application of TRA2A in glioma therapy.

SUPPA2: fast, accurate, and uncertainty-aware differential splicing analysis across multiple conditions

Despite the many approaches to study differential splicing from RNA-seq, many challenges remain unsolved, including computing capacity and sequencing depth requirements. Here we present SUPPA2, a new method that addresses these challenges, and enables streamlined analysis across multiple conditions taking into account biological variability. Using experimental and simulated data, we show that SUPPA2 achieves higher accuracy compared to other methods, especially at low sequencing depth and short read length. We use SUPPA2 to identify novel Transformer2-regulated exons, novel microexons induced during differentiation of bipolar neurons, and novel intron retention events during erythroblast differentiation.

SRSF10 Connects DNA Damage to the Alternative Splicing of Transcripts Encoding Apoptosis, Cell-Cycle Control, and DNA Repair Factors

RNA binding proteins and signaling components control the production of pro-death and pro-survival splice variants of Bcl-x. DNA damage promoted by oxaliplatin increases the level of pro-apoptotic Bcl-xS in an ATM/CHK2-dependent manner, but how this shift is enforced is not known. Here, we show that in normally growing cells, when the 5′ splice site of Bcl-xS is largely repressed, SRSF10 partially relieves repression and interacts with repressor hnRNP K and stimulatory hnRNP F/H proteins. Oxaliplatin abrogates the interaction of SRSF10 with hnRNP F/H and decreases the association of SRSF10 and hnRNP K with the Bcl-x pre-mRNA. Dephosphorylation of SRSF10 is linked with these changes. A broader analysis reveals that DNA damage co-opts SRSF10 to control splicing decisions in transcripts encoding components involved in DNA repair, cell-cycle control, and apoptosis. DNA damage therefore alters the interactions between splicing regulators to elicit a splicing response that determines cell fate.

Proteomic analysis of polyribosomes identifies splicing factors as potential regulators of translation during mitosis

Precise regulation of mRNA translation is critical for proper cell division, but little is known about the factors that mediate it. To identify mRNA-binding proteins that regulate translation during mitosis, we analyzed the composition of polysomes from interphase and mitotic cells using unbiased quantitative mass-spectrometry (LC–MS/MS). We found that mitotic polysomes are enriched with a subset of proteins involved in RNA processing, including alternative splicing and RNA export. To demonstrate that these may indeed be regulators of translation, we focused on heterogeneous nuclear ribonucleoprotein C (hnRNP C) as a test case and confirmed that it is recruited to elongating ribosomes during mitosis. Then, using a combination of pulsed SILAC, metabolic labeling and ribosome profiling, we showed that knockdown of hnRNP C affects both global and transcript-specific translation rates and found that hnRNP C is specifically important for translation of mRNAs that encode ribosomal proteins and translation factors. Taken together, our results demonstrate how proteomic analysis of polysomes can provide insight into translation regulation under various cellular conditions of interest and suggest that hnRNP C facilitates production of translation machinery components during mitosis to provide daughter cells with the ability to efficiently synthesize proteins as they enter G1 phase.

RNA splicing and splicing regulator changes in prostate cancer pathology

Changes in mRNA splice patterns have been associated with key pathological mechanisms in prostate cancer progression. The androgen receptor (abbreviated AR) transcription factor is a major driver of prostate cancer pathology and activated by androgen steroid hormones. Selection of alternative promoters by the activated AR can critically alter gene function by switching mRNA isoform production, including creating a pro-oncogenic isoform of the normally tumour suppressor gene TSC2. A number of androgen-regulated genes generate alternatively spliced mRNA isoforms, including a prostate-specific splice isoform of ST6GALNAC1 mRNA. ST6GALNAC1 encodes a sialyltransferase that catalyses the synthesis of the cancer-associated sTn antigen important for cell mobility. Genetic rearrangements occurring early in prostate cancer development place ERG oncogene expression under the control of the androgen-regulated TMPRSS2 promoter to hijack cell behaviour. This TMPRSS2-ERG fusion gene shows different patterns of alternative splicing in invasive versus localised prostate cancer. Alternative AR mRNA isoforms play a key role in the generation of prostate cancer drug resistance, by providing a mechanism through which prostate cancer cells can grow in limited serum androgen concentrations. A number of splicing regulator proteins change expression patterns in prostate cancer and may help drive key stages of disease progression. Up-regulation of SRRM4 establishes neuronal splicing patterns in neuroendocrine prostate cancer. The splicing regulators Sam68 and Tra2β increase expression in prostate cancer. The SR protein kinase SRPK1 that modulates the activity of SR proteins is up-regulated in prostate cancer and has already given encouraging results as a potential therapeutic target in mouse models.

Modulation of the splicing regulatory function of SRSF10 by a novel compound that impairs HIV-1 replication

We recently identified the 4-pyridinone-benzisothiazole carboxamide compound 1C8 as displaying strong anti-HIV-1 potency against a variety of clinical strains in vitro. Here we show that 1C8 decreases the expression of HIV-1 and alters splicing events involved in the production of HIV-1 mRNAs. Although 1C8 was designed to be a structural mimic of the fused tetracyclic indole compound IDC16 that targets SRSF1, it did not affect the splice site shifting activity of SRSF1. Instead, 1C8 altered splicing regulation mediated by SRSF10. Depleting SRSF10 by RNA interference affected viral splicing and, like 1C8, decreased expression of Tat, Gag and Env. Incubating cells with 1C8 promoted the dephosphorylation of SRSF10 and increased its interaction with hTra2β, a protein previously implicated in the control of HIV-1 RNA splicing. While 1C8 affects the alternative splicing of cellular transcripts controlled by SRSF10 and hTra2β, concentrations greater than those needed to inhibit HIV-1 replication were required to elicit significant alterations. Thus, the ability of 1C8 to alter the SRSF10-dependent splicing of HIV-1 transcripts, with minor effects on cellular splicing, supports the view that SRSF10 may be used as a target for the development of new anti-viral agents.

Identifying subtype-specific associations between gene expression and DNA methylation profiles in breast cancer

BACKGROUND

Breast cancer is a complex disease in which different genomic patterns exists depending on different subtypes. Recent researches present that multiple subtypes of breast cancer occur at different rates, and play a crucial role in planning treatment. To better understand underlying biological mechanisms on breast cancer subtypes, investigating the specific gene regulatory system via different subtypes is desirable.

METHODS

Gene expression, as an intermediate phenotype, is estimated based on methylation profiles to identify the impact of epigenomic features on transcriptomic changes in breast cancer. We propose a kernel weighted l1-regularized regression model to incorporate tumor subtype information and further reveal gene regulations affected by different breast cancer subtypes. For the proper control of subtype-specific estimation, samples from different breast cancer subtype are learned at different rate based on target estimates. Kolmogorov Smirnov test is conducted to determine learning rate of each sample from different subtype.

RESULTS

It is observed that genes that might be sensitive to breast cancer subtype show prediction improvement when estimated using our proposed method. Comparing to a standard method, overall performance is also enhanced by incorporating tumor subtypes. In addition, we identified subtype-specific network structures based on the associations between gene expression and DNA methylation.

CONCLUSIONS

In this study, kernel weighted lasso model is proposed for identifying subtype-specific associations between gene expressions and DNA methylation profiles. Identification of subtype-specific gene expression associated with epigenomic changes might be helpful for better planning treatment and developing new therapies.

Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment

While intron retention (IR) is considered a widely conserved and distinct mechanism of gene expression control, its regulation is poorly understood. Here we show that DNA methylation directly regulates IR. We also find reduced occupancy of MeCP2 near the splice junctions of retained introns, mirroring the reduced DNA methylation at these sites. Accordingly, MeCP2 depletion in tissues and cells enhances IR. By analysing the MeCP2 interactome using mass spectrometry and RNA co-precipitation, we demonstrate that decreased MeCP2 binding near splice junctions facilitates IR via reduced recruitment of splicing factors, including Tra2b, and increased RNA polymerase II stalling. These results suggest an association between IR and a slower rate of transcription elongation, which reflects inefficient splicing factor recruitment. In summary, our results reinforce the interdependency between alternative splicing involving IR and epigenetic controls of gene expression.

Intron retention is a conserved mechanism that controls gene expression but its regulation is poorly understood. Here, the authors provide evidence that DNA methylation regulates intron retention and find reduced MeCP2 occupancy and splicing factor recruitment near affected splice junctions.

TRA2A Promoted Paclitaxel Resistance and Tumor Progression in Triple-Negative Breast Cancers via Regulating Alternative Splicing

Treatment of triple-negative breast cancer (TNBC) has been challenging, and paclitaxel resistance is one of the major obstacles to the better prognosis. Deregulation of alternative splicing (AS) may contribute to tumor progression and chemotherapy resistance. Human AS factor TRA2 has two separate gene paralogs encoding TRA2A and TRA2B proteins. TRA2B is associated with cancer cell survival and therapeutic sensitivity. However, the individual role of TRA2A in cancer progression has not been reported. Here we report that TRA2A facilitates proliferation and survival and migration and invasion of TNBC cells. In addition, TRA2A promotes paclitaxel resistance of TNBC by specifically controlling cancer-related splicing, which is independent of other splicing factors. TRA2A overexpression could promote AS of CALU, RSRC2, and PALM during paclitaxel treatment of TNBC cells. The isoform shift of RSRC2 from RSRC2s to RSRC2l leads to a decreased RSRC2 protein expression, which could contribute to TNBC paclitaxel resistance. TRA2A can regulate RSRC2 AS by specifically binding upstream intronic sequence of exon4. Strikingly, TRA2A expression is increased dramatically in patients with TNBC, and has a close relationship with decreased RSRC2 expression; both are associated with poor survival of TNBC. Collectively, our findings suggest that paclitaxel targets the TRA2A-RSRC2 splicing pathway, and deregulated TRA2A and RSRC2 expression may confer paclitaxel resistance. In addition to providing a novel molecular mechanism of cancer-related splicing dysregulation, our study demonstrates that expression of TRA2A in conjunction with RSRC2 may provide valuable molecular biomarker evidence for TNBC clinical treatment decisions and patient outcome. Mol Cancer Ther; 16(7); 1377-88. ©2017 AACR.

Myelodysplasia-associated mutations in serine/arginine-rich splicing factor SRSF2 lead to alternative splicing of CDC25C

BACKGROUND

Serine-arginine rich splicing factor 2 (SRSF2) is a protein known for its role in RNA splicing and genome stability. It has been recently discovered that SRSF2, along with other splicing regulators, is frequently mutated in patients with myelodysplastic syndrome (MDS). The most common MDS mutations in SRSF2 occur at proline 95; the mutant proteins are shown to have different RNA binding preferences, which may contribute to splicing changes detected in mutant cells. However, the influence of these SRSF2 MDS-associated mutations on specific splicing events remains poorly understood.

RESULTS

A tetracycline-inducible TF-1 erythroleukemia cell line was transduced with retroviruses to create cell lines expressing HA-tagged wildtype SRSF2, SRSF2 with proline 95 point mutations found in MDS, or SRSF2 with a deletion of one of the four major domains of the protein. Effects of these mutants on apoptosis and specific alternative splicing events were evaluated. Cells were also treated with DNA damaging drugs for comparison. MDS-related P95 point mutants of SRSF2 were expressed and phosphorylated at similar levels as wildtype SRSF2. However, cells expressing mutant SRSF2 exhibited higher levels of apoptosis than cells expressing wildtype SRSF2. Regarding alternative splicing events, in nearly all examined cases, SRSF2 P95 mutants acted in a similar fashion as the wildtype SRSF2. However, cells expressing SRSF2 P95 mutants had a percent increase in the C5 spliced isoform of cell division cycle 25C (CDC25C). The same alternative splicing of CDC25C was detected by treating cells with DNA damaging drugs, such as cisplatin, camptothecin, and trichostatin A at appropriate dosage. However, unlike DNA damaging drugs, SRSF2 P95 mutants did not activate the Ataxia telangiectasia mutated (ATM) pathway.

CONCLUSION

SRSF2 P95 mutants lead to alternative splicing of CDC25C in a manner that is not dependent on the DNA damage response.

Role of Bmznf-2, a Bombyx mori CCCH zinc finger gene, in masculinisation and differential splicing of Bmtra-2

Deciphering the regulatory factors involved in Bombyx mori sex determination has been a puzzle, challenging researchers for nearly a century now. The pre-mRNA of B. mori doublesex (Bmdsx), a master regulator gene of sexual differentiation, is differentially spliced, producing Bmdsxm and Bmdsxf transcripts in males and females respectively. The putative proteins encoded by these differential transcripts orchestrate antagonistic functions, which lead to sexual differentiation. A recent study in B. mori illustrated the role of a W-derived fem piRNA in conferring femaleness. In females, the fem piRNA was shown to suppress the activity of a Z-linked CCCH type zinc finger (znf) gene, Masculiniser (masc), which indirectly promotes the Bmdsxm type of splicing. In this study, we report a novel autosomal (Chr 25) CCCH type znf motif encoding gene Bmznf-2 as one of the potential factors in the Bmdsx sex specific differential splicing, and we also provide insights into its role in the alternative splicing of Bmtra2 by using ovary derived BmN cells. Over-expression of Bmznf-2 induced Bmdsxm type of splicing (masculinisation) with a correspondingly reduced expression of Bmdsxf type isoform in BmN cells. Further, the site-directed mutational studies targeting the tandem CCCH znf motifs revealed their indispensability in the observed phenotype of masculinisation. Additionally, the dual luciferase assays in BmN cells using 5' UTR region of the Bmznf-2 strongly implied the existence of a translational repression over this gene. From these findings, we propose Bmznf-2 to be one of the potential factors of masculinisation similar to Masc. From the growing number of Bmdsx splicing regulators, we assume that the sex determination cascade of B. mori is quite intricate in nature; hence, it has to be further investigated for its comprehensive understanding.

SURVIV for survival analysis of mRNA isoform variation

The rapid accumulation of clinical RNA-seq data sets has provided the opportunity to associate mRNA isoform variations to clinical outcomes. Here we report a statistical method SURVIV (Survival analysis of mRNA Isoform Variation), designed for identifying mRNA isoform variation associated with patient survival time. A unique feature and major strength of SURVIV is that it models the measurement uncertainty of mRNA isoform ratio in RNA-seq data. Simulation studies suggest that SURVIV outperforms the conventional Cox regression survival analysis, especially for data sets with modest sequencing depth. We applied SURVIV to TCGA RNA-seq data of invasive ductal carcinoma as well as five additional cancer types. Alternative splicing-based survival predictors consistently outperform gene expression-based survival predictors, and the integration of clinical, gene expression and alternative splicing profiles leads to the best survival prediction. We anticipate that SURVIV will have broad utilities for analysing diverse types of mRNA isoform variation in large-scale clinical RNA-seq projects.

Large-scale analysis of genome and transcriptome alterations in multiple tumors unveils novel cancer-relevant splicing networks

Alternative splicing is regulated by multiple RNA-binding proteins and influences the expression of most eukaryotic genes. However, the role of this process in human disease, and particularly in cancer, is only starting to be unveiled. We systematically analyzed mutation, copy number, and gene expression patterns of 1348 RNA-binding protein (RBP) genes in 11 solid tumor types, together with alternative splicing changes in these tumors and the enrichment of binding motifs in the alternatively spliced sequences. Our comprehensive study reveals widespread alterations in the expression of RBP genes, as well as novel mutations and copy number variations in association with multiple alternative splicing changes in cancer drivers and oncogenic pathways. Remarkably, the altered splicing patterns in several tumor types recapitulate those of undifferentiated cells. These patterns are predicted to be mainly controlled by MBNL1 and involve multiple cancer drivers, including the mitotic gene NUMA1 We show that NUMA1 alternative splicing induces enhanced cell proliferation and centrosome amplification in nontumorigenic mammary epithelial cells. Our study uncovers novel splicing networks that potentially contribute to cancer development and progression.

A tale of two paralogs: human Transformer2 proteins with differential RNA-binding affinities

The Transformer2 (Tra2) proteins in humans are homologues of the Drosophila Tra2 protein. One of the two RNA-binding paralogs, Tra2β, has been very well-studied over the past decade, but not much is known about Tra2α. It was very recently shown that the two proteins demonstrate the phenomenon of paralog compensation. Here, we provide a structural basis for this genetic backup circuit, using molecular modelling and dynamics studies. We show that the two proteins display similar binding specificities, but differential affinities to a short GAA-rich RNA stretch. Starting from the 6-nucleotide RNA in the solution structure, close to 4000 virtual mutations were modelled on RNA and the domain-RNA interactions were studied after energy minimisation to convergence. Separately, another known 13-nucleotide stretch was docked and the domain-RNA interactions were observed through a 100-ns dynamics trajectory. We have also demonstrated the 'compensatory' mechanism at the level of domains in one of the domain repeat-containing RNA-binding proteins.

The RNA Splicing Response to DNA Damage

The number of factors known to participate in the DNA damage response (DDR) has expanded considerably in recent years to include splicing and alternative splicing factors. While the binding of splicing proteins and ribonucleoprotein complexes to nascent transcripts prevents genomic instability by deterring the formation of RNA/DNA duplexes, splicing factors are also recruited to, or removed from, sites of DNA damage. The first steps of the DDR promote the post-translational modification of splicing factors to affect their localization and activity, while more downstream DDR events alter their expression. Although descriptions of molecular mechanisms remain limited, an emerging trend is that DNA damage disrupts the coupling of constitutive and alternative splicing with the transcription of genes involved in DNA repair, cell-cycle control and apoptosis. A better understanding of how changes in splice site selection are integrated into the DDR may provide new avenues to combat cancer and delay aging.

Identification of new biomarkers for human papillary thyroid carcinoma employing NanoString analysis

We previously reported an upregulation of the clock transcript BMAL1, correlating with TIMP1 expression in fresh-frozen samples from papillary thyroid carcinoma (PTC). Since frozen postoperative biopsy samples are difficult to obtain, we aimed to validate the application of high-precision NanoString analysis for formalin-fixed paraffin-embedded (FFPE) thyroid nodule samples and to screen for potential biomarkers associated with PTC. No significant differences were detected between fresh-frozen and FFPE samples. NanoString analysis of 51 transcripts in 17 PTC and 17 benign nodule samples obtained from different donors and in 24 pairs of benign and PTC nodules, obtained from the same donor (multinodular goiters), confirmed significant alterations in the levels of BMAL1, c-MET, c-KIT, TIMP1, and other transcripts. Moreover, we identified for the first time alterations in CHEK1 and BCL2 levels in PTC. A predictive score was established for each sample, based on the combined expression levels of BMAL1, CHEK1, c-MET, c-KIT and TIMP1. In combination with BRAF mutation analysis, this predictive score closely correlated with the clinicopathological characteristics of the analyzed thyroid nodules. Our study identified new thyroid transcripts with altered levels in PTC using the NanoString approach. A predictive score correlation coefficient might contribute to improve the preoperative diagnosis of thyroid nodules.

Deregulation of splicing factors and breast cancer development

It is well known that many genes implicated in the development and progression of breast cancer undergo aberrant alternative splicing events to produce proteins with pro-cancer properties. These changes in alternative splicing can arise from mutations or single-nucleotide polymorphisms (SNPs) within the DNA sequences of cancer-related genes, which can strongly affect the activity of splicing factors and influence the splice site choice. However, it is important to note that absence of mutations is not sufficient to prevent misleading choices in splice site selection. There is now increasing evidence to demonstrate that the expression profile of ten splicing factors (including SRs and hnRNPs) and eight RNA-binding proteins changes in breast cancer cells compared with normal cells. These modifications strongly influence the alternative splicing pattern of many cancer-related genes despite the absence of any detrimental mutations within their DNA sequences. Thus, a comprehensive assessment of the splicing factor status in breast cancer is important to provide insights into the mechanisms that lead to breast cancer development and metastasis. Whilst most studies focus on mutations that affect alternative splicing in cancer-related genes, this review focuses on splicing factors and RNA-binding proteins that are themselves deregulated in breast cancer and implicated in cancer-related alternative splicing events.