Efficient and Accurate Quantitative Profiling of Alternative Splicing Patterns of Any Complexity on a Laptop.

Alternative splicing (AS) is a widespread process underlying the generation of transcriptomic and proteomic diversity and is frequently misregulated in human disease. Accordingly, an important goal of biomedical research is the development of tools capable of comprehensively, accurately, and efficiently profiling AS. Here, we describe Whippet, an easy-to-use RNA-seq analysis method that rapidly-with hardware requirements compatible with a laptop-models and quantifies AS events of any complexity without loss of accuracy. Using an entropic measure of splicing complexity, Whippet reveals that one-third of human protein coding genes produce transcripts with complex AS events involving co-expression of two or more principal splice isoforms. We observe that high-entropy AS events are more prevalent in tumor relative to matched normal tissues and correlate with increased expression of proto-oncogenic splicing factors. Whippet thus affords the rapid and accurate analysis of AS events of any complexity, and as such will facilitate future biomedical research.

Misregulation of an Activity-Dependent Splicing Network as a Common Mechanism Underlying Autism Spectrum Disorders.

A key challenge in understanding and ultimately treating autism is to identify common molecular mechanisms underlying this genetically heterogeneous disorder. Transcriptomic profiling of autistic brains has revealed correlated misregulation of the neuronal splicing regulator nSR100/SRRM4 and its target microexon splicing program in more than one-third of analyzed individuals. To investigate whether nSR100 misregulation is causally linked to autism, we generated mutant mice with reduced levels of this protein and its target splicing program. Remarkably, these mice display multiple autistic-like features, including altered social behaviors, synaptic density, and signaling. Moreover, increased neuronal activity, which is often associated with autism, results in a rapid decrease in nSR100 and splicing of microexons that significantly overlap those misregulated in autistic brains. Collectively, our results provide evidence that misregulation of an nSR100-dependent splicing network controlled by changes in neuronal activity is causally linked to a substantial fraction of autism cases.

A global regulatory mechanism for activating an exon network required for neurogenesis.

The vertebrate and neural-specific Ser/Arg (SR)-related protein nSR100/SRRM4 regulates an extensive program of alternative splicing with critical roles in nervous system development. However, the mechanism by which nSR100 controls its target exons is poorly understood. We demonstrate that nSR100-dependent neural exons are associated with a unique configuration of intronic cis-elements that promote rapid switch-like regulation during neurogenesis. A key feature of this configuration is the insertion of specialized intronic enhancers between polypyrimidine tracts and acceptor sites that bind nSR100 to potently activate exon inclusion in neural cells while weakening 3' splice site recognition and contributing to exon skipping in nonneural cells. nSR100 further operates by forming multiple interactions with early spliceosome components bound proximal to 3' splice sites. These multifaceted interactions achieve dominance over neural exon silencing mediated by the splicing regulator PTBP1. The results thus illuminate a widespread mechanism by which a critical neural exon network is activated during neurogenesis.

An atlas of alternative splicing profiles and functional associations reveals new regulatory programs and genes that simultaneously express multiple major isoforms.

Alternative splicing (AS) generates remarkable regulatory and proteomic complexity in metazoans. However, the functions of most AS events are not known, and programs of regulated splicing remain to be identified. To address these challenges, we describe the Vertebrate Alternative Splicing and Transcription Database (VastDB), the largest resource of genome-wide, quantitative profiles of AS events assembled to date. VastDB provides readily accessible quantitative information on the inclusion levels and functional associations of AS events detected in RNA-seq data from diverse vertebrate cell and tissue types, as well as developmental stages. The VastDB profiles reveal extensive new intergenic and intragenic regulatory relationships among different classes of AS and previously unknown and conserved landscapes of tissue-regulated exons. Contrary to recent reports concluding that nearly all human genes express a single major isoform, VastDB provides evidence that at least 48% of multiexonic protein-coding genes express multiple splice variants that are highly regulated in a cell/tissue-specific manner, and that >18% of genes simultaneously express multiple major isoforms across diverse cell and tissue types. Isoforms encoded by the latter set of genes are generally coexpressed in the same cells and are often engaged by translating ribosomes. Moreover, they are encoded by genes that are significantly enriched in functions associated with transcriptional control, implying they may have an important and wide-ranging role in controlling cellular activities. VastDB thus provides an unprecedented resource for investigations of AS function and regulation.

Frac-seq reveals isoform-specific recruitment to polyribosomes.

Pre-mRNA splicing is required for the accurate expression of virtually all human protein coding genes. However, splicing also plays important roles in coordinating subsequent steps of pre-mRNA processing such as polyadenylation and mRNA export. Here, we test the hypothesis that nuclear pre-mRNA processing influences the polyribosome association of alternative mRNA isoforms. By comparing isoform ratios in cytoplasmic and polyribosomal extracts, we determined that the alternative products of ∼30% (597/1954) of mRNA processing events are differentially partitioned between these subcellular fractions. Many of the events exhibiting isoform-specific polyribosome association are highly conserved across mammalian genomes, underscoring their possible biological importance. We find that differences in polyribosome association may be explained, at least in part by the observation that alternative splicing alters the cis-regulatory landscape of mRNAs isoforms. For example, inclusion or exclusion of upstream open reading frames (uORFs) in the 5'UTR as well as Alu-elements and microRNA target sites in the 3'UTR have a strong influence on polyribosome association of alternative mRNA isoforms. Taken together, our data demonstrate for the first time the potential link between alternative splicing and translational control of the resultant mRNA isoforms.

