Abstract:Many popular single-cell and spatial transcriptomics platforms exhibit 3′ bias, making it challenging to resolve all transcripts but potentially more feasible to resolve alternative polyadenylation (APA) events. Despite the development of several tools for identifying APA events in scRNA-seq data, a neutral benchmark is lacking, complicating the choice for biologists. We categorized existing APA analysis tools into three main classes, with the alignment-based class being the largest and we further divided this category into four subtypes. We compared the performance of methods from each algorithmic subtype in terms of site identification, quantification, and differential expression analysis across four single-cell and spatial transcriptomic datasets, using matched nanopore data as ground truth. No single method showed absolute superiority in all comparisons. Therefore, we selected representative methods (Sierra, scAPAtrap, and SCAPE) to deeply analyze the impact of different algorithmic choices on performance. SCAPE which is based on the distance estimation demonstrated less sensitivity to changes in read length and sequencing depth. It identified the most sites and achieved high recall but does not account for the impact of alternative splicing on site identification, leading to a loss in precision. Sierra that fits a coverage distribution is sensitive to changes in sequencing depth and identifies relatively fewer sites, but it considers the impact of junction reads on site identification and this results in relatively high precision. scAPAtrap combines peak calling and soft clipping, both of which are sensitive to sequencing depth. Moreover, soft clipping is particularly sensitive to read length, with increased read length leading to more false positive sites. Quantification consistency was affected by Cell Ranger versions and parameters, influencing downstream analysis but having less effect on differential expression between cell types. Each method has unique strengths. SCAPE is recommended for low-coverage data, scAPAtrap for moderate read lengths including intergenic sites, and Sierra for high-depth data with alternative splicing considerations. Filtering low-confidence sites, choosing appropriate mapping tools, and optimizing window size can improve performance.

Mining alternative splicing patterns in scRNA-seq data using scASfind

Survey of Programs Used to Detect Alternative Splicing Isoforms from Deep Sequencing Data In Silico

ASAS-EGB: A statistical framework for estimating allele-specific alternative splicing events using transcriptome data.

CASH: a constructing comprehensive splice site method for detecting alternative splicing events

SAW: a Method to Identify Splicing Events from RNA-Seq Data Based on Splicing Fingerprints.

Interrogations of single-cell RNA splicing landscapes with SCASL define new cell identities with physiological relevance

Prospects of Identifying Alternative Splicing Events from Single-Cell RNA Sequencing Data

Comparisons of computational methods for differential alternative splicing detection using RNA-seq in plant systems

MkcDBGAS: a reference-free approach to identify comprehensive alternative splicing events in a transcriptome

How to choose the optimal RNA-Seq library characteristics for alternative splicing analysis

Machine learning-optimized targeted detection of alternative splicing

Comprehensive Assessment of Isoform Detection Methods for Third-Generation Sequencing Data

Identification of Novel Splice Sites and Alternative Splicing Events in Drosophila Melanogaster Using RNA-seq Data

A review of alternative splicing event detection from RNA-Seq data

Identification Exon Skipping Events from High-Throughput RNA Sequencing Data

ASpedia-R: a package to retrieve junction-incorporating features and knowledge-based functions of human alternative splicing events

Guidelines for alternative polyadenylation identification tools using single-cell and spatial transcriptomics data

ASTool: An Easy-to-Use Tool to Accurately Identify Alternative Splicing Events from Plant RNA-Seq Data

scAPAtrap: identification and quantification of alternative polyadenylation sites from single-cell RNA-seq data

junctionCounts: comprehensive alternative splicing analysis and prediction of isoform-level impacts to the coding sequence

Computational Identification of Alternative Splicing in Next-Generation Sequencing Data