Long-read transcriptome sequencing of CLL and MDS patients uncovers molecular effects of SF3B1 mutations.

Mutations in splicing factor 3B subunit 1 (SF3B1) frequently occur in patients with chronic lymphocytic leukemia (CLL) and myelodysplastic syndromes (MDS). These mutations have different effects on the disease prognosis with beneficial effect in MDS and worse prognosis in CLL patients. A full-length transcriptome approach can expand our knowledge on SF3B1 mutation effects on RNA splicing and its contribution to patient survival and treatment options. We applied long-read transcriptome sequencing (LRTS) to 44 MDS and CLL patients, as well as two pairs of isogenic cell lines with and without SF3B1 mutations, and found >60% of novel isoforms. Splicing alterations were largely shared between cancer types and specifically affected the usage of introns and 3' splice sites. Our data highlighted a constrained window at canonical 3' splice sites in which dynamic splice site switches occurred in SF3B1-mutated patients. Using transcriptome-wide RNA binding maps and molecular dynamics simulations, we showed multimodal SF3B1 binding at 3' splice sites and predicted reduced RNA binding at the second binding pocket of SF3B1K700E Our work presents the hitherto most complete LRTS study of the SF3B1 mutation in CLL and MDS and provides a resource to study aberrant splicing in cancer. Moreover, we showed that different disease prognosis most likely results from the different cell types expanded during carcinogenesis rather than different mechanisms of action of the mutated SF3B1 These results have important implications for understanding the role of SF3B1 mutations in hematological malignancies and other related diseases.

Heterogenous nuclear ribonucleoprotein D-like controls endothelial cell functions.

HnRNPs are ubiquitously expressed RNA-binding proteins, tightly controlling posttranscriptional gene regulation. Consequently, hnRNP networks are essential for cellular homeostasis and their dysregulation is associated with cancer and other diseases. However, the physiological function of hnRNPs in non-cancerous cell systems are poorly understood. We analyzed the importance of HNRNPDL in endothelial cell functions. Knockdown of HNRNPDL led to impaired proliferation, migration and sprouting of spheroids. Transcriptome analysis identified cyclin D1 (CCND1) and tropomyosin 4 (TPM4) as targets of HNRNPDL, reflecting the phenotypic changes after knockdown. Our findings underline the importance of HNRNPDL for the homeostasis of physiological processes in endothelial cells.

FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns.

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

Poison cassette exon splicing of SRSF6 regulates nuclear speckle dispersal and the response to hypoxia.

Hypoxia induces massive changes in alternative splicing (AS) to adapt cells to the lack of oxygen. Here, we identify the splicing factor SRSF6 as a key factor in the AS response to hypoxia. The SRSF6 level is strongly reduced in acute hypoxia, which serves a dual purpose: it allows for exon skipping and triggers the dispersal of nuclear speckles. Our data suggest that cells use dispersal of nuclear speckles to reprogram their gene expression during hypoxic adaptation and that SRSF6 plays an important role in cohesion of nuclear speckles. Down-regulation of SRSF6 is achieved through inclusion of a poison cassette exon (PCE) promoted by SRSF4. Removing the PCE 3' splice site using CRISPR/Cas9 abolishes SRSF6 reduction in hypoxia. Aberrantly high SRSF6 levels in hypoxia attenuate hypoxia-mediated AS and impair dispersal of nuclear speckles. As a consequence, proliferation and genomic instability are increased, while the stress response is suppressed. The SRSF4-PCE-SRSF6 hypoxia axis is active in different cancer types, and high SRSF6 expression in hypoxic tumors correlates with a poor prognosis. We propose that the ultra-conserved PCE of SRSF6 acts as a tumor suppressor and that its inclusion in hypoxia is crucial to reduce SRSF6 levels. This may prevent tumor cells from entering the metastatic route of hypoxia adaptation.

High-throughput mutagenesis identifies mutations and RNA-binding proteins controlling CD19 splicing and CART-19 therapy resistance.

