flexiMAP: a regression-based method for discovering differential alternative polyadenylation events in standard RNA-seq data.

MOTIVATION We present flexiMAP (flexible Modeling of Alternative PolyAdenylation), a new beta-regression-based method implemented in R, for discovering differential alternative polyadenylation events in standard RNA-seq data. RESULTS We show, using both simulated and real data, that flexiMAP exhibits a good balance between specificity and sensitivity and compares favourably to existing methods, especially at low fold changes. In addition, the tests on simulated data reveal some hitherto unrecognised caveats of existing methods. Importantly, flexiMAP allows modeling of multiple known covariates that often confound the results of RNA-seq data analysis. AVAILABILITY The flexiMAP R package is available at: https://github.com/kszkop/flexiMAPScripts and data to reproduce the analysis in this paper are available at: https://doi.org/10.5281/zenodo.3689788. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.