Abstract
Intronic polyadenylation (IPA) is a key mechanism driving transcriptome diversity, yet its detection and functional characterization remain challenging owing to complex splicing patterns and the complexity of intronic regions. Here, we introduce IPAseek, a dynamic programming based computational framework that leverages the pruned exact linear time (PELT) algorithm and changepoints over a range of penalties (CROPS) to enable de novo identification of IPA events from bulk RNA-seq data. IPAseek robustly detects both composite and skipped IPA isoforms. Applying IPAseek to bulk RNA-seq of hematopoietic cell types reveals lineage and stage-specific IPA signatures, with lymphoid cells exhibiting higher IPA site usage compared with myeloid cells. Temporal profiling during megakaryocyte differentiation uncovers dynamic, gene-specific IPA regulation linked to functional pathways including peroxisomal metabolism and autophagy, which are known to play a crucial role in megakaryocytic differentiation, impacting the development and maturation of megakaryocytes. Further, integrative analysis demonstrates that IPA site usage is associated with lower DNA methylation within introns, supporting a regulatory axis connecting epigenetic state and IPA. This finding aligns with emerging evidence that DNA methylation modulates alternative polyadenylation via CTCF-mediated chromatin looping. Thus, IPAseek provides a platform to characterize IPA across physiological systems and disease contexts using widely available bulk RNA-seq data. These IPA events can be further integrated with other regulatory data sets to elucidate their interplay and functional significance.