RNA splicing programs define tissue compartments and cell types at single-cell resolution

Neoantigens in cancer immunotherapy: focusing on alternative splicing

Alternative splicing (AS) functions as a crucial program in transcriptional modulation, leading to proteomic diversity and functional alterations of proteins. These splicing actions induce various neoantigens that hold prognostic significance and contribute to various aspects of cancer progression, including immune responses against cancer. The advent of immunotherapy has remarkably revolutionized tumor therapy. In this regard, AS-derived neoantigens are potent targets for cancer vaccines and chimeric antigen receptor (CAR) T cell therapies. In this review, we outline that AS-derived neoantigens serve as promising immunotherapeutic targets and guide immunotherapy strategies. This evidence contributes to a deeper comprehension of the complexity of proteomic diversity and provides novel perspectives and techniques for precision medicine in immunotherapy. Moreover, we underscore the obstacles that are awaited to be addressed for this novel approach to become clinically applicable.

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.

Regulation of endocrine cell alternative splicing revealed by single-cell RNA sequencing in type 2 diabetes pathogenesis

The prevalent RNA alternative splicing (AS) contributes to molecular diversity, which has been demonstrated in cellular function regulation and disease pathogenesis. However, the contribution of AS in pancreatic islets during diabetes progression remains unclear. Here, we reanalyze the full-length single-cell RNA sequencing data from the deposited database to investigate AS regulation across human pancreatic endocrine cell types in non-diabetic (ND) and type 2 diabetic (T2D) individuals. Our analysis demonstrates the significant association between transcriptomic AS profiles and cell-type-specificity, which could be applied to distinguish the clustering of major endocrine cell types. Moreover, AS profiles are enabled to clearly define the mature subset of β-cells in healthy controls, which is completely lost in T2D. Further analysis reveals that RNA-binding proteins (RBPs), heterogeneous nuclear ribonucleoproteins (hnRNPs) and FXR1 family proteins are predicted to induce the functional impairment of β-cells through regulating AS profiles. Finally, trajectory analysis of endocrine cells suggests the β-cell identity shift through dedifferentiation and transdifferentiation of β-cells during the progression of T2D. Together, our study provides a mechanism for regulating β-cell functions and suggests the significant contribution of AS program during diabetes pathogenesis.

Histamine H3 Receptor Isoforms: Insights from Alternative Splicing to Functional Complexity

Alternative splicing significantly enhances the diversity of the G protein-coupled receptor (GPCR) family, including the histamine H3 receptor (H3R). This post-transcriptional modification generates multiple H3R isoforms with potentially distinct pharmacological and physiological profiles. H3R is primarily involved in the presynaptic inhibition of neurotransmitter release in the central nervous system. Despite the approval of pitolisant for narcolepsy (Wakix®) and daytime sleepiness in adults with obstructive sleep apnea (Ozawade®) and ongoing clinical trials for other H3R antagonists/inverse agonists, the functional significance of the numerous H3R isoforms remains largely enigmatic. Recent publicly available RNA sequencing data have confirmed the expression of multiple H3R isoforms in the brain, with some isoforms exhibiting unique tissue-specific distribution patterns hinting at isoform-specific functions and interactions within neural circuits. In this review, we discuss the complexity of H3R isoforms with a focus on their potential roles in central nervous system (CNS) function. Comparative analysis across species highlights evolutionary conservation and divergence in H3R splicing, suggesting species-specific regulatory mechanisms. Understanding the functionality of H3R isoforms is crucial for the development of targeted therapeutics. This knowledge will inform the design of more precise pharmacological interventions, potentially enhancing therapeutic efficacy and reducing adverse effects in the treatment of neurological and psychiatric disorders.

RPS24 alternative splicing is a marker of cancer progression and epithelial-mesenchymal transition

Although alternative splicing (AS) is a major mechanism that adds diversity to gene expression patterns, its precise role in generating variability in ribosomal proteins, known as ribosomal heterogeneity, remains unclear. The ribosomal protein S24 (RPS24) gene, encoding a ribosomal component, undergoes AS; however, in-depth studies have been challenging because of three microexons between exons 4 and 6. We conducted a detailed analysis of RPS24 AS isoforms using a direct approach to investigate the splicing junctions related to these microexons, focusing on four AS isoforms. Each of these isoforms showed tissue specificity and relative differences in expression among cancer types. Significant differences in the proportions of these RPS24 AS isoforms between cancerous and normal tissues across diverse cancer types were also observed. Our study highlighted a significant correlation between the expression levels of a specific RPS24 AS isoform and the epithelial-mesenchymal transition process in lung and breast cancers. Our research contributes to a better understanding of the intricate regulatory mechanisms governing AS of ribosomal protein genes and highlights the biological implications of RPS24 AS isoforms in tissue development and tumorigenesis.

