Machine learning-optimized targeted detection of alternative splicing.

RNA-sequencing (RNA-seq) is widely adopted for transcriptome analysis but has inherent biases which hinder the comprehensive detection and quantification of alternative splicing. To address this, we present an efficient targeted RNA-seq method that greatly enriches for splicing-informative junction-spanning reads. Local Splicing Variation sequencing (LSV-seq) utilizes multiplexed reverse transcription from highly scalable pools of primers anchored near splicing events of interest. Primers are designed using Optimal Prime, a novel machine learning algorithm trained on the performance of thousands of primer sequences. In experimental benchmarks, LSV-seq achieves high on-target capture rates and concordance with RNA-seq, while requiring significantly lower sequencing depth. Leveraging deep learning splicing code predictions, we used LSV-seq to target events with low coverage in GTEx RNA-seq data and newly discover hundreds of tissue-specific splicing events. Our results demonstrate the ability of LSV-seq to quantify splicing of events of interest at high-throughput and with exceptional sensitivity.

An anciently diverged family of RNA binding proteins maintain correct splicing of a class of ultra-long exons through cryptic splice site repression.

Previously, we showed that the germ cell-specific nuclear protein RBMXL2 represses cryptic splicing patterns during meiosis and is required for male fertility (Ehrmann et al., 2019). Here, we show that in somatic cells the similar yet ubiquitously expressed RBMX protein has similar functions. RBMX regulates a distinct class of exons that exceed the median human exon size. RBMX protein-RNA interactions are enriched within ultra-long exons, particularly within genes involved in genome stability, and repress the selection of cryptic splice sites that would compromise gene function. The RBMX gene is silenced during male meiosis due to sex chromosome inactivation. To test whether RBMXL2 might replace the function of RBMX during meiosis we induced expression of RBMXL2 and the more distantly related RBMY protein in somatic cells, finding each could rescue aberrant patterns of RNA processing caused by RBMX depletion. The C-terminal disordered domain of RBMXL2 is sufficient to rescue proper splicing control after RBMX depletion. Our data indicate that RBMX and RBMXL2 have parallel roles in somatic tissues and the germline that must have been conserved for at least 200 million years of mammalian evolution. We propose RBMX family proteins are particularly important for the splicing inclusion of some ultra-long exons with increased intrinsic susceptibility to cryptic splice site selection.

Trametinib Sensitivity is Defined by a Myeloid Differentiation Profile in Acute Myeloid Leukemia.

BACKGROUND AND OBJECTIVE: Acute myelogenous leukemia (AML) is a common blood cancer marked by heterogeneity in disease and diverse genetic abnormalities. Additional therapies are needed as the 5-year survival remains below 30%. Trametinib is a mitogen-activated extracellular signal-regulated kinase (MEK) inhibitor that is widely used in solid tumors and also in tumors with activating RAS mutations. A subset of patients with AML carry activating RAS mutations; however, a small-scale clinical trial with trametinib showed little efficacy. Here, we sought to identify transcriptomic determinants of trametinib sensitivity in AML. METHODS: We tested the activity of trametinib against a panel of tumor cells from patients with AML ex vivo and compared this with RNA sequencing (RNA-Seq) data from untreated blasts from the same patient samples. We then used a correlation analysis between gene expression and trametinib sensitivity to identify potential biomarkers predictive of drug response. RESULTS: We found that a subset of AML tumor cells were sensitive to trametinib ex vivo, only a fraction of which (3/10) carried RAS mutations. On the basis of our RNA-Seq analysis we found that markers of trametinib sensitivity are associated with a myeloid differentiation profile that includes high expression of CD14 and CLEC7A (Dectin-1), similar to the gene expression profile of monocytes. Further characterization confirmed that trametinib-sensitive samples display features of monocytic differentiation with high CD14 surface expression and were enriched for the M4 subtypes of the FAB classification. CONCLUSIONS: Our study identifies additional molecular markers that can be used with molecular features including RAS status to identify patients with AML that may benefit from trametinib treatment.

A Deep Dive into Statistical Modeling of RNA Splicing QTLs Reveals New Variants that Explain Neurodegenerative Disease.

