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

RNA modulation in asthma: unraveling the role of splicing and non-coding RNAs in disease pathogenesis

OBJECTIVE

To synthesize the current understanding of RNA-based regulatory mechanisms, focusing on how RNA splicing and non-coding RNAs shape immune responses and airway remodeling in asthma, with the aim of exploring their potential as therapeutic targets for asthma treatment.

DATASOURCE

Recent advances and emerging research in molecular biology and immunology related to RNA splicing, non-coding RNAs (lncRNAs, circRNAs), and N6-methyladenosine (m6A) RNA methylation in asthma pathogenesis.

STUDY SELECTIONS

The review incorporates studies highlighting the roles of alternative RNA splicing, non-coding RNAs (lncRNAs and circRNAs), and RNA methylation (m6A) in regulating immune and inflammatory pathways involved in asthma.

RESULTS

RNA splicing events, non-coding RNAs, and m6A RNA methylation are critical in modulating immune dysregulation, airway remodeling, and inflammation in asthma. These mechanisms influence key inflammatory pathways, mRNA stability, and the overall immune landscape of the disease.

CONCLUSION

RNA splicing and non-coding RNAs represent promising areas of research for understanding asthma's immune pathology and hold potential as novel therapeutic targets for more effective treatment strategies.

LincRNA-ASAO promotes dental pulp repair through interacting with PTBP1 to increase ALPL alternative splicing

BACKGROUND

Alternative splicing not only expands the genetic encoding of genes but also determines cellular activities. This study aimed to elucidate the regulation mechanism and biological functions of lincRNA-ASAO in the process of odontogenesis-related genes alternative splicing mediated odontogenic differentiation of hDPSCs.

METHODS

RACE, RNA-seq, FISH and bioinformatics techniques were used to identify novel lincRNA-ASAO. ALP staining, alizarin red staining, qRT-PCR and western blot were used to identify the role of lincRNA-ASAO in regulating the odontoblast differentiation of hDPSCs. The binding protein PTBP1 of lincRNA-ASAO was screened by RNA-Pulldown, protein profiling and bioinformatics. The target gene ALPL of lincRNA-ASAO/PTBP1 was identified by RNA-seq, bioinformatics technology and DNA agarose gel electrophoresis. FISH, IF, PAR-CLIP and bioinformatics techniques were used to determine the roles of lincRNA-ASAO, PTBP1 and ALPL pre-mRNA in the odontoblast differentiation of hDPSCs.

RESULTS

We identified a novel lincRNA-ASAO that could promote the odontogenic differentiation of human Dental Pulp Stem Cells (hDPSCs). And, the interaction between lincRNA-ASAO and alternative splicing factor PTBP1 promoted the odontoblast differentiation of hDPSCs. In addition, lincRNA-ASAO forms duplexes with ALPL pre-mRNA, targeting PTBP1 to exonic splicing silencer (ESS) of ALPL and regulating exon 2 skipping. Notably, lincRNA-ASAO/PTBP1 regulated ALPL production to increase the type 2 splice variant, which promoted the odontoblast differentiation of hDPSCs.

CONCLUSIONS

We have identified the novel lincRNA-ASAO, which can promote the odontoblast differentiation of hDPSCs. The mechanism study found that lincRNA-ASAO/PTBP1 mediated the exon 2 skipping of ALPL pre-mRNA, resulting in the type 2 splice variant of ALPL. Our results enrich the understanding of lncRNAs and alternative splicing in regulating the odontoblast differentiation of hDPSCs, and provide clues to improve the clinical therapeutic potential of hDPSCs for dental pulp restoration.

A subcellular map of translational machinery composition and regulation at the single-molecule level