TM7SF3 controls TEAD1 splicing to prevent MASH-induced liver fibrosis

The mechanisms of hepatic stellate cell (HSC) activation and the development of liver fibrosis are not fully understood. Here, we show that deletion of a nuclear seven transmembrane protein, TM7SF3, accelerates HSC activation in liver organoids, primary human HSCs, and in vivo in metabolic-dysfunction-associated steatohepatitis (MASH) mice, leading to activation of the fibrogenic program and HSC proliferation. Thus, TM7SF3 knockdown promotes alternative splicing of the Hippo pathway transcription factor, TEAD1, by inhibiting the splicing factor heterogeneous nuclear ribonucleoprotein U (hnRNPU). This results in the exclusion of the inhibitory exon 5, generating a more active form of TEAD1 and triggering HSC activation. Furthermore, inhibiting TEAD1 alternative splicing with a specific antisense oligomer (ASO) deactivates HSCs in vitro and reduces MASH diet-induced liver fibrosis. In conclusion, by inhibiting TEAD1 alternative splicing, TM7SF3 plays a pivotal role in mitigating HSC activation and the progression of MASH-related fibrosis.

Alternative splicing in prostate cancer progression and therapeutic resistance

Prostate cancer (CaP) remains the second leading cause of cancer deaths in western men. CaP mortality results from diverse molecular mechanisms that mediate resistance to the standard of care treatments for metastatic disease. Recently, alternative splicing has been recognized as a hallmark of CaP aggressiveness. Alternative splicing events cause treatment resistance and aggressive CaP behavior and are determinants of the emergence of the two major types of late-stage treatment-resistant CaP, namely castration-resistant CaP (CRPC) and neuroendocrine CaP (NEPC). Here, we review recent multi-omics data that are uncovering the complicated landscape of alternative splicing events during CaP progression and the impact that different gene transcript isoforms can have on CaP cell biology and behavior. We discuss renewed insights in the molecular machinery by which alternative splicing occurs and contributes to the failure of systemic CaP therapies. The potential for alternative splicing events to serve as diagnostic markers and/or therapeutic targets is explored. We conclude by considering current challenges and promises associated with splicing-modulating therapies, and their potential for clinical translation into CaP patient care.

Counter-regulation of RNA stability by UPF1 and TDP43

RNA quality control is crucial for proper regulation of gene expression. Disruption of nonsense mediated mRNA decay (NMD), the primary RNA decay pathway responsible for the degradation of transcripts containing premature termination codons (PTCs), can disrupt development and lead to multiple diseases in humans and other animals. Similarly, therapies targeting NMD may have applications in hematological, neoplastic and neurological disorders. As such, tools capable of accurately quantifying NMD status could be invaluable for investigations of disease pathogenesis and biomarker identification. Toward this end, we assemble, validate, and apply a next-generation sequencing approach (NMDq) for identifying and measuring the abundance of PTC-containing transcripts. After validating NMDq performance and confirming its utility for tracking RNA surveillance, we apply it to determine pathway activity in two neurodegenerative diseases, amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) characterized by RNA misprocessing and abnormal RNA stability. Despite the genetic and pathologic evidence implicating dysfunctional RNA metabolism, and NMD in particular, in these conditions, we detected no significant differences in PTC-encoding transcripts in ALS models or disease. Contrary to expectations, overexpression of the master NMD regulator UPF1 had little effect on the clearance of transcripts with PTCs, but rather restored RNA homeostasis through differential use and decay of alternatively poly-adenylated isoforms. Together, these data suggest that canonical NMD is not a significant contributor to ALS/FTD pathogenesis, and that UPF1 promotes neuronal survival by regulating transcripts with abnormally long 3'UTRs.

Predicting the Structural Impact of Human Alternative Splicing

