Most eukaryotic genes express alternative polyadenylation (APA) isoforms. A growing number of RNA sequencing methods, especially those used for single-cell transcriptome analysis, generate reads close to the polyadenylation site (PAS), termed nearSite reads, hence inherently containing information about APA isoform abundance. Here, we present a probabilistic model-based method named MAAPER to utilize nearSite reads for APA analysis. MAAPER predicts PASs with high accuracy and sensitivity and examines different types of APA events with robust statistics. We show MAAPER's performance with both bulk and single-cell data and its applicability in unpaired or paired experimental designs.

The proper balance of hematopoietic stem cell (HSC) self-renewal and differentiation is critical for normal hematopoiesis and is disrupted in hematologic malignancy. Among regulators of HSC fate, transcription factors have a well-defined central role, and mutations promote malignant transformation. More recently, studies have illuminated the importance of posttranscriptional regulation by RNA-binding proteins (RBPs) in hematopoiesis and leukemia development. However, the RBPs involved and the breadth of regulation are only beginning to be elucidated. Furthermore, the intersection between posttranscriptional regulation and hematopoietic transcription factor function is poorly understood. Here, we studied the posttranscriptional regulation of RUNX1, a key hematopoietic transcription factor. Alternative polyadenylation (APA) of RUNX1 produces functionally antagonistic protein isoforms (RUNX1a vs RUNX1b/c) that mediate HSC self-renewal vs differentiation, an RNA-processing event that is dysregulated in malignancy. Consequently, RBPs that regulate this event directly contribute to healthy and aberrant hematopoiesis. We modeled RUNX1 APA using a split GFP minigene reporter and confirmed the sensitivity of our model to detect changes in RNA processing. We used this reporter in a clustered regularly interspaced short palindromic repeats (CRISPR) screen consisting of single guide RNAs exclusively targeting RBPs and uncovered HNRNPA1 and KHDRBS1 as antagonistic regulators of RUNX1a isoform generation. Overall, our study provides mechanistic insight into the posttranscriptional regulation of a key hematopoietic transcription factor and identifies RBPs that may have widespread and important functions in hematopoiesis.

Long-lasting forms of synaptic plasticity that underlie learning and memory require new transcription and translation for their persistence. The remarkable polarity and compartmentalization of neurons raises questions about the spatial and temporal regulation of gene expression within neurons. Alternative cleavage and polyadenylation (APA) generates mRNA isoforms with different 3' untranslated regions (3'UTRs) and/or coding sequences. Changes in the 3'UTR composition of mRNAs can alter gene expression by regulating transcript localization, stability and/or translation, while changes in the coding sequences lead to mRNAs encoding distinct proteins. Using specialized 3' end deep sequencing methods, we undertook a comprehensive analysis of APA following induction of long-term potentiation (LTP) of mouse hippocampal CA3-CA1 synapses. We identified extensive LTP-induced APA changes, including a general trend of 3'UTR shortening and activation of intronic APA isoforms. Comparison with transcriptome profiling indicated that most APA regulatory events were uncoupled from changes in transcript abundance. We further show that specific APA regulatory events can impact expression of two molecules with known functions during LTP, including 3'UTR APA of Notch1 and intronic APA of Creb1. Together, our results reveal that activity-dependent APA provides an important layer of gene regulation during learning and memory.