Following CART-19 immunotherapy for B-cell acute lymphoblastic leukaemia (B-ALL), many patients relapse due to loss of the cognate CD19 epitope. Since epitope loss can be caused by aberrant CD19 exon 2 processing, we herein investigate the regulatory code that controls CD19 splicing. We combine high-throughput mutagenesis with mathematical modelling to quantitatively disentangle the effects of all mutations in the region comprising CD19 exons 1-3. Thereupon, we identify ~200 single point mutations that alter CD19 splicing and thus could predispose B-ALL patients to developing CART-19 resistance. Furthermore, we report almost 100 previously unknown splice isoforms that emerge from cryptic splice sites and likely encode non-functional CD19 proteins. We further identify cis-regulatory elements and trans-acting RNA-binding proteins that control CD19 splicing (e.g., PTBP1 and SF3B4) and validate that loss of these factors leads to pervasive CD19 mis-splicing. Our dataset represents a comprehensive resource for identifying predictive biomarkers for CART-19 therapy.

Modulation of pre-mRNA structure by hnRNP proteins regulates alternative splicing of MALT1.

Alternative splicing plays key roles for cell type-specific regulation of protein function. It is controlled by cis-regulatory RNA elements that are recognized by RNA binding proteins (RBPs). The MALT1 paracaspase is a key factor of signaling pathways that mediate innate and adaptive immune responses. Alternative splicing of MALT1 is critical for controlling optimal T cell activation. We demonstrate that MALT1 splicing depends on RNA structural elements that sequester the splice sites of the alternatively spliced exon7. The RBPs hnRNP U and hnRNP L bind competitively to stem-loop RNA structures that involve the 5' and 3' splice sites flanking exon7. While hnRNP U stabilizes RNA stem-loop conformations that maintain exon7 skipping, hnRNP L disrupts these RNA elements to facilitate recruitment of the essential splicing factor U2AF2, thereby promoting exon7 inclusion. Our data represent a paradigm for the control of splice site selection by differential RBP binding and modulation of pre-mRNA structure.

Disruption of Prostaglandin E2 Signaling in Cancer-Associated Fibroblasts Limits Mammary Carcinoma Growth but Promotes Metastasis.

The activation and differentiation of cancer-associated fibroblasts (CAF) are involved in tumor progression. Here, we show that the tumor-promoting lipid mediator prostaglandin E2 (PGE2) plays a paradoxical role in CAF activation and tumor progression. Restricting PGE2 signaling via knockout of microsomal prostaglandin E synthase-1 (mPGES-1) in PyMT mice or of the prostanoid E receptor 3 (EP3) in CAFs stunted mammary carcinoma growth associated with strong CAF proliferation. CAF proliferation upon EP3 inhibition required p38 MAPK signaling. Mechanistically, TGFß-activated kinase-like protein (TAK1L), which was identified as a negative regulator of p38 MAPK activation, was decreased following ablation of mPGES-1 or EP3. In contrast with its effects on primary tumor growth, disruption of PGE2 signaling in CAFs induced epithelial-to-mesenchymal transition in cancer organoids and promoted metastasis in mice. Moreover, TAK1L expression in CAFs was associated with decreased CAF activation, reduced metastasis, and prolonged survival in human breast cancer. These data characterize a new pathway of regulating inflammatory CAF activation, which affects breast cancer progression. SIGNIFICANCE: The inflammatory lipid prostaglandin E2 suppresses cancer-associated fibroblast expansion and activation to limit primary mammary tumor growth while promoting metastasis.

Deep and accurate detection of m6A RNA modifications using miCLIP2 and m6Aboost machine learning.

N6-methyladenosine (m6A) is the most abundant internal RNA modification in eukaryotic mRNAs and influences many aspects of RNA processing. miCLIP (m6A individual-nucleotide resolution UV crosslinking and immunoprecipitation) is an antibody-based approach to map m6A sites with single-nucleotide resolution. However, due to broad antibody reactivity, reliable identification of m6A sites from miCLIP data remains challenging. Here, we present miCLIP2 in combination with machine learning to significantly improve m6A detection. The optimized miCLIP2 results in high-complexity libraries from less input material. Importantly, we established a robust computational pipeline to tackle the inherent issue of false positives in antibody-based m6A detection. The analyses were calibrated with Mettl3 knockout cells to learn the characteristics of m6A deposition, including m6A sites outside of DRACH motifs. To make our results universally applicable, we trained a machine learning model, m6Aboost, based on the experimental and RNA sequence features. Importantly, m6Aboost allows prediction of genuine m6A sites in miCLIP2 data without filtering for DRACH motifs or the need for Mettl3 depletion. Using m6Aboost, we identify thousands of high-confidence m6A sites in different murine and human cell lines, which provide a rich resource for future analysis. Collectively, our combined experimental and computational methodology greatly improves m6A identification.

