# Deciphering splicing heterogeneity at single-cell resolution by SCSES

**Authors:** Xiao Wen, Xuan Lv, Dan Guo, Nan Han, Lei Zhou, Peizhuo Wang, Zhaoqi Liu

PMC · DOI: 10.1038/s41467-025-64517-5 · 2025-10-27

## TL;DR

SCSES is a new tool that improves the analysis of alternative splicing in single-cell RNA sequencing data, revealing hidden splicing patterns and cell subgroups.

## Contribution

SCSES introduces a novel computational framework for estimating splicing heterogeneity by sharing information across similar cells and events.

## Key findings

- SCSES outperforms existing methods in recovering splicing changes and diversity in simulated data.
- SCSES identifies splicing heterogeneity and cell subgroups with unique splicing patterns not detected by standard clustering.
- The tool is versatile and applicable across different species and sequencing platforms.

## Abstract

Alternative splicing (AS) plays a critical role in generating cellular transcriptomic heterogeneity. While single-cell RNA sequencing (scRNA-seq) has become a standard approach for exploring this heterogeneity, it remains challenging to accurately characterize splicing changes at the single-cell level due to high dropout rates, inevitable noise, and limited coverage. To address this, we developed SCSES (Single-Cell Splicing EStimation), a computational framework designed to enhance the AS profiles. SCSES infers and completes the missing splicing changes by sharing information across similar cells and events with data diffusion. Through systematic simulation studies, SCSES outperforms existing algorithms in recovering percent spliced-in (PSI) values and diversity across cell populations. When applied to various datasets, SCSES uncovers substantial splicing heterogeneity and cell subgroups with exclusive splicing patterns, which cannot be captured by conventional single-cell gene expression clustering. Together, our study provides SCSES as a valuable tool in deciphering splicing heterogeneity and is widely capable of handling different biological scenarios, species and sequencing platforms.

Alternative splicing (AS) is a source of transcriptomic heterogeneity, but analysing splicing changes at single-cell level is challenging due to dropout, noise, and limited coverage in scRNA-seq. Here, the authors propose SCSES, a tool to identify AS events and estimate splicing intensity at single-cell resolution. SCSES reveals high splicing heterogeneity beyond conventional expression-based clustering.

## Full-text entities

- **Genes:** Runx1 (runt related transcription factor 1) [NCBI Gene 12394] {aka AML1, CBF-alpha-2, Cbfa2, Pebp2a2, Pebpa2b}, VPS29 (VPS29 retromer complex component) [NCBI Gene 51699] {aka DC15, DC7, PEP11}, SOX17 (SRY-box transcription factor 17) [NCBI Gene 64321] {aka PPH7, VUR3}, TCF15 (transcription factor 15) [NCBI Gene 6939] {aka EC2, PARAXIS, bHLHa40}, MCM7 (minichromosome maintenance complex component 7) [NCBI Gene 4176] {aka CDC47, MCM2, P1.1-MCM3, P1CDC47, P85MCM, PNAS146}, EPS15 (epidermal growth factor receptor pathway substrate 15) [NCBI Gene 2060] {aka AF-1P, AF1P, MLLT5}, PIK3CB (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta) [NCBI Gene 5291] {aka P110BETA, PI3K, PI3KBETA, PIK3C1}, RBP4 (retinol binding protein 4) [NCBI Gene 5950] {aka MCOPCB10, RDCCAS}, TCF19 (transcription factor 19) [NCBI Gene 6941] {aka SC1, TCF-19}, Becn1 (beclin 1, autophagy related) [NCBI Gene 56208] {aka Atg6}, CD79A (CD79a molecule) [NCBI Gene 973] {aka IGA, IGAlpha, MB-1, MB1}, Lamp2 (lysosomal-associated membrane protein 2) [NCBI Gene 16784] {aka CD107b, LGP-B, Lamp II, Lamp-2, Lamp-2a, Lamp-2b}, EIF4B (eukaryotic translation initiation factor 4B) [NCBI Gene 1975] {aka EIF-4B, PRO1843}, USP8 (ubiquitin specific peptidase 8) [NCBI Gene 9101] {aka HumORF8, PITA4, SPG59, UBPY}, HUWE1 (HECT, UBA and WWE domain containing E3 ubiquitin protein ligase 1) [NCBI Gene 10075] {aka ARF-BP1, HECTH9, HSPC272, Ib772, LASU1, MRXST}, RBM24 (RNA binding motif protein 24) [NCBI Gene 221662] {aka RNPC6, dJ259A10.1}, Hspa8 (heat shock protein family A (Hsp70) member 8) [NCBI Gene 15481] {aka 2410008N15Rik, Hsc70, Hsc71, Hsc73, Hsp73, Hspa10}, NKS1 (natural killer cell susceptibility 1) [NCBI Gene 4819] {aka EC-1, EC1}, CASC3 (CASC3 exon junction complex subunit) [NCBI Gene 22794] {aka BTZ, MLN51}, Sh3glb1 (SH3-domain GRB2-like B1 (endophilin)) [NCBI Gene 54673] {aka Bif-1, mKIAA0491}, NUMB (NUMB endocytic adaptor protein) [NCBI Gene 8650] {aka C14orf41, S171, c14_5527}, SRRM2 (serine/arginine repetitive matrix 2) [NCBI Gene 23524] {aka 300-KD, CWF21, Cwc21, HSPC075, MRD72, SRL300}, CDC25B (cell division cycle 25B) [NCBI Gene 994] {aka MPIP2}, SRSF3 (serine and arginine rich splicing factor 3) [NCBI Gene 6428] {aka SFRS3, SRp20}, DNMT3B (DNA methyltransferase 3 beta) [NCBI Gene 1789] {aka FSHD4, ICF, ICF1, M.HsaIIIB}, LINC02605 (long intergenic non-protein coding RNA 2605) [NCBI Gene 112935892] {aka AS, IL-7, IL-7-AS}, ADAM10 (ADAM metallopeptidase domain 10) [NCBI Gene 102] {aka AD10, AD18, CD156c, CDw156, HsT18717, MADM}, Prtn3 (proteinase 3) [NCBI Gene 19152] {aka PR-3, PR3, mPR3}, ESRP1 (epithelial splicing regulatory protein 1) [NCBI Gene 54845] {aka DFNB109, RBM35A, RMB35A}, ERBB4 (erb-b2 receptor tyrosine kinase 4) [NCBI Gene 2066] {aka ALS19, HER4, p180erbB4}, POU5F1 (POU class 5 homeobox 1) [NCBI Gene 5460] {aka OCT3, OCT4, OCT4Borf1, OTF-3, OTF3, OTF4}, NOTCH2 (notch receptor 2) [NCBI Gene 4853] {aka AGS2, HJCYS, hN2}, NFE2L3 (NFE2 like bZIP transcription factor 3) [NCBI Gene 9603] {aka NRF3}, Ctsg (cathepsin G) [NCBI Gene 13035] {aka CatG, VSP}, TECR (trans-2,3-enoyl-CoA reductase) [NCBI Gene 9524] {aka GPSN2, MRT14, SC2, TER}, Elane (elastase, neutrophil expressed) [NCBI Gene 50701] {aka Ela2, F430011M15Rik, NE}, CADM1 (cell adhesion molecule 1) [NCBI Gene 23705] {aka BL2, IGSF4, IGSF4A, NECL2, Necl-2, RA175}, Uvrag (UV radiation resistance associated gene) [NCBI Gene 78610] {aka 9530039D02Rik, Uvrag1, Uvragl}, MRPL20 (mitochondrial ribosomal protein L20) [NCBI Gene 55052] {aka L20mt, MRP-L20, bL20m}, CHEK1 (checkpoint kinase 1) [NCBI Gene 1111] {aka CHK1, OZEMA21}, EDEM1 (ER degradation enhancing alpha-mannosidase like protein 1) [NCBI Gene 9695] {aka EDEM}, AURKA (aurora kinase A) [NCBI Gene 6790] {aka AIK, ARK1, AURA, BTAK, PPP1R47, STK15}, JMJD1C (jumonji domain containing 1C) [NCBI Gene 221037] {aka KDM3C, TRIP-8, TRIP8}, FOXO3 (forkhead box O3) [NCBI Gene 2309] {aka AF6q21, FKHRL1, FKHRL1P2, FOXO2, FOXO3A}, PTBP1 (polypyrimidine tract binding protein 1) [NCBI Gene 5725] {aka HNRNP-I, HNRNPI, HNRPI, PTB, PTB-1, PTB-T}, CXCR4 (C-X-C motif chemokine receptor 4) [NCBI Gene 7852] {aka CD184, D2S201E, FB22, HM89, HSY3RR, LCR1}
- **Diseases:** AS (MESH:C536589), Cancer (MESH:D009369), Ovarian Cancer (MESH:D010051), MM (MESH:D009101), ascites (MESH:D001201), DE (MESH:D003635), blood tumor (MESH:D009383), PD (MESH:D010300), BD (MESH:D021081), DSEs (MESH:D002318), metastasis (MESH:D009362), ND (MESH:C580335)
- **Chemicals:** Gln (MESH:D005973), thalidomide (MESH:D013792), bortezomib (MESH:D000069286)
- **Species:** Homo sapiens (human, species) [taxon 9606], Mus musculus (house mouse, species) [taxon 10090]
- **Mutations:** Lys650Glu
- **Cell lines:** HL-60 — Homo sapiens (Human), Adult acute myeloid leukemia with maturation, Cancer cell line (CVCL_0002), 758 — Homo sapiens (Human), Transformed cell line (CVCL_F569), ESCs — Mus musculus (Mouse), Embryonic stem cell (CVCL_9108), HCT116 — Homo sapiens (Human), Colon carcinoma, Cancer cell line (CVCL_0291), NiPSC — Homo sapiens (Human), Induced pluripotent stem cell (CVCL_C2VK), HCC1954 — Homo sapiens (Human), Breast ductal carcinoma, Cancer cell line (CVCL_1259), NMN-C2 — Homo sapiens (Human), Neuroblastoma, Cancer cell line (CVCL_0529), MN — Mus musculus (Mouse), Hybridoma (CVCL_G564), MN-C1 — Homo sapiens (Human), Rett syndrome, Finite cell line (CVCL_EQ52), Mono3 — Homo sapiens (Human), Adult acute myelomonocytic leukemia, Transformed cell line (CVCL_WN74), H9 — Homo sapiens (Human), Sezary syndrome, Cancer cell line (CVCL_1240), NMN-C1 — Pan troglodytes (Chimpanzee), Induced pluripotent stem cell (CVCL_1G30), SC4 — Mus musculus (Mouse), Hybridoma (CVCL_C4H4), SC3 — Homo sapiens (Human), Embryonic stem cell (CVCL_6F21), HepG2 — Homo sapiens (Human), Hepatoblastoma, Cancer cell line (CVCL_0027), H9 hESC — Homo sapiens (Human), Embryonic stem cell (CVCL_WN07), MM16 — Homo sapiens (Human), Pleural biphasic mesothelioma, Cancer cell line (CVCL_L789)

## Figures

6 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12559242/full.md

---
Source: https://tomesphere.com/paper/PMC12559242