Millions of ribosomes are packed within mammalian cells, yet we lack tools to visualize them in toto and characterize their subcellular composition. In this study, we present ribosome expansion microscopy (RiboExM) to visualize individual ribosomes and an optogenetic proximity-labeling technique (ALIBi) to probe their composition. We generated a super-resolution ribosomal map, revealing subcellular translational hotspots and enrichment of 60S subunits near polysomes at the endoplasmic reticulum (ER). We found that Lsg1 tethers 60S to the ER and regulates translation of select proteins. Additionally, we discovered ribosome heterogeneity at mitochondria guiding translation of metabolism-related transcripts. Lastly, we visualized ribosomes in neurons, revealing a dynamic switch between monosomes and polysomes in neuronal translation. Together, these approaches enable exploration of ribosomal localization and composition at unprecedented resolution.

Exon inclusion signatures enable accurate estimation of splicing factor activity

Splicing factors control exon inclusion in messenger RNA, shaping transcriptome and proteome diversity. Their catalytic activity is regulated by multiple layers, making single-omic measurements on their own fall short in identifying which splicing factors underlie a phenotype. Here, we propose splicing factor activity can be estimated by interpreting changes in exon inclusion. We benchmark methods to construct splicing factor→exon networks and calculate activity. Combining RNA-seq perturbation-based networks with VIPER (virtual inference of protein activity by enriched regulon analysis) accurately captures splicing factor activation modulated by different regulatory layers. This approach consolidates splicing factor regulation into a single score derived solely from exon inclusion signatures, allowing functional interpretation of heterogeneous conditions. As a proof of concept, we identify recurrent cancer splicing programs, revealing oncogenic- and tumor suppressor-like splicing factors missed by conventional methods. These programs correlate with patient survival and key cancer hallmarks: initiation, proliferation, and immune evasion. Altogether, we show splicing factor activity can be accurately estimated from exon inclusion changes, enabling comprehensive analyses of splicing regulation with minimal data requirements.

Monocytes serve as Shiga toxin carriers during the development of hemolytic uremic syndrome

Shiga toxin (Stx)-induced hemolytic uremic syndrome (HUS) poses a life-threatening complication for which a definitive treatment remains elusive. To exert its cytotoxic effect on renal cells, Stx must be delivered from the infected intestines to the kidney. However, the mechanism underlying Stx delivery remains unclear. Here we pinpoint monocytes as the primary carriers responsible for transporting Stx2 to the renal region. Through single-cell sequencing analysis of Stx2-B-bound peripheral white blood cells sorted by flow cytometry, we observe that nearly all monocytes exhibit strong Stx2-B binding, whereas less than 10% of neutrophils are associated with Stx2-B, albeit with a lower affinity. Further examination of the single-cell dataset and cell binding assays suggest that monocytes likely bind to Stx2-B through the Toll-like receptor 4. Remarkably, Stx-laden monocytes demonstrate their ability to transport Stx2 to human renal glomerular endothelial cells (HRGEC), subsequently inducing apoptosis in HRGEC. In a mouse model of Stx1/2-positive EDL933 infection-induced HUS, the presence of Stx2-positive monocytes in peripheral blood and infiltrated kidney tissues was observed. Finally, depleting monocytes through the usage of a CD14 neutralizing antibody or blocking monocyte chemotaxis via inhibition of CCL2 notably mitigates kidney injury and dysfunction caused by lipopolysaccharide (LPS)/Stx2 treatment. Our findings unveil the pivotal role of monocytes in Stx delivery during STEC infection and offer a promising therapeutic approach for Stx-induced HUS.

mRNA Transcript Variants Expressed in Mammalian Cells

Gene expression or gene regulation studies often assume one gene expresses one mRNA. However, contrary to the conventional idea, a single gene in mammalian cells can express multiple transcript variants translated into several different proteins. The transcript variants are generated through transcription from alternative start sites and alternative post-transcriptional processing of the precursor mRNA (pre-mRNA). In addition, gene mutations and RNA editing further enhance the diversity of the transcript variants. The transcript variants can encode proteins with various domains, expanding the functional repertoire of a single gene. Some transcript variants may not encode proteins but function as non-coding RNAs and regulate gene expression. The expression level of the transcript variants may vary between cell types or within the same cells under different biological conditions. Transcript variants are characteristic of cell differentiation in a particular tissue, and the variants may play a key role in normal development and aging. Studies also reported that some transcript variants may have roles in disease pathogenesis. The biological significances urge studying the complexity of gene expression at the transcript level. This article updates the molecular basis of transcript variants in mammalian cells, including the formation mechanisms and potential roles in host biology. Gaining insight into the transcript variants will not only identify novel mechanisms of gene regulation but also unravel the role of the variants in health and disease.