Direct long-read RNA sequencing identifies a subset of questionable exitrons likely arising from reverse transcription artifacts.

Resistance to CD19-directed immunotherapies in lymphoblastic leukemia has been attributed, among other factors, to several aberrant CD19 pre-mRNA splicing events, including recently reported excision of a cryptic intron embedded within CD19 exon 2. While "exitrons" are known to exist in hundreds of human transcripts, we discovered, using reporter assays and direct long-read RNA sequencing (dRNA-seq), that the CD19 exitron is an artifact of reverse transcription. Extending our analysis to publicly available datasets, we identified dozens of questionable exitrons, dubbed "falsitrons," that appear only in cDNA-seq, but never in dRNA-seq. Our results highlight the importance of dRNA-seq for transcript isoform validation.

Exon Definition Facilitates Reliable Control of Alternative Splicing in the RON Proto-Oncogene.

Alternative splicing is a key step in eukaryotic gene expression that allows for the production of multiple transcript and protein isoforms from the same gene. Even though splicing is perturbed in many diseases, we currently lack insights into regulatory mechanisms promoting its precision and efficiency. We analyze high-throughput mutagenesis data obtained for an alternatively spliced exon in the proto-oncogene RON and determine the functional units that control this splicing event. Using mathematical modeling of distinct splicing mechanisms, we show that alternative splicing is based in RON on a so-called "exon definition" mechanism. Here, the recognition of the adjacent exons by the spliceosome is required for removal of an intron. We use our model to analyze the differences between the exon and intron definition scenarios and find that exon definition prevents the accumulation of deleterious, partially spliced retention products during alternative splicing regulation. Furthermore, it modularizes splicing control, as multiple regulatory inputs are integrated into a common net input, irrespective of the location and nature of the corresponding cis-regulatory elements in the pre-messenger RNA. Our analysis suggests that exon definition promotes robust and reliable splicing outcomes in RON splicing.

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.

Polypyrimidine tract-binding proteins are essential for B cell development.

Polypyrimidine tract-binding protein 1 (PTBP1) is a RNA-binding protein (RBP) expressed throughout B cell development. Deletion of Ptbp1 in mouse pro-B cells results in upregulation of PTBP2 and normal B cell development. We show that PTBP2 compensates for PTBP1 in B cell ontogeny as deletion of both Ptbp1 and Ptbp2 results in a complete block at the pro-B cell stage and a lack of mature B cells. In pro-B cells PTBP1 ensures precise synchronisation of the activity of cyclin dependent kinases at distinct stages of the cell cycle, suppresses S-phase entry and promotes progression into mitosis. PTBP1 controls mRNA abundance and alternative splicing of important cell cycle regulators including CYCLIN-D2, c-MYC, p107 and CDC25B. Our results reveal a previously unrecognised mechanism mediated by a RBP that is essential for B cell ontogeny and integrates transcriptional and post-translational determinants of progression through the cell cycle.

The RNA-binding ubiquitin ligase MKRN1 functions in ribosome-associated quality control of poly(A) translation.

BACKGROUND: Cells have evolved quality control mechanisms to ensure protein homeostasis by detecting and degrading aberrant mRNAs and proteins. A common source of aberrant mRNAs is premature polyadenylation, which can result in non-functional protein products. Translating ribosomes that encounter poly(A) sequences are terminally stalled, followed by ribosome recycling and decay of the truncated nascent polypeptide via ribosome-associated quality control. RESULTS: Here, we demonstrate that the conserved RNA-binding E3 ubiquitin ligase Makorin Ring Finger Protein 1 (MKRN1) promotes ribosome stalling at poly(A) sequences during ribosome-associated quality control. We show that MKRN1 directly binds to the cytoplasmic poly(A)-binding protein (PABPC1) and associates with polysomes. MKRN1 is positioned upstream of poly(A) tails in mRNAs in a PABPC1-dependent manner. Ubiquitin remnant profiling and in vitro ubiquitylation assays uncover PABPC1 and ribosomal protein RPS10 as direct ubiquitylation substrates of MKRN1. CONCLUSIONS: We propose that MKRN1 mediates the recognition of poly(A) tails to prevent the production of erroneous proteins from prematurely polyadenylated transcripts, thereby maintaining proteome integrity.

