Disentangling the splicing factor programs underlying complex molecular phenotypes.

The regulation of exon inclusion through alternative splicing tunes the cell's behavior by increasing the functional diversity of the transcriptome and the proteome. Splicing factors work in concert to generate gene isoform pools that contribute to cell phenotypes yet their activity is controlled by multiple regulatory and signaling layers. This hinders identification of functional, phenotype-specific splicing factors using traditional single-omic measurements, such as their mutational state or expression. To address this challenge, we propose repurposing the virtual inference of protein activity by enriched regulon analysis (VIPER) to measure splicing factor activity solely from their downstream exon transcriptomic inclusion signatures. This approach is effective in assessing the effect of co-occurring splicing factor perturbations, as well as their post-translational regulation. As proof of concept, we dissect recurrent splicing factor programs underlying tumorigenesis including aberrantly activated factors acting as oncogenes and inactivated ones acting as tumor suppressors, which are undetectable by more conventional methodologies. Activation and inactivation of these cancer splicing programs effectively stratifies overall survival, as well as cancer hallmarks such as proliferation and immune evasion. Altogether, repurposing network-based inference of protein activity for splicing factor networks distills common, functionally relevant splicing programs in otherwise heterogeneous molecular contexts.

Oncogenic NOVA1 expression dysregulates alternative splicing in breast cancer.

Neuro-Oncological Ventral Antigen 1 (NOVA1) is best known for its role in mediating an alternative splicing (AS) program in neurons, yet was first discovered as an antigen expressed in breast tumors, causing rare autoimmune reactions and paraneoplastic neurological disorders (PNDs). The PND model suggests a plausible role of the tumor antigen expression in tumor suppression, whereas it has emerged that NOVA may function as an oncogene in a variety of cancers. In addition, whether NOVA mediates AS in breast cancer remains unanswered. Here we examine the AS profiles of breast invasive carcinoma (BRCA) tumor samples and demonstrate that ectopic NOVA1 expression led to the activation of neuron-like splicing patterns in many genes, including exons targeted by NOVA in the brain. The splicing dysregulation is especially prevalent in cell periphery and cytoskeleton genes related to cell-cell communication, actin-based movement, and neuronal functions. We find that NOVA1-mediated AS is most prominent in Luminal A tumors and high NOVA1 expression in this subtype is associated with poorer prognosis. Our results suggest that ectopic NOVA1 in tumors has regulatory activity affecting pathways with high relevance to tumor progression and that this might be a more general mechanism for PND antigens.

Reverse engineering neuron type-specific and type-orthogonal splicing-regulatory networks using single-cell transcriptomes.

Cell type-specific alternative splicing (AS) enables differential gene isoform expression between diverse neuron types with distinct identities and functions. Current studies linking individual RNA-binding proteins (RBPs) to AS in a few neuron types underscore the need for holistic modeling. Here, we use network reverse engineering to derive a map of the neuron type-specific AS regulatory landscape from 133 mouse neocortical cell types defined by single-cell transcriptomes. This approach reliably inferred the regulons of 350 RBPs and their cell type-specific activities. Our analysis revealed driving factors delineating neuronal identities, among which we validated Elavl2 as a key RBP for MGE-specific splicing in GABAergic interneurons using an in vitro ESC differentiation system. We also identified a module of exons and candidate regulators specific for long- and short-projection neurons across multiple neuronal classes. This study provides a resource for elucidating splicing regulatory programs that drive neuronal molecular diversity, including those that do not align with gene expression-based classifications.

Complexity and graded regulation of neuronal cell-type-specific alternative splicing revealed by single-cell RNA sequencing.

The enormous cellular diversity in the mammalian brain, which is highly prototypical and organized in a hierarchical manner, is dictated by cell-type-specific gene-regulatory programs at the molecular level. Although prevalent in the brain, the contribution of alternative splicing (AS) to the molecular diversity across neuronal cell types is just starting to emerge. Here, we systematically investigated AS regulation across over 100 transcriptomically defined neuronal types of the adult mouse cortex using deep single-cell RNA-sequencing data. We found distinct splicing programs between glutamatergic and GABAergic neurons and between subclasses within each neuronal class. These programs consist of overlapping sets of alternative exons showing differential splicing at multiple hierarchical levels. Using an integrative approach, our analysis suggests that RNA-binding proteins (RBPs) Celf1/2, Mbnl2, and Khdrbs3 are preferentially expressed and more active in glutamatergic neurons, while Elavl2 and Qk are preferentially expressed and more active in GABAergic neurons. Importantly, these and additional RBPs also contribute to differential splicing between neuronal subclasses at multiple hierarchical levels, and some RBPs contribute to splicing dynamics that do not conform to the hierarchical structure defined by the transcriptional profiles. Thus, our results suggest graded regulation of AS across neuronal cell types, which may provide a molecular mechanism to specify neuronal identity and function that are orthogonal to established classifications based on transcriptional regulation.

Pharmacological Profiling of Purified Human Stem Cell-Derived and Primary Mouse Motor Neurons.

Directed differentiation of human pluripotent stem cells (hPSCs) has enabled the generation of specific neuronal subtypes that approximate the intended primary mammalian cells on both the RNA and protein levels. These cells offer unique opportunities, including insights into mechanistic understanding of the early driving events in neurodegenerative disease, replacement of degenerating cell populations, and compound identification and evaluation in the context of precision medicine. However, whether the derived neurons indeed recapitulate the physiological features of the desired bona fide neuronal subgroups remains an unanswered question and one important for validating stem cell models as accurate functional representations of the primary cell types. Here, we purified both hPSC-derived and primary mouse spinal motor neurons in parallel and used extracellular multi-electrode array (MEA) recording to compare the pharmacological sensitivity of neuronal excitability and network function. We observed similar effects for most receptor and channel agonists and antagonists, supporting the consistency between human PSC-derived and mouse primary spinal motor neuron models from a physiological perspective.

The Splicing Code Goes Deep.

The importance of genomic sequence context in generating transcriptome diversity through RNA splicing is independently unmasked by two studies in this issue (Jaganathan et al., 2019; Baeza-Centurion et al., 2019).

Transcriptome and in Vitro Differentiation Profile of Human Embryonic Stem Cell Derived NKX2.1-Positive Neural Progenitors.

The generation of inhibitory interneuron progenitors from human embryonic stem cells (ESCs) is of great interest due to their potential use in transplantation therapies designed to treat central nervous system disorders. The medial ganglionic eminence (MGE) is a transient embryonic structure in the ventral telencephalon that is a major source of cortical GABAergic inhibitory interneuron progenitors. These progenitors migrate tangentially to sites in the cortex and differentiate into a variety of interneuron subtypes, forming local synaptic connections with excitatory projection neurons to modulate activity of the cortical circuitry. The homeobox domain-containing transcription factor NKX2.1 is highly expressed in the MGE and pre-optic area of the ventral subpallium and is essential for specifying cortical interneuron fate. Using a combination of growth factor agonists and antagonists to specify ventral telencephalic fates, we previously optimized a protocol for the efficient generation of NKX2.1-positive MGE-like neural progenitors from human ESCs. To establish their identity, we now characterize the transcriptome of these MGE-like neural progenitors using RNA sequencing and demonstrate the capacity of these cells to differentiate into inhibitory interneurons in vitro using a neuron-astrocyte co-culture system. These data provide information on the potential origin of interneurons in the human brain.