Joint representation of spliced and unspliced RNA expands information content in single-cell expression data

MotivationSingle-cell RNA-seq captures both mature (spliced) and nascent (unspliced) transcripts, yet standard preprocessing typically collapses these signals into a single expression matrix, obscuring recoverable transcriptional structure and introducing splicing-dependent artifacts in downstream embeddings in foundation models. We asked whether treating spliced and unspliced counts as distinct feature sets increases recoverable latent structure in unsupervised dimension reduction. ResultsAcross four datasets with paired spliced and unspliced count layers, we applied non-negative matrix factorization (NMF) to spliced-only, unspliced-only, summed, and concatenated representations. Using cross-validation, we estimated the revealable rank k* (largest stable rank achieving near-optimal held-out error). In all datasets, concatenated representations consistently supported higher k* than spliced-only or unspliced-only inputs and typically exceeded the summed representation, indicating that preserving splicing status exposes additional stable latent programs. Minimal overlap between spliced- and unspliced-derived factors corresponded to modality-specific gene set enrichment patterns consistent with complementary steady-state and nascent regulatory activity. These results establish splicing status as an informative axis of variation and motivate preserving spliced and unspliced layers as first-class features to maximize information recovery in single-cell analysis. Availability and ImplementationThe singlet R package is available at github.com/zdebruine/singlet. Code to reproduce all figures is available at github.com/Arya-86/nmf-spliced-unspliced. Supplementary informationSupplementary data are appended to this manuscript.