Protein structure prediction with neural networks is a powerful new method for linking protein sequence, structure, and function, but structures have generally been predicted for only a single isoform of each gene, neglecting splice variants. To investigate the structural implications of alternative splicing, we used AlphaFold2 to predict the structures of more than 11,000 human isoforms. We employed multiple metrics to identify splicing-induced structural alterations, including template matching score, secondary structure composition, surface charge distribution, radius of gyration, accessibility of post-translational modification sites, and structure-based function prediction. We identified examples of how alternative splicing induced clear changes in each of these properties. Structural similarity between isoforms largely correlated with degree of sequence identity, but we identified a subset of isoforms with low structural similarity despite high sequence similarity. Exon skipping and alternative last exons tended to increase the surface charge and radius of gyration. Splicing also buried or exposed numerous post-translational modification sites, most notably among the isoforms of BAX. Functional prediction nominated numerous functional differences among isoforms of the same gene, with loss of function compared to the reference predominating. Finally, we used single-cell RNA-seq data from the Tabula Sapiens to determine the cell types in which each structure is expressed. Our work represents an important resource for studying the structure and function of splice isoforms across the cell types of the human body.

SPLASH: A statistical, reference-free genomic algorithm unifies biological discovery

Today's genomics workflows typically require alignment to a reference sequence, which limits discovery. We introduce a unifying paradigm, SPLASH (Statistically Primary aLignment Agnostic Sequence Homing), which directly analyzes raw sequencing data, using a statistical test to detect a signature of regulation: sample-specific sequence variation. SPLASH detects many types of variation and can be efficiently run at scale. We show that SPLASH identifies complex mutation patterns in SARS-CoV-2, discovers regulated RNA isoforms at the single-cell level, detects the vast sequence diversity of adaptive immune receptors, and uncovers biology in non-model organisms undocumented in their reference genomes: geographic and seasonal variation and diatom association in eelgrass, an oceanic plant impacted by climate change, and tissue-specific transcripts in octopus. SPLASH is a unifying approach to genomic analysis that enables expansive discovery without metadata or references.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Cardiac splicing as a diagnostic and therapeutic target

Despite advances in therapeutics for heart failure and arrhythmias, a substantial proportion of patients with cardiomyopathy do not respond to interventions, indicating a need to identify novel modifiable myocardial pathobiology. Human genetic variation associated with severe forms of cardiomyopathy and arrhythmias has highlighted the crucial role of alternative splicing in myocardial health and disease, given that it determines which mature RNA transcripts drive the mechanical, structural, signalling and metabolic properties of the heart. In this Review, we discuss how the analysis of cardiac isoform expression has been facilitated by technical advances in multiomics and long-read and single-cell sequencing technologies. The resulting insights into the regulation of alternative splicing - including the identification of cardiac splice regulators as therapeutic targets and the development of a translational pipeline to evaluate splice modulators in human engineered heart tissue, animal models and clinical trials - provide a basis for improved diagnosis and therapy. Finally, we consider how the medical and scientific communities can benefit from facilitated acquisition and interpretation of splicing data towards improved clinical decision-making and patient care.

[WITHDRAWN] SPLASH: a statistical, reference-free genomic algorithm unifies biological discovery

The authors have withdrawn this manuscript due to a duplicate posting of manuscript number BIORXIV/2022/497555. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author. The correct preprint can be found at doi: https://doi.org/10.1101/2022.06.24.497555.

SPLASH: a statistical, reference-free genomic algorithm unifies biological discovery

Today's genomics workflows typically require alignment to a reference sequence, which limits discovery. We introduce a new unifying paradigm, SPLASH (Statistically Primary aLignment Agnostic Sequence Homing), an approach that directly analyzes raw sequencing data to detect a signature of regulation: sample-specific sequence variation. The approach, which includes a new statistical test, is computationally efficient and can be run at scale. SPLASH unifies detection of myriad forms of sequence variation. We demonstrate that SPLASH identifies complex mutation patterns in SARS-CoV-2 strains, discovers regulated RNA isoforms at the single cell level, documents the vast sequence diversity of adaptive immune receptors, and uncovers biology in non-model organisms undocumented in their reference genomes: geographic and seasonal variation and diatom association in eelgrass, an oceanic plant impacted by climate change, and tissue-specific transcripts in octopus. SPLASH is a new unifying approach to genomic analysis that enables an expansive scope of discovery without metadata or references.

Primary cilia as dynamic and diverse signalling hubs in development and disease