Phenotypic Plasticity of Fibroblasts during Mammary Carcinoma Development.

: Cancer-associated fibroblasts (CAFs) in the tumor microenvironment contribute to all stages of tumorigenesis and are usually considered to be tumor-promoting cells. CAFs show a remarkable degree of heterogeneity, which is attributed to developmental origin or to local environmental niches, resulting in distinct CAF subsets within individual tumors. While CAF heterogeneity is frequently investigated in late-stage tumors, data on longitudinal CAF development in tumors are lacking. To this end, we used the transgenic polyoma middle T oncogene-induced mouse mammary carcinoma model and performed whole transcriptome analysis in FACS-sorted fibroblasts from early- and late-stage tumors. We observed a shift in fibroblast populations over time towards a subset previously shown to negatively correlate with patient survival, which was confirmed by multispectral immunofluorescence analysis. Moreover, we identified a transcriptomic signature distinguishing CAFs from early- and late-stage tumors. Importantly, the signature of early-stage CAFs correlated well with tumor stage and survival in human mammary carcinoma patients. A random forest analysis suggested predictive value of the complete set of differentially expressed genes between early- and late-stage CAFs on bulk tumor patient samples, supporting the clinical relevance of our findings. In conclusion, our data show transcriptome alterations in CAFs during tumorigenesis in the mammary gland, which suggest that CAFs are educated by the tumor over time to promote tumor development. Moreover, we show that murine CAF gene signatures can harbor predictive value for human cancer.

A combined computational pipeline to detect circular RNAs in human cancer cells under hypoxic stress.

Hypoxia is associated with several diseases, including cancer. Cells that are deprived of adequate oxygen supply trigger transcriptional and post-transcriptional responses, which control cellular pathways such as angiogenesis, proliferation, and metabolic adaptation. Circular RNAs (circRNAs) are a novel class of mainly non-coding RNAs, which have been implicated in multiple cancers and attract increasing attention as potential biomarkers. Here, we characterize the circRNA signatures of three different cancer cell lines from cervical (HeLa), breast (MCF-7), and lung (A549) cancer under hypoxia. In order to reliably detect circRNAs, we integrate available tools with custom approaches for quantification and statistical analysis. Using this consolidated computational pipeline, we identify ~12000 circRNAs in the three cancer cell lines. Their molecular characteristics point to an involvement of complementary RNA sequences as well as trans-acting factors in circRNA biogenesis, such as the RNA-binding protein HNRNPC. Notably, we detect a number of circRNAs that are more abundant than their linear counterparts. In addition, 64 circRNAs significantly change in abundance upon hypoxia, in most cases in a cell type-specific manner. In summary, we present a comparative circRNA profiling in human cancer cell lines, which promises novel insights into the biogenesis and function of circRNAs under hypoxic stress.

Membrane-Associated RNA-Binding Proteins Orchestrate Organelle-Coupled Translation.

Proteins are positioned and act at defined subcellular locations. This is particularly important in eukaryotic cells that deliver proteins to membrane-bound organelles such as the endoplasmic reticulum (ER), mitochondria, or endosomes. It is axiomatic that organelle targeting depends mainly on polypeptide signals. However, recent results demonstrate that targeting elements within the encoding transcripts are essential for efficient protein localisation. Key readers of these elements are membrane-associated RNA-binding proteins (memRBPs) that orchestrate organelle-coupled translation. The translation products then either cross the membrane for organelle entry or hitchhike on organelle surfaces for complex assembly and co-transport. Understanding the interaction of protein- and RNA-based targeting signals is essential to decipher the molecular basis for mutant phenotypes in disease.

Decoding a cancer-relevant splicing decision in the RON proto-oncogene using high-throughput mutagenesis.

