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.

Benchmarking alternative polyadenylation detection in single-cell and spatial transcriptomes

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

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

Accurate transcriptome-wide identification and quantification of alternative polyadenylation from RNA-seq data with APAIQ

Alternative Polyadenylation: Methods, Findings, and Impacts

Analysis of alternative polyadenylation from single-cell RNA-seq using scDaPars reveals cell subpopulations invisible to gene expression

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

Comprehensive Mapping of Alternative Polyadenylation Site Usage and Its Dynamics at Single-Cell Resolution

A comprehensive benchmarking of differential splicing tools for RNA-seq analysis at the event level

Distinct Polyadenylation Landscapes of Diverse Human Tissues Revealed by A Modified Pa-Seq Strategy

Using synthetic RNA to benchmark poly(A) length inference from direct RNA sequencing.

vizAPA: visualizing dynamics of alternative polyadenylation from bulk and single-cell data

SCAPE: a mixture model revealing single-cell polyadenylation diversity and cellular dynamics during cell differentiation and reprogramming

scAPAatlas: an atlas of alternative polyadenylation across cell types in human and mouse

Systematic comparative analysis of single cell RNA-sequencing methods

A comprehensive benchmarking with practical guidelines for cellular deconvolution of spatial transcriptomics

Using Population-Scale Transcriptomic and Genomic Data to Map 3′ UTR Alternative Polyadenylation Quantitative Trait Loci

Benchmarking spatial and single-cell transcriptomics integration methods for transcript distribution prediction and cell type deconvolution

Genome-wide Identification and Predictive Modeling of Tissue-Specific Alternative Polyadenylation.

Benchmarking algorithms for spatially variable gene identification in spatial transcriptomics

A multicenter study benchmarking single-cell RNA sequencing technologies using reference samples