Genome-wide association studies (GWAS) have identified thousands of putative disease causing variants with unknown regulatory effects. Efforts to connect these variants with splicing quantitative trait loci (sQTLs) have provided functional insights, yet sQTLs reported by existing methods cannot explain many GWAS signals. We show current sQTL modeling approaches can be improved by considering alternative splicing representation, model calibration, and covariate integration. We then introduce MAJIQTL, a new pipeline for sQTL discovery. MAJIQTL includes two new statistical methods: a weighted multiple testing approach for sGene discovery and a model for sQTL effect size inference to improve variant prioritization. By applying MAJIQTL to GTEx, we find significantly more sGenes harboring sQTLs with functional significance. Notably, our analysis implicates the novel variant rs582283 in Alzheimer's disease. Using antisense oligonucleotides, we validate this variant's effect by blocking the implicated YBX3 binding site, leading to exon skipping in the gene MS4A3.

Contrasting and combining transcriptome complexity captured by short and long RNA sequencing reads.

Mapping transcriptomic variations using either short- or long-read RNA sequencing is a staple of genomic research. Long reads are able to capture entire isoforms and overcome repetitive regions, whereas short reads still provide improved coverage and error rates. Yet, open questions remain, such as how to quantitatively compare the technologies, can we combine them, and what is the benefit of such a combined view? We tackle these questions by first creating a pipeline to assess matched long- and short-read data using a variety of transcriptome statistics. We find that across data sets, algorithms, and technologies, matched short-read data detects ∼30% more splice junctions, such that ∼10%-30% of the splice junctions included at ≥20% by short reads are missed by long reads. In contrast, long reads detect many more intron-retention events and can detect full isoforms, pointing to the benefit of combining the technologies. We introduce MAJIQ-L, an extension of the MAJIQ software, to enable a unified view of transcriptome variations from both technologies and demonstrate its benefits. Our software can be used to assess any future long-read technology or algorithm and can be combined with short-read data for improved transcriptome analysis.

An Alternatively Spliced Gain-of-Function NT5C2 Isoform Contributes to Chemoresistance in Acute Lymphoblastic Leukemia.

Relapsed or refractory B-cell acute lymphoblastic leukemia (B-ALL) is a major cause of pediatric cancer-related deaths. Relapse-specific mutations do not account for all chemotherapy failures in B- ALL patients, suggesting additional mechanisms of resistance. By mining RNA-seq datasets of paired diagnostic/relapse pediatric B-ALL samples, we discovered pervasive alternative splicing (AS) patterns linked to relapse and affecting drivers of resistance to glucocorticoids, anti-folates, and thiopurines. Most splicing variations represented cassette exon skipping, "poison" exon inclusion, and intron retention, phenocopying well-documented loss-of-function mutations. In contrast, relapse-associated AS of NT5C2 mRNA yielded an isoform with the functionally uncharacterized in-frame exon 6a. Incorporation of the 8-amino acid sequence SQVAVQKR into this enzyme created a putative phosphorylation site and resulted in elevated nucleosidase activity, which is a known consequence of gain-of-function mutations in NT5C2 and a common determinant of 6-mercaptopurine (6-MP) resistance. Consistent with this finding, NT5C2ex6a and the R238W hotspot variant conferred comparable levels of resistance to 6-MP in B-ALL cells both in vitro and in vivo. Furthermore, both the NT5C2ex6a and R238W variants induced collateral sensitivity to the inosine monophosphate dehydrogenase (IMPDH) inhibitor mizoribine. These results ascribe an important role for splicing perturbations in chemotherapy resistance in relapsed B-ALL and suggest that IMPDH inhibitors, including the commonly used immunosuppressive agent mycophenolate mofetil, could be a valuable therapeutic option for treating thiopurine-resistant leukemias.

An Alternatively Spliced Gain-of-Function NT5C2 Isoform Contributes to Chemoresistance in Acute Lymphoblastic Leukemia.

