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.

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

scAPAmod: Profiling Alternative Polyadenylation Modalities in Single Cells from Single-Cell RNA-Seq Data

Scapadb: a Comprehensive Database of Alternative Polyadenylation at Single-Cell Resolution.

Discovery of alternative polyadenylation dynamics from single cell types

scDAPA: detection and visualization of dynamic alternative polyadenylation from single cell RNA-seq data.

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

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

Integrative Analysis of Gene Expression and Alternative Polyadenylation from Single-Cell RNA-seq Data.

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

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

APAtrap: Identification and Quantification of Alternative Polyadenylation Sites from RNA-seq Data.

Alternative Polyadenylation: Methods, Findings, and Impacts

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

Dynamics of alternative polyadenylation in single root cells of Arabidopsis thaliana

A Survey on Methods for Predicting Polyadenylation Sites from DNA Sequences, Bulk RNA-seq, and Single-cell RNA-seq

Learning Association for Single-Cell Transcriptomics by Integrating Profiling of Gene Expression and Alternative Polyadenylation

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

Alternative polyadenylation by sequential activation of distal and proximal PolyA sites

by PAS-Seq Complex and dynamic landscape of RNA polyadenylation revealed Material Supplemental

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

Complex and dynamic landscape of RNA polyadenylation revealed by PAS-Seq.