Alternative pre-mRNA splicing has long been proposed to contribute greatly to proteome complexity. However, the extent to which mature mRNA isoforms are successfully translated into protein remains controversial. Here, we used high-throughput RNA sequencing and mass spectrometry (MS)-based proteomics to better evaluate the translation of alternatively spliced mRNAs. To increase proteome coverage and improve protein quantitation, we optimized cell fractionation and sample processing steps at both the protein and peptide level. Furthermore, we generated a custom peptide database trained on analysis of RNA-seq data with MAJIQ, an algorithm optimized to detect and quantify differential and unannotated splice junction usage. We matched tandem mass spectra acquired by data-dependent acquisition (DDA) against our custom RNA-seq based database, as well as SWISS-PROT and RefSeq databases to improve identification of splicing-derived proteoforms by 28% compared with use of the SWISS-PROT database alone. Altogether, we identified peptide evidence for 554 alternate proteoforms corresponding to 274 genes. Our increased depth and detection of proteins also allowed us to track changes in the transcriptome and proteome induced by T-cell stimulation, as well as fluctuations in protein subcellular localization. In sum, our data here confirm that use of generic databases in proteomic studies underestimates the number of spliced mRNA isoforms that are translated into protein and provides a workflow that improves isoform detection in large-scale proteomic experiments.

The effects of confounding factors on gene expression analysis have been extensively studied following the introduction of high-throughput microarrays and subsequently RNA sequencing. In contrast, there is a lack of equivalent analysis and tools for RNA splicing. Here we first assess the effect of confounders on both expression and splicing quantifications in two large public RNA-Seq datasets (TARGET, ENCODE). We show quantification of splicing variations are affected at least as much as those of gene expression, revealing unwanted sources of variations in both datasets. Next, we develop MOCCASIN, a method to correct the effect of both known and unknown confounders on RNA splicing quantification and demonstrate MOCCASIN's effectiveness on both synthetic and real data. Code, synthetic and corrected datasets are all made available as resources.

Alternative splicing contributes to gene expression dynamics in many tissues, yet its role in auditory development remains unclear. We performed whole-exome sequencing in individuals with sensorineural hearing loss (SNHL) and identified pathogenic mutations in Epithelial Splicing-Regulatory Protein 1 (ESRP1). Patient-derived induced pluripotent stem cells showed alternative splicing defects that were restored upon repair of an ESRP1 mutant allele. To determine how ESRP1 mutations cause hearing loss, we evaluated Esrp1-/- mouse embryos and uncovered alterations in cochlear morphogenesis, auditory hair cell differentiation, and cell fate specification. Transcriptome analysis revealed impaired expression and splicing of genes with essential roles in cochlea development and auditory function. Aberrant splicing of Fgfr2 blocked stria vascularis formation due to erroneous ligand usage, which was corrected by reducing Fgf9 gene dosage. These findings implicate mutations in ESRP1 as a cause of SNHL and demonstrate the complex interplay between alternative splicing, inner ear development, and auditory function.

Alternative splicing occurs in the vast majority of human genes, giving rise to distinct mRNA and protein isoforms. We, and others, have previously identified hundreds of genes that change their isoform expression upon T cell activation via alternative splicing; however, how these changes link activation input with functional output remains largely unknown. Here, we investigate how costimulation of T cells through the CD28 receptor impacts alternative splicing in T cells activated through the T cell receptor (TCR, CD3) and find that while CD28 signaling alone has minimal impact on splicing, it enhances the extent of change for up to 20% of TCR-induced alternative splicing events. Interestingly, a set of CD28-enhanced splicing events occur within genes encoding key components of the apoptotic signaling pathway; namely caspase-9, Bax, and Bim. Using both CRISPR-edited cells and antisense oligos to force expression of specific isoforms, we show for all three of these genes that the isoform induced by CD3/CD28 costimulation promotes resistance to apoptosis, and that changes in all three genes together function combinatorially to further promote cell viability. Finally, we show that the JNK signaling pathway, induced downstream of CD3/CD28 costimulation, is required for each of these splicing events, further highlighting their co-regulation. Together, these findings demonstrate that alternative splicing is a key mechanism by which costimulation of CD28 promotes viability of activated T cells.

Cells adjust to hypoxic stress within the tumor microenvironment by downregulating energy-consuming processes including translation. To delineate mechanisms of cellular adaptation to hypoxia, we performed RNA-Seq of normoxic and hypoxic head and neck cancer cells. These data revealed a significant down regulation of genes known to regulate RNA processing and splicing. Exon-level analyses classified > 1,000 mRNAs as alternatively spliced under hypoxia and uncovered a unique retained intron (RI) in the master regulator of translation initiation, EIF2B5. Notably, this intron was expressed in solid tumors in a stage-dependent manner. We investigated the biological consequence of this RI and demonstrate that its inclusion creates a premature termination codon (PTC), that leads to a 65kDa truncated protein isoform that opposes full-length eIF2Bε to inhibit global translation. Furthermore, expression of 65kDa eIF2Bε led to increased survival of head and neck cancer cells under hypoxia, providing evidence that this isoform enables cells to adapt to conditions of low oxygen. Additional work to uncover -cis and -trans regulators of EIF2B5 splicing identified several factors that influence intron retention in EIF2B5: a weak splicing potential at the RI, hypoxia-induced expression and binding of the splicing factor SRSF3, and increased binding of total and phospho-Ser2 RNA polymerase II specifically at the intron retained under hypoxia. Altogether, these data reveal differential splicing as a previously uncharacterized mode of translational control under hypoxia and are supported by a model in which hypoxia-induced changes to cotranscriptional processing lead to selective retention of a PTC-containing intron in EIF2B5.

Human CD4+ T cells are often subdivided into distinct subtypes, including Th1, Th2, Th17, and Treg cells, that are thought to carry out distinct functions in the body. Typically, these T-cell subpopulations are defined by the expression of distinct gene repertoires; however, there is variability between studies regarding the methods used for isolation and the markers used to define each T-cell subtype. Therefore, how reliably studies can be compared to one another remains an open question. Moreover, previous analysis of gene expression in CD4+ T-cell subsets has largely focused on gene expression rather than alternative splicing. Here we take a meta-analysis approach, comparing eleven independent RNA-seq studies of human Th1, Th2, Th17, and/or Treg cells to determine the consistency in gene expression and splicing within each subtype across studies. We find that known master-regulators are consistently enriched in the appropriate subtype; however, cytokines and other genes often used as markers are more variable. Importantly, we also identify previously unknown transcriptomic markers that appear to consistently differentiate between subsets, including a few Treg-specific splicing patterns. Together this work highlights the heterogeneity in gene expression between samples designated as the same subtype, but also suggests additional markers that can be used to define functional groupings.

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.