Loss of exon identity is a common mechanism of human inherited disease.

It is widely accepted that at least 10% of all mutations causing human inherited disease disrupt splice-site consensus sequences. In contrast to splice-site mutations, the role of auxiliary cis-acting elements such as exonic splicing enhancers (ESE) and exonic splicing silencers (ESS) in human inherited disease is still poorly understood. Here we use a top-down approach to determine rates of loss or gain of known human exonic splicing regulatory (ESR) sequences associated with either disease-causing mutations or putatively neutral single nucleotide polymorphisms (SNPs). We observe significant enrichment toward loss of ESEs and gain of ESSs among inherited disease-causing variants relative to neutral polymorphisms, indicating that exon skipping may play a prominent role in aberrant gene regulation. Both computational and biochemical approaches underscore the relevance of exonic splicing enhancer loss and silencer gain in inherited disease. Additionally, we provide direct evidence that both SRp20 (SRSF3) and possibly PTB (PTBP1) are involved in the function of a splicing silencer that is created de novo by a total of 83 different inherited disease mutations in 67 different disease genes. Taken together, we find that ~25% (7154/27,681) of known mis-sense and nonsense disease-causing mutations alter functional splicing signals within exons, suggesting a much more widespread role for aberrant mRNA processing in causing human inherited disease than has hitherto been appreciated.

Accurate long-read transcript discovery and quantification at single-cell, pseudo-bulk and bulk resolution with Isosceles.

Accurate detection and quantification of mRNA isoforms from nanopore long-read sequencing remains challenged by technical noise, particularly in single cells. To address this, we introduce Isosceles, a computational toolkit that outperforms other methods in isoform detection sensitivity and quantification accuracy across single-cell, pseudo-bulk and bulk resolution levels, as demonstrated using synthetic and biologically-derived datasets. Here we show Isosceles improves the fidelity of single-cell transcriptome quantification at the isoform-level, and enables flexible downstream analysis. As a case study, we apply Isosceles, uncovering coordinated splicing within and between neuronal differentiation lineages. Isosceles is suitable to be applied in diverse biological systems, facilitating studies of cellular heterogeneity across biomedical research applications.

Single-cell long-read targeted sequencing reveals transcriptional variation in ovarian cancer.

Single-cell RNA sequencing predominantly employs short-read sequencing to characterize cell types, states and dynamics; however, it is inadequate for comprehensive characterization of RNA isoforms. Long-read sequencing technologies enable single-cell RNA isoform detection but are hampered by lower throughput and unintended sequencing of artifacts. Here we develop Single-cell Targeted Isoform Long-Read Sequencing (scTaILoR-seq), a hybridization capture method which targets over a thousand genes of interest, improving the median number of on-target transcripts per cell by 29-fold. We use scTaILoR-seq to identify and quantify RNA isoforms from ovarian cancer cell lines and primary tumors, yielding 10,796 single-cell transcriptomes. Using long-read variant calling we reveal associations of expressed single nucleotide variants (SNVs) with alternative transcript structures. Phasing of SNVs across transcripts enables the measurement of allelic imbalance within distinct cell populations. Overall, scTaILoR-seq is a long-read targeted RNA sequencing method and analytical framework for exploring transcriptional variation at single-cell resolution.

Differential contribution of transcriptomic regulatory layers in the definition of neuronal identity.

Previous transcriptomic profiling studies have typically focused on separately analyzing mRNA expression, alternative splicing and alternative polyadenylation differences between cell and tissue types. However, the relative contribution of these three transcriptomic regulatory layers to cell type specification is poorly understood. This question is particularly relevant to neurons, given their extensive heterogeneity associated with brain location, morphology and function. In the present study, we generated profiles for the three regulatory layers from developmentally and regionally distinct subpopulations of neurons from the mouse hippocampus and broader nervous system. Multi-omics factor analyses revealed differing contributions of each transcriptomic layer in the discrimination of neurons based on their stage of development, region, and function. Importantly, profiles of differential alternative splicing and polyadenylation better discriminated specific neuronal subtype populations than gene expression patterns. These results provide evidence for differential relative contributions of coordinated gene regulatory layers in the specification of neuronal subtypes.

The ribosome-engaged landscape of alternative splicing.

High-throughput RNA sequencing (RNA-seq) has revealed an enormous complexity of alternative splicing (AS) across diverse cell and tissue types. However, it is currently unknown to what extent repertoires of splice-variant transcripts are translated into protein products. Here, we surveyed AS events engaged by the ribosome. Notably, at least 75% of human exon-skipping events detected in transcripts with medium-to-high abundance in RNA-seq data were also detected in ribosome profiling data. Furthermore, relatively small subsets of functionally related splice variants are engaged by ribosomes at levels that do not reflect their absolute abundance, thus indicating a role for AS in modulating translational output. This mode of regulation is associated with control of the mammalian cell cycle. Our results thus suggest that a major fraction of splice variants is translated and that specific cellular functions including cell-cycle control are subject to AS-dependent modulation of translation output.

RNA-binding protein IGF2BP3 targeting of oncogenic transcripts promotes hematopoietic progenitor proliferation.

Posttranscriptional control of gene expression is important for defining both normal and pathological cellular phenotypes. In vitro, RNA-binding proteins (RBPs) have recently been shown to play important roles in posttranscriptional regulation; however, the contribution of RBPs to cell specification is not well understood. Here, we determined that the RBP insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3) is specifically overexpressed in mixed lineage leukemia-rearranged (MLL-rearranged) B-acute lymphoblastic leukemia (B-ALL), which constitutes a subtype of this malignancy associated with poor prognosis and high risk of relapse. IGF2BP3 was required for the survival of B-ALL cell lines, as knockdown led to decreased proliferation and increased apoptosis. Enforced expression of IGF2BP3 provided murine BM cells with a strong survival advantage, led to proliferation of hematopoietic stem and progenitor cells, and skewed hematopoietic development to the B cell/myeloid lineage. Cross-link immunoprecipitation and high throughput sequencing uncovered the IGF2BP3-regulated transcriptome, which includes oncogenes MYC and CDK6 as direct targets. IGF2BP3 regulated transcripts via targeting elements within 3' untranslated regions (3'UTR), and enforced IGF2BP3 expression in mice resulted in enhanced expression of Myc and Cdk6 in BM. Together, our data suggest that IGF2BP3-mediated targeting of oncogenic transcripts may represent a critical pathogenetic mechanism in MLL-rearranged B-ALL and support IGF2BP3 and its cognate RNA-binding partners as potential therapeutic targets in this disease.

IGF2BP3 Modulates the Interaction of Invasion-Associated Transcripts with RISC.

Insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3) expression correlates with malignancy, but its role(s) in pathogenesis remains enigmatic. We interrogated the IGF2BP3-RNA interaction network in pancreatic ductal adenocarcinoma (PDAC) cells. Using a combination of genome-wide approaches, we have identified 164 direct mRNA targets of IGF2BP3. These transcripts encode proteins enriched for functions such as cell migration, proliferation, and adhesion. Loss of IGF2BP3 reduced PDAC cell invasiveness and remodeled focal adhesion junctions. Individual nucleotide resolution crosslinking immunoprecipitation (iCLIP) revealed significant overlap of IGF2BP3 and microRNA (miRNA) binding sites. IGF2BP3 promotes association of the RNA-induced silencing complex (RISC) with specific transcripts. Our results show that IGF2BP3 influences a malignancy-associated RNA regulon by modulating miRNA-mRNA interactions.

Exon identity crisis: disease-causing mutations that disrupt the splicing code.

Cis-acting RNA elements control the accurate expression of human multi-exon protein coding genes. Single nucleotide variants altering the fidelity of this regulatory code and, consequently, pre-mRNA splicing are expected to contribute to the etiology of numerous human diseases.

MutPred Splice: machine learning-based prediction of exonic variants that disrupt splicing.

We have developed a novel machine-learning approach, MutPred Splice, for the identification of coding region substitutions that disrupt pre-mRNA splicing. Applying MutPred Splice to human disease-causing exonic mutations suggests that 16% of mutations causing inherited disease and 10 to 14% of somatic mutations in cancer may disrupt pre-mRNA splicing. For inherited disease, the main mechanism responsible for the splicing defect is splice site loss, whereas for cancer the predominant mechanism of splicing disruption is predicted to be exon skipping via loss of exonic splicing enhancers or gain of exonic splicing silencer elements. MutPred Splice is available at http://mutdb.org/mutpredsplice.

Structural RNAs of known and unknown function identified in malaria parasites by comparative genomics and RNA analysis.

As the genomes of more eukaryotic pathogens are sequenced, understanding how molecular differences between parasite and host might be exploited to provide new therapies has become a major focus. Central to cell function are RNA-containing complexes involved in gene expression, such as the ribosome, the spliceosome, snoRNAs, RNase P, and telomerase, among others. In this article we identify by comparative genomics and validate by RNA analysis numerous previously unknown structural RNAs encoded by the Plasmodium falciparum genome, including the telomerase RNA, U3, 31 snoRNAs, as well as previously predicted spliceosomal snRNAs, SRP RNA, MRP RNA, and RNAse P RNA. Furthermore, we identify six new RNA coding genes of unknown function. To investigate the relationships of the RNA coding genes to other genomic features in related parasites, we developed a genome browser for P. falciparum (http://areslab.ucsc.edu/cgi-bin/hgGateway). Additional experiments provide evidence supporting the prediction that snoRNAs guide methylation of a specific position on U4 snRNA, as well as predicting an snRNA promoter element particular to Plasmodium sp. These findings should allow detailed structural comparisons between the RNA components of the gene expression machinery of the parasite and its vertebrate hosts.