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.

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

Genome-wide Identification and Predictive Modeling of Polyadenylation Sites in Eukaryotes

Alternative Polyadenylation: Methods, Findings, and Impacts

Discovery of alternative polyadenylation dynamics from single cell types

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

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

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

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

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

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

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

Benchmarking alternative polyadenylation detection in single-cell and spatial transcriptomes

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

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

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

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

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

A Classification-Based Prediction Model of Messenger RNA Polyadenylation Sites

Alternative polyadenylation by sequential activation of distal and proximal PolyA sites

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

PolyAMiner-Bulk is a deep learning-based algorithm that decodes alternative polyadenylation dynamics from bulk RNA-seq data