Mutations causing aberrant splicing are frequently implicated in human diseases including cancer. Here, we establish a high-throughput screen of randomly mutated minigenes to decode the cis-regulatory landscape that determines alternative splicing of exon 11 in the proto-oncogene MST1R (RON). Mathematical modelling of splicing kinetics enables us to identify more than 1000 mutations affecting RON exon 11 skipping, which corresponds to the pathological isoform RON∆165. Importantly, the effects correlate with RON alternative splicing in cancer patients bearing the same mutations. Moreover, we highlight heterogeneous nuclear ribonucleoprotein H (HNRNPH) as a key regulator of RON splicing in healthy tissues and cancer. Using iCLIP and synergy analysis, we pinpoint the functionally most relevant HNRNPH binding sites and demonstrate how cooperative HNRNPH binding facilitates a splicing switch of RON exon 11. Our results thereby offer insights into splicing regulation and the impact of mutations on alternative splicing in cancer.

Insights into the design and interpretation of iCLIP experiments.

BACKGROUND: Ultraviolet (UV) crosslinking and immunoprecipitation (CLIP) identifies the sites on RNAs that are in direct contact with RNA-binding proteins (RBPs). Several variants of CLIP exist, which require different computational approaches for analysis. This variety of approaches can create challenges for a novice user and can hamper insights from multi-study comparisons. Here, we produce data with multiple variants of CLIP and evaluate the data with various computational methods to better understand their suitability. RESULTS: We perform experiments for PTBP1 and eIF4A3 using individual-nucleotide resolution CLIP (iCLIP), employing either UV-C or photoactivatable 4-thiouridine (4SU) combined with UV-A crosslinking and compare the results with published data. As previously noted, the positions of complementary DNA (cDNA)-starts depend on cDNA length in several iCLIP experiments and we now find that this is caused by constrained cDNA-ends, which can result from the sequence and structure constraints of RNA fragmentation. These constraints are overcome when fragmentation by RNase I is efficient and when a broad cDNA size range is obtained. Our study also shows that if RNase does not efficiently cut within the binding sites, the original CLIP method is less capable of identifying the longer binding sites of RBPs. In contrast, we show that a broad size range of cDNAs in iCLIP allows the cDNA-starts to efficiently delineate the complete RNA-binding sites. CONCLUSIONS: We demonstrate the advantage of iCLIP and related methods that can amplify cDNAs that truncate at crosslink sites and we show that computational analyses based on cDNAs-starts are appropriate for such methods.

The RNA-binding protein HuR is essential for the B cell antibody response.

Post-transcriptional regulation of mRNA by the RNA-binding protein HuR (encoded by Elavl1) is required in B cells for the germinal center reaction and for the production of class-switched antibodies in response to thymus-independent antigens. Transcriptome-wide examination of RNA isoforms and their abundance and translation in HuR-deficient B cells, together with direct measurements of HuR-RNA interactions, revealed that HuR-dependent splicing of mRNA affected hundreds of transcripts, including that encoding dihydrolipoamide S-succinyltransferase (Dlst), a subunit of the 2-oxoglutarate dehydrogenase (α-KGDH) complex. In the absence of HuR, defective mitochondrial metabolism resulted in large amounts of reactive oxygen species and B cell death. Our study shows how post-transcriptional processes control the balance of energy metabolism required for the proliferation and differentiation of B cells.

Analysis of alternative splicing associated with aging and neurodegeneration in the human brain.

Age is the most important risk factor for neurodegeneration; however, the effects of aging and neurodegeneration on gene expression in the human brain have most often been studied separately. Here, we analyzed changes in transcript levels and alternative splicing in the temporal cortex of individuals of different ages who were cognitively normal, affected by frontotemporal lobar degeneration (FTLD), or affected by Alzheimer's disease (AD). We identified age-related splicing changes in cognitively normal individuals and found that these were present also in 95% of individuals with FTLD or AD, independent of their age. These changes were consistent with increased polypyrimidine tract binding protein (PTB)-dependent splicing activity. We also identified disease-specific splicing changes that were present in individuals with FTLD or AD, but not in cognitively normal individuals. These changes were consistent with the decreased neuro-oncological ventral antigen (NOVA)-dependent splicing regulation, and the decreased nuclear abundance of NOVA proteins. As expected, a dramatic down-regulation of neuronal genes was associated with disease, whereas a modest down-regulation of glial and neuronal genes was associated with aging. Whereas our data indicated that the age-related splicing changes are regulated independently of transcript-level changes, these two regulatory mechanisms affected expression of genes with similar functions, including metabolism and DNA repair. In conclusion, the alternative splicing changes identified in this study provide a new link between aging and neurodegeneration.