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

Qian Zhao,Magnus Rattray
DOI: https://doi.org/10.1101/2024.11.29.626111
2024-12-04
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.
Bioinformatics
What problem does this paper attempt to address?