Primary cilia, antenna-like sensory organelles protruding from the surface of most vertebrate cell types, are essential for regulating signalling pathways during development and adult homeostasis. Mutations in genes affecting cilia cause an overlapping spectrum of >30 human diseases and syndromes, the ciliopathies. Given the immense structural and functional diversity of the mammalian cilia repertoire, there is a growing disconnect between patient genotype and associated phenotypes, with variable severity and expressivity characteristic of the ciliopathies as a group. Recent technological developments are rapidly advancing our understanding of the complex mechanisms that control biogenesis and function of primary cilia across a range of cell types and are starting to tackle this diversity. Here, we examine the structural and functional diversity of primary cilia, their dynamic regulation in different cellular and developmental contexts and their disruption in disease.

Evidence-Based Guide to Using Artificial Introns for Tissue-Specific Knockout in Mice

Up until recently, methods for generating floxed mice either conventionally or by CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-Cas9 (CRISPR-associated protein 9) editing have been technically challenging, expensive and error-prone, or time-consuming. To circumvent these issues, several labs have started successfully using a small artificial intron to conditionally knockout (KO) a gene of interest in mice. However, many other labs are having difficulty getting the technique to work. The key problem appears to be either a failure in achieving correct splicing after the introduction of the artificial intron into the gene or, just as crucial, insufficient functional KO of the gene's protein after Cre-induced removal of the intron's branchpoint. Presented here is a guide on how to choose an appropriate exon and where to place the recombinase-regulated artificial intron (rAI) in that exon to prevent disrupting normal gene splicing while maximizing mRNA degradation after recombinase treatment. The reasoning behind each step in the guide is also discussed. Following these recommendations should increase the success rate of this easy, new, and alternative technique for producing tissue-specific KO mice.

Approaches for single-cell RNA sequencing across tissues and cell types

Single-cell sequencing of RNA (scRNA-seq) has advanced our understanding of cellular heterogeneity and signaling in developmental biology and disease. A large number of complementary assays have been developed to profile transcriptomes of individual cells, also in combination with other readouts, such as chromatin accessibility or antibody-based analysis of protein surface markers. As scRNA-seq technologies are advancing fast, it is challenging to establish robust workflows and up-to-date protocols that are best suited to address the large range of research questions. Here, we review scRNA-seq techniques from mRNA end-counting to total RNA in relation to their specific features and outline the necessary sample preparation steps and quality control measures. Based on our experience in dealing with the continuously growing portfolio from the perspective of a central single-cell facility, we aim to provide guidance on how workflows can be best automatized and share our experience in coping with the continuous expansion of scRNA-seq techniques.

Nanoparticle Spikes Enhance Cellular Uptake via Regulating Myosin IIA Recruitment

Spike-like nanostructures are omnipresent in natural and artificial systems. Although biorecognition of nanostructures to cellular receptors has been indicated as the primary factor for virus infection pathways, how the spiky morphology of DNA-modified nanoparticles affects their cellular uptake and intracellular fate remains to be explored. Here, we design dually emissive gold nanoparticles with varied spikiness (from 0 to 2) to probe the interactions of spiky nanoparticles with cells. We discovered that nanospikes at the nanoparticle regulated myosin IIA recruitment at the cell membrane during cellular uptake, thereby enhancing cellular uptake efficiency, as revealed by dual-modality (plasmonic and fluorescence) imaging. Furthermore, the spiky nanoparticles also exhibited facilitated endocytosis dynamics, as revealed by real-time dark-field microscopy (DFM) imaging and colorimetry-based classification algorithms. These findings highlight the crucial role of the spiky morphology in regulating the intracellular fate of nanoparticles, which may shed light on engineering theranostic nanocarriers.

Regulation of adult stem cell quiescence and its functions in the maintenance of tissue integrity

Adult stem cells are important for mammalian tissues, where they act as a cell reserve that supports normal tissue turnover and can mount a regenerative response following acute injuries. Quiescent stem cells are well established in certain tissues, such as skeletal muscle, brain, and bone marrow. The quiescent state is actively controlled and is essential for long-term maintenance of stem cell pools. In this Review, we discuss the importance of maintaining a functional pool of quiescent adult stem cells, including haematopoietic stem cells, skeletal muscle stem cells, neural stem cells, hair follicle stem cells, and mesenchymal stem cells such as fibro-adipogenic progenitors, to ensure tissue maintenance and repair. We discuss the molecular mechanisms that regulate the entry into, maintenance of, and exit from the quiescent state in mice. Recent studies revealed that quiescent stem cells have a discordance between RNA and protein levels, indicating the importance of post-transcriptional mechanisms, such as alternative polyadenylation, alternative splicing, and translation repression, in the control of stem cell quiescence. Understanding how these mechanisms guide stem cell function during homeostasis and regeneration has important implications for regenerative medicine.