Relapsed or refractory B-cell acute lymphoblastic leukemia (B-ALL) is a major cause of pediatric cancer-related deaths. Relapse-specific mutations do not account for all chemotherapy failures in B-ALL patients, suggesting additional mechanisms of resistance. By mining RNA sequencing datasets of paired diagnostic/relapse pediatric B-ALL samples, we discovered pervasive alternative splicing (AS) patterns linked to relapse and affecting drivers of resistance to glucocorticoids, antifolates, and thiopurines. Most splicing variations represented cassette exon skipping, "poison" exon inclusion, and intron retention, phenocopying well-documented loss-of-function mutations. In contrast, relapse-associated AS of NT5C2 mRNA yielded an isoform with the functionally uncharacterized in-frame exon 6a. Incorporation of the 8-amino acid sequence SQVAVQKR into this enzyme created a putative phosphorylation site and resulted in elevated nucleosidase activity, which is a known consequence of gain-of-function mutations in NT5C2 and a common determinant of 6-mercaptopurine resistance. Consistent with this finding, NT5C2ex6a and the R238W hotspot variant conferred comparable levels of resistance to 6-mercaptopurine in B-ALL cells both in vitro and in vivo. Furthermore, both NT5C2ex6a and the R238W variant induced collateral sensitivity to the inosine monophosphate dehydrogenase inhibitor mizoribine. These results ascribe to splicing perturbations an important role in chemotherapy resistance in relapsed B-ALL and suggest that inosine monophosphate dehydrogenase inhibitors, including the commonly used immunosuppressive agent mycophenolate mofetil, could be a valuable therapeutic option for treating thiopurine-resistant leukemias. Significance: Alternative splicing is a potent mechanism of acquired drug resistance in relapsed/refractory acute lymphoblastic leukemias that has diagnostic and therapeutic implications for patients who lack mutations in known chemoresistance genes.

MAJIQlopedia: an encyclopedia of RNA splicing variations in human tissues and cancer.

Quantification of RNA splicing variations based on RNA-Sequencing can reveal tissue- and disease-specific splicing patterns. To study such splicing variations, we introduce MAJIQlopedia, an encyclopedia of splicing variations that encompasses 86 human tissues and 41 cancer datasets. MAJIQlopedia reports annotated and unannotated splicing events for a total of 486 175 alternative splice junctions in normal tissues and 338 317 alternative splice junctions in cancer. This database, available at https://majiq.biociphers.org/majiqlopedia/, includes a user-friendly interface that provides graphical representations of junction usage quantification for each junction across all tissue or cancer types. To demonstrate case usage of MAJIQlopedia, we review splicing variations in genes WT1, MAPT and BIN1, which all have known tissue or cancer-specific splicing variations. We also use MAJIQlopedia to highlight novel splicing variations in FDX1 and MEGF9 in normal tissues, and we uncover a novel exon inclusion event in RPS6KA6 that only occurs in two cancer types. Users can download the database, request the addition of data to the webtool, or install a MAJIQlopedia server to integrate proprietary data. MAJIQlopedia can serve as a reference database for researchers seeking to understand what splicing variations exist in genes of interest, and those looking to understand tissue- or cancer-specific splice isoform usage.

Contrasting and Combining Transcriptome Complexity Captured by Short and Long RNA Sequencing Reads.

Mapping transcriptomic variations using either short or long reads RNA sequencing is a staple of genomic research. Long reads are able to capture entire isoforms and overcome repetitive regions, while short reads still provides improved coverage and error rates. Yet how to quantitatively compare the technologies, can we combine those, and what may be the benefit of such a combined view remain open questions. We tackle these questions by first creating a pipeline to assess matched long and short reads data using a variety of transcriptome statistics. We find that across datasets, algorithms and technologies, matched short reads data detects roughly 50% more splice junctions, with 10-30% of the splice junctions included at 20% or more are missed by long reads. In contrast, long reads detect many more intron retention events, pointing to the benefit of combining the technologies. We introduce MAJIQ-L, an extension of the MAJIQ software to enable a unified view of transcriptome variations from both technologies and demonstrate its benefits. Our software can be used to assess any future long reads technology or algorithm, and combine it with short reads data for improved transcriptome analysis.

Extensible benchmarking of methods that identify and quantify polyadenylation sites from RNA-seq data.

The tremendous rate with which data is generated and analysis methods emerge makes it increasingly difficult to keep track of their domain of applicability, assumptions, limitations, and consequently, of the efficacy and precision with which they solve specific tasks. Therefore, there is an increasing need for benchmarks, and for the provision of infrastructure for continuous method evaluation. APAeval is an international community effort, organized by the RNA Society in 2021, to benchmark tools for the identification and quantification of the usage of alternative polyadenylation (APA) sites from short-read, bulk RNA-sequencing (RNA-seq) data. Here, we reviewed 17 tools and benchmarked eight on their ability to perform APA identification and quantification, using a comprehensive set of RNA-seq experiments comprising real, synthetic, and matched 3'-end sequencing data. To support continuous benchmarking, we have incorporated the results into the OpenEBench online platform, which allows for continuous extension of the set of methods, metrics, and challenges. We envisage that our analyses will assist researchers in selecting the appropriate tools for their studies, while the containers and reproducible workflows could easily be deployed and extended to evaluate new methods or data sets.

Alternative splicing of its 5'-UTR limits CD20 mRNA translation and enables resistance to CD20-directed immunotherapies.

Aberrant skipping of coding exons in CD19 and CD22 compromises the response to immunotherapy in B-cell malignancies. Here, we showed that the MS4A1 gene encoding human CD20 also produces several messenger RNA (mRNA) isoforms with distinct 5' untranslated regions. Four variants (V1-4) were detected using RNA sequencing (RNA-seq) at distinct stages of normal B-cell differentiation and B-lymphoid malignancies, with V1 and V3 being the most abundant. During B-cell activation and Epstein-Barr virus infection, redirection of splicing from V1 to V3 coincided with increased CD20 positivity. Similarly, in diffuse large B-cell lymphoma, only V3, but not V1, correlated with CD20 protein levels, suggesting that V1 might be translation-deficient. Indeed, the longer V1 isoform contained upstream open reading frames and a stem-loop structure, which cooperatively inhibited polysome recruitment. By modulating CD20 isoforms with splice-switching morpholino oligomers, we enhanced CD20 expression and anti-CD20 antibody rituximab-mediated cytotoxicity in a panel of B-cell lines. Furthermore, reconstitution of CD20-knockout cells with V3 mRNA led to the recovery of CD20 positivity, whereas V1-reconstituted cells had undetectable levels of CD20 protein. Surprisingly, in vitro CD20-directed chimeric antigen receptor T cells were able to kill both V3- and V1-expressing cells, but the bispecific T-cell engager mosunetuzumab was only effective against V3-expressing cells. To determine whether CD20 splicing is involved in immunotherapy resistance, we performed RNA-seq on 4 postmosunetuzumab follicular lymphoma relapses and discovered that in 2 of them, the downregulation of CD20 was accompanied by a V3-to-V1 shift. Thus, splicing-mediated mechanisms of epitope loss extend to CD20-directed immunotherapies.

Parkin Deficiency Suppresses Antigen Presentation to Promote Tumor Immune Evasion and Immunotherapy Resistance.

Parkin is an E3 ubiquitin ligase, which plays a key role in the development of Parkinson disease. Parkin defects also occur in numerous cancers, and a growing body of evidence indicates that Parkin functions as a tumor suppressor that impedes a number of cellular processes involved in tumorigenesis. Here, we generated murine and human models that closely mimic the advanced-stage tumors where Parkin deficiencies are found to provide deeper insights into the tumor suppressive functions of Parkin. Loss of Parkin expression led to aggressive tumor growth, which was associated with poor tumor antigen presentation and limited antitumor CD8+ T-cell infiltration and activation. The effect of Parkin deficiency on tumor growth was lost following depletion of CD8+ T cells. In line with previous findings, Parkin deficiency was linked with mitochondria-associated metabolic stress, PTEN degradation, and enhanced Akt activation. Increased Akt signaling led to dysregulation of antigen presentation, and treatment with the Akt inhibitor MK2206-2HCl restored antigen presentation in Parkin-deficient tumors. Analysis of data from patients with clear cell renal cell carcinoma indicated that Parkin expression was downregulated in tumors and that low expression correlated with reduced overall survival. Furthermore, low Parkin expression correlated with reduced patient response to immunotherapy. Overall, these results identify a role for Parkin deficiency in promoting tumor immune evasion that may explain the poor prognosis associated with loss of Parkin across multiple types of cancer. SIGNIFICANCE: Parkin prevents immune evasion by regulating tumor antigen processing and presentation through the PTEN/Akt network, which has important implications for immunotherapy treatments in patients with Parkin-deficient tumors.

Alternative splicing of its 5'-UTR limits CD20 mRNA translation and enables resistance to CD20-directed immunotherapies.

Aberrant skipping of coding exons in CD19 and CD22 compromises responses to immunotherapy for B-cell malignancies. Here, we show that the MS4A1 gene encoding human CD20 also produces several mRNA isoforms with distinct 5' untranslated regions (5'-UTR). Four variants (V1-4) were detectable by RNA-seq in distinct stages of normal B-cell differentiation and B-lymphoid malignancies, with V1 and V3 being the most abundant by far. During B-cell activation and Epstein-Barr virus infection, redirection of splicing from V1 to V3 coincided with increased CD20 positivity. Similarly, in diffuse large B-cell lymphoma only V3, but not V1, correlated with CD20 protein levels, suggesting that V1 might be translation-deficient. Indeed, the longer V1 isoform was found to contain upstream open reading frames (uORFs) and a stem-loop structure, which cooperatively inhibited polysome recruitment. By modulating CD20 isoforms with splice-switching Morpholino oligomers, we enhanced CD20 expression and anti-CD20 antibody rituximab-mediated cytotoxicity in a panel of B-cell lines. Furthermore, reconstitution of CD20-knockout cells with V3 mRNA led to the recovery of CD20 positivity, while V1-reconstituted cells had undetectable levels of CD20 protein. Surprisingly, in vitro CD20-directed CAR T cells were able to kill both V3- and V1-expressing cells, but the bispecific T cell engager mosunetuzumab was only effective against V3-expressing cells. To determine whether CD20 splicing is involved in immunotherapy resistance, we performed RNA-seq on four post-mosunetuzumab follicular lymphoma relapses and discovered that in two of them downregulation of CD20 was accompanied by the V3-to-V1 shift. Thus, splicing-mediated mechanisms of epitope loss extend to CD20-directed immunotherapies. Key Points: In normal & malignant human B cells, CD20 mRNA is alternatively spliced into four 5'-UTR isoforms, some of which are translation-deficient.The balance between translation-deficient and -competent isoforms modulates CD20 protein levels & responses to CD20-directed immunotherapies. Explanation of Novelty: We discovered that in normal and malignant B-cells, CD20 mRNA is alternatively spliced to generate four distinct 5'-UTRs, including the longer translation-deficient V1 variant. Cells predominantly expressing V1 were still sensitive to CD20-targeting chimeric antigen receptor T-cells. However, they were resistant to the bispecific anti-CD3/CD20 antibody mosunetuzumab, and the shift to V1 were observed in CD20-negative post-mosunetuzumab relapses of follicular lymphoma.

Extensible benchmarking of methods that identify and quantify polyadenylation sites from RNA-seq data.

The tremendous rate with which data is generated and analysis methods emerge makes it increasingly difficult to keep track of their domain of applicability, assumptions, and limitations and consequently, of the efficacy and precision with which they solve specific tasks. Therefore, there is an increasing need for benchmarks, and for the provision of infrastructure for continuous method evaluation. APAeval is an international community effort, organized by the RNA Society in 2021, to benchmark tools for the identification and quantification of the usage of alternative polyadenylation (APA) sites from short-read, bulk RNA-sequencing (RNA-seq) data. Here, we reviewed 17 tools and benchmarked eight on their ability to perform APA identification and quantification, using a comprehensive set of RNA-seq experiments comprising real, synthetic, and matched 3'-end sequencing data. To support continuous benchmarking, we have incorporated the results into the OpenEBench online platform, which allows for seamless extension of the set of methods, metrics, and challenges. We envisage that our analyses will assist researchers in selecting the appropriate tools for their studies. Furthermore, the containers and reproducible workflows generated in the course of this project can be seamlessly deployed and extended in the future to evaluate new methods or datasets.

The Tumor Suppressor Adenomatous Polyposis Coli (apc) Is Required for Neural Crest-Dependent Craniofacial Development in Zebrafish.

Neural crest (NC) is a unique vertebrate cell type arising from the border of the neural plate and epidermis that gives rise to diverse tissues along the entire body axis. Roberto Mayor and colleagues have made major contributions to our understanding of NC induction, delamination, and migration. We report that a truncating mutation of the classical tumor suppressor Adenomatous Polyposis Coli (apc) disrupts craniofacial development in zebrafish larvae, with a marked reduction in the cranial neural crest (CNC) cells that contribute to mandibular and hyoid pharyngeal arches. While the mechanism is not yet clear, the altered expression of signaling molecules that guide CNC migration could underlie this phenotype. For example, apcmcr/mcr larvae express substantially higher levels of complement c3, which Mayor and colleagues showed impairs CNC cell migration when overexpressed. However, we also observe reduction in stroma-derived factor 1 (sdf1/cxcl12), which is required for CNC migration into the head. Consistent with our previous work showing that APC directly enhances the activity of glycogen synthase kinase 3 (GSK-3) and, independently, that GSK-3 phosphorylates multiple core mRNA splicing factors, we identify 340 mRNA splicing variations in apc mutant zebrafish, including a splice variant that deletes a conserved domain in semaphorin 3f (sema3f), an axonal guidance molecule and a known regulator of CNC migration. Here, we discuss potential roles for apc in CNC development in the context of some of the seminal findings of Mayor and colleagues.

RNA splicing analysis using heterogeneous and large RNA-seq datasets.

The ubiquity of RNA-seq has led to many methods that use RNA-seq data to analyze variations in RNA splicing. However, available methods are not well suited for handling heterogeneous and large datasets. Such datasets scale to thousands of samples across dozens of experimental conditions, exhibit increased variability compared to biological replicates, and involve thousands of unannotated splice variants resulting in increased transcriptome complexity. We describe here a suite of algorithms and tools implemented in the MAJIQ v2 package to address challenges in detection, quantification, and visualization of splicing variations from such datasets. Using both large scale synthetic data and GTEx v8 as benchmark datasets, we assess the advantages of MAJIQ v2 compared to existing methods. We then apply MAJIQ v2 package to analyze differential splicing across 2,335 samples from 13 brain subregions, demonstrating its ability to offer insights into brain subregion-specific splicing regulation.

A Bayesian model for unsupervised detection of RNA splicing based subtypes in cancers.

Identification of cancer sub-types is a pivotal step for developing personalized treatment. Specifically, sub-typing based on changes in RNA splicing has been motivated by several recent studies. We thus develop CHESSBOARD, an unsupervised algorithm tailored for RNA splicing data that captures "tiles" in the data, defined by a subset of unique splicing changes in a subset of patients. CHESSBOARD allows for a flexible number of tiles, accounts for uncertainty of splicing quantification, and is able to model missing values as additional signals. We first apply CHESSBOARD to synthetic data to assess its domain specific modeling advantages, followed by analysis of several leukemia datasets. We show detected tiles are reproducible in independent studies, investigate their possible regulatory drivers and probe their relation to known AML mutations. Finally, we demonstrate the potential clinical utility of CHESSBOARD by supplementing mutation based diagnostic assays with discovered splicing profiles to improve drug response correlation.

FBXW7ß isoform drives transcriptional activation of the proinflammatory TNF cluster in human pro-B cells.

Noncanonical exon usage plays many important roles in cellular phenotypes, but its contribution to human B-cell development remains sketchily understood. To fill this gap, we collected various B-cell fractions from bone marrow (BM) and tonsil donors, performed RNA sequencing, and examined transcript variants. We identified 150 genes that harbor local splicing variations in all pairwise comparisons. One of them encodes FBXW7, an E3 ubiquitin ligase implicated as a driver in several blood cancers. Surprisingly, we discovered that in normal human pro-B cells, the predominant transcript used an alternative first exon to produce the poorly characterized FBXW7ß isoform, previously thought to be restricted to neural tissues. The FBXW7ß transcript was also abundant in cell lines and primary samples of pediatric B-cell acute lymphoblastic leukemia (B-ALL), which originates in the BM. When overexpressed in a heterologous cell system, this transcript yielded the expected protein product, as judged by anti-FLAG immunoblotting and mass spectrometry. Furthermore, in REH B-ALL cells, FBXW7ß mRNA was the only FBXW7 isoform enriched in the polyribosome fraction. To shed light on possible functions of FBXW7ß, we used gain- and loss-of-function approaches and identified an FBXW7-dependent inflammatory gene signature, apparent in a subset of B-ALL with high FBXW7ß expression. This signature contained several members of the tumor necrosis factor superfamily, including those comprising the HLA Class III cluster (LTB, LST1, NCR3, LTA, and NFKBIL1). Our findings suggest that FBXW7ß expression drives proinflammatory responses, which could contribute to normal B-cell development, leukemogenesis, and responses to anticancer therapies.

The global Protein-RNA interaction map of ESRP1 defines a post-transcriptional program that is essential for epithelial cell function.

The epithelial splicing regulatory proteins, ESRP1 and ESRP2, are essential for mammalian development through the regulation of a global program of alternative splicing of genes involved in the maintenance of epithelial cell function. To further inform our understanding of the molecular functions of ESRP1, we performed enhanced crosslinking immunoprecipitation coupled with high-throughput sequencing (eCLIP) in epithelial cells of mouse epidermis. The genome-wide binding sites of ESRP1 were integrated with RNA-Seq analysis of alterations in splicing and total gene expression that result from epidermal ablation of Esrp1 and Esrp2. These studies demonstrated that ESRP1 functions in splicing regulation occur primarily through direct binding in a position-dependent manner to promote either exon inclusion or skipping. In addition, we also identified widespread binding of ESRP1 in 3' and 5' untranslated regions (UTRs) of genes involved in epithelial cell function, suggesting that its post-transcriptional functions extend beyond splicing regulation.

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.