sc-SPLASH provides ultra-efficient reference-free discovery in barcoded single-cell sequencing

Typical high-throughput single-cell RNA-sequencing (scRNA-seq) analyses are primarily conducted by (pseudo)alignment, through the lens of annotated gene models, and aimed at detecting differential gene expression. This misses diversity generated by other mechanisms that diversify the transcriptome such as splicing and V(D)J recombination, and is blind to sequences missing from imperfect reference genomes. Here, we present sc-SPLASH, a highly efficient pipeline that extends our SPLASH framework for statistics-first, reference-free discovery to barcoded scRNA-seq (10x Chromium) and spatial transcriptomics (10x Visium); we also provide its optimized module for preprocessing and k-mer counting in barcoded data, BKC, as a standalone tool. sc-SPLASH rediscovers known biology including V(D)J recombination and cell-type-specific alternative splicing in human and trans-splicing in tunicate (Ciona) and when applied to spatial datasets, detects sequence variation including tumor-specific somatic mutation. In sponge (Spongilla) and tunicate (Ciona), we uncover secreted repeat proteins expressed in immune-type cells and regulated during development; the sponge genes were absent from the reference assembly. sc-SPLASH provides a powerful alternative tool for exploring transcriptomes that is applicable to the breadth of life's diversity.

Statistical analysis supports pervasive RNA subcellular localization and alternative 3' UTR regulation

Targeted low-throughput studies have previously identified subcellular RNA localization as necessary for cellular functions including polarization, and translocation. Furthermore, these studies link localization to RNA isoform expression, especially 3' Untranslated Region (UTR) regulation. The recent introduction of genome-wide spatial transcriptomics techniques enables the potential to test if subcellular localization is regulated in situ pervasively. In order to do this, robust statistical measures of subcellular localization and alternative poly-adenylation (APA) at single-cell resolution are needed. Developing a new statistical framework called SPRAWL, we detect extensive cell-type specific subcellular RNA localization regulation in the mouse brain and to a lesser extent mouse liver. We integrated SPRAWL with a new approach to measure cell-type specific regulation of alternative 3' UTR processing and detected examples of significant correlations between 3' UTR length and subcellular localization. Included examples, Timp3, Slc32a1, Cxcl14, and Nxph1 have subcellular localization in the mouse brain highly correlated with regulated 3' UTR processing that includes the use of unannotated, but highly conserved, 3' ends. Together, SPRAWL provides a statistical framework to integrate multi-omic single-cell resolved measurements of gene-isoform pairs to prioritize an otherwise impossibly large list of candidate functional 3' UTRs for functional prediction and study. In these studies of data from mice, SPRAWL predicts that 3' UTR regulation of subcellular localization may be more pervasive than currently known.

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.

Spatial enrichment and genomic analyses reveal the link of NOMO1 with amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a severe motor neuron disease with uncertain genetic predisposition in most sporadic cases. The spatial architecture of cell types and gene expression are the basis of cell-cell interactions, biological function and disease pathology, but are not well investigated in the human motor cortex, a key ALS-relevant brain region. Recent studies indicated single nucleus transcriptomic features of motor neuron vulnerability in ALS motor cortex. However, the brain regional vulnerability of ALS-associated genes and the genetic link between region-specific genes and ALS risk remain largely unclear. Here, we developed an entropy-weighted differential gene expression matrix-based tool (SpatialE) to identify the spatial enrichment of gene sets in spatial transcriptomics. We benchmarked SpatialE against another enrichment tool (multimodal intersection analysis) using spatial transcriptomics data from both human and mouse brain tissues. To investigate regional vulnerability, we analysed three human motor cortex and two dorsolateral prefrontal cortex tissues for spatial enrichment of ALS-associated genes. We also used Cell2location to estimate the abundance of cell types in ALS-related cortex layers. To dissect the link of regionally expressed genes and ALS risk, we performed burden analyses of rare loss-of-function variants detected by whole-genome sequencing in ALS patients and controls, then analysed differential gene expression in the TargetALS RNA-sequencing dataset. SpatialE showed more accurate and specific spatial enrichment of regional cell type markers than multimodal intersection analysis in both mouse brain and human dorsolateral prefrontal cortex. Spatial transcriptomic analyses of human motor cortex showed heterogeneous cell types and spatial gene expression profiles. We found that 260 manually curated ALS-associated genes are significantly enriched in layer 5 of the motor cortex, with abundant expression of upper motor neurons and layer 5 excitatory neurons. Burden analyses of rare loss-of-function variants in Layer 5-associated genes nominated NOMO1 as a novel ALS-associated gene in a combined sample set of 6814 ALS patients and 3324 controls (P = 0.029). Gene expression analyses in CNS tissues revealed downregulation of NOMO1 in ALS, which is consistent with a loss-of-function disease mechanism. In conclusion, our integrated spatial transcriptomics and genomic analyses identified regional brain vulnerability in ALS and the association of a layer 5 gene (NOMO1) with ALS risk.

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.

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.

Bioinformatics integration reveals key genes associated with mitophagy in myocardial ischemia-reperfusion injury

BACKGROUND

Myocardial ischemia is a prevalent cardiovascular disorder associated with significant morbidity and mortality. While prompt restoration of blood flow is essential for improving patient outcomes, the subsequent reperfusion process can result in myocardial ischemia-reperfusion injury (MIRI). Mitophagy, a specialized autophagic mechanism, has consistently been implicated in various cardiovascular disorders. However, the specific connection between ischemia-reperfusion and mitophagy remains elusive. This study aims to elucidate and validate central mitophagy-related genes associated with MIRI through comprehensive bioinformatics analysis.

METHODS

We acquired the microarray expression profile dataset (GSE108940) from the Gene Expression Omnibus (GEO) and identified differentially expressed genes (DEGs) using GEO2R. Subsequently, these DEGs were cross-referenced with the mitophagy database, and differential nucleotide sequence analysis was performed through enrichment analysis. Protein-protein interaction (PPI) network analysis was employed to identify hub genes, followed by clustering of these hub genes using cytoHubba and MCODE within Cytoscape software. Gene set enrichment analysis (GSEA) was conducted on central genes. Additionally, Western blotting, immunofluorescence, and quantitative polymerase chain reaction (qPCR) analyses were conducted to validate the expression patterns of pivotal genes in MIRI rat model and H9C2 cardiomyocytes.

RESULTS

A total of 2719 DEGs and 61 mitophagy-DEGs were identified, followed by enrichment analyses and the construction of a PPI network. HSP90AA1, RPS27A, EEF2, EIF4A1, EIF2S1, HIF-1α, and BNIP3 emerged as the seven hub genes identified by cytoHubba and MCODE of Cytoscape software. Functional clustering analysis of HIF-1α and BNIP3 yielded a score of 9.647, as determined by Cytoscape (MCODE). In our MIRI rat model, Western blot and immunofluorescence analyses confirmed a significant elevation in the expression of HIF-1α and BNIP3, accompanied by a notable increase in the ratio of LC3II to LC3I. Subsequently, qPCR confirmed a significant upregulation of HIF-1α, BNIP3, and LC3 mRNA in the MIRI group. Activation of the HIF-1α/BNIP3 pathway mediates the regulation of the degree of Mitophagy, thereby effectively reducing apoptosis in rat H9C2 cardiomyocytes.

CONCLUSIONS

This study has identified seven central genes among mitophagy-related DEGs that may play a pivotal role in MIRI, suggesting a correlation between the HIF-1α/BNIP3 pathway of mitophagy and the pathogenesis of MIRI. The findings highlight the potential importance of mitophagy in MIRI and provide valuable insights into underlying mechanisms and potential therapeutic targets for further exploration in future studies.

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.