Analysis of RNA processing directly from spatial transcriptomics data reveals previously unknown regulation

Technical advances have led to an explosion in the amount of biological data available in recent years, especially in the field of RNA sequencing. Specifically, spatial transcriptomics (ST) datasets, which allow each RNA molecule to be mapped to the 2D location it originated from within a tissue, have become readily available. Due to computational challenges, ST data has rarely been used to study RNA processing such as splicing or differential UTR usage. We apply the ReadZS and the SpliZ, methods developed to analyze RNA process in scRNA-seq data, to analyze spatial localization of RNA processing directly from ST data for the first time. Using Moran's I metric for spatial autocorrelation, we identify genes with spatially regulated RNA processing in the mouse brain and kidney, re-discovering known spatial regulation in Myl6 and identifying previously-unknown spatial regulation in genes such as Rps24, Gng13, Slc8a1, Gpm6a, Gpx3, ActB, Rps8, and S100A9. The rich set of discoveries made here from commonly used reference datasets provides a small taste of what can be learned by applying this technique more broadly to the large quantity of Visium data currently being created.

Decision level integration of unimodal and multimodal single cell data with scTriangulate

Decisively delineating cell identities from uni- and multimodal single-cell datasets is complicated by diverse modalities, clustering methods, and reference atlases. We describe scTriangulate, a computational framework to mix-and-match multiple clustering results, modalities, associated algorithms, and resolutions to achieve an optimal solution. Rather than ensemble approaches which select the "consensus", scTriangulate picks the most stable solution through coalitional iteration. When evaluated on diverse multimodal technologies, scTriangulate outperforms alternative approaches to identify high-confidence cell-populations and modality-specific subtypes. Unlike existing integration strategies that rely on modality-specific joint embedding or geometric graphs, scTriangulate makes no assumption about the distributions of raw underlying values. As a result, this approach can solve unprecedented integration challenges, including the ability to automate reference cell-atlas construction, resolve clonal architecture within molecularly defined cell-populations and subdivide clusters to discover splicing-defined disease subtypes. scTriangulate is a flexible strategy for unified integration of single-cell or multimodal clustering solutions, from nearly unlimited sources.

The Tabula Sapiens: A multiple-organ, single-cell transcriptomic atlas of humans

Molecular characterization of cell types using single-cell transcriptome sequencing is revolutionizing cell biology and enabling new insights into the physiology of human organs. We created a human reference atlas comprising nearly 500,000 cells from 24 different tissues and organs, many from the same donor. This atlas enabled molecular characterization of more than 400 cell types, their distribution across tissues, and tissue-specific variation in gene expression. Using multiple tissues from a single donor enabled identification of the clonal distribution of T cells between tissues, identification of the tissue-specific mutation rate in B cells, and analysis of the cell cycle state and proliferative potential of shared cell types across tissues. Cell type-specific RNA splicing was discovered and analyzed across tissues within an individual.

The SpliZ generalizes 'percent spliced in' to reveal regulated splicing at single-cell resolution

Detecting single-cell-regulated splicing from droplet-based technologies is challenging. Here, we introduce the splicing Z score (SpliZ), an annotation-free statistical method to detect regulated splicing in single-cell RNA sequencing. We applied the SpliZ to human lung cells, discovering hundreds of genes with cell-type-specific splicing patterns including ones with potential implications for basic and translational biology.

RNA splicing programs define tissue compartments and cell types at single-cell resolution

The extent splicing is regulated at single-cell resolution has remained controversial due to both available data and methods to interpret it. We apply the SpliZ, a new statistical approach, to detect cell-type-specific splicing in >110K cells from 12 human tissues. Using 10X Chromium data for discovery, 9.1% of genes with computable SpliZ scores are cell-type-specifically spliced, including ubiquitously expressed genes MYL6 and RPS24. These results are validated with RNA FISH, single-cell PCR, and Smart-seq2. SpliZ analysis reveals 170 genes with regulated splicing during human spermatogenesis, including examples conserved in mouse and mouse lemur. The SpliZ allows model-based identification of subpopulations indistinguishable based on gene expression, illustrated by subpopulation-specific splicing of classical monocytes involving an ultraconserved exon in SAT1. Together, this analysis of differential splicing across multiple organs establishes that splicing is regulated cell-type-specifically.