A framework for designing splice-junction experiments in deep 3' single-cell RNA sequencing

Alternative splicing is cell-type-specific and disease-relevant, but single-cell RNA sequencing is optimized for gene-level quantification, and how sequencing depth governs splice-junction recovery in the dominant 3'-biased chemistries has not been quantified, so experiments cannot be designed to a depth target and datasets cannot be interpreted against expected recovery. We generated deep single-cell RNA sequencing of human airway epithelial cells (12 samples, 102 to 803 million uniquely mapped reads per sample) and benchmarked junction recovery in five annotated samples against matched bulk RNA sequencing from a 190-donor cohort. Junction detection did not saturate with depth, and validity was governed by read support rather than annotation status. Recovering 90% of the observed full-depth junction complement required approximately 52,250 reads per cell, several-fold more than gene-level analysis. Quantifying splicing was more depth-limited than detecting it, and rare cell types were limited by cell number rather than depth. Against matched bulk, single-cell recall plateaued at a practical ceiling set by sample size and sparse sampling, and was independent of distance from the 3' end. These results establish a framework for designing splicing experiments in deep 3'-biased single-cell RNA sequencing, and provide an airway-epithelial dataset with matched bulk for benchmarking splicing methods.