ALS/FTD-associated protein FUS induces mitochondrial dysfunction by preferentially sequestering respiratory chain complex mRNAs.

Dysregulation of the DNA/RNA-binding protein FUS causes certain subtypes of ALS/FTD by largely unknown mechanisms. Recent evidence has shown that FUS toxic gain of function due either to mutations or to increased expression can disrupt critical cellular processes, including mitochondrial functions. Here, we demonstrate that in human cells overexpressing wild-type FUS or expressing mutant derivatives, the protein associates with multiple mRNAs, and these are enriched in mRNAs encoding mitochondrial respiratory chain components. Notably, this sequestration leads to reduced levels of the encoded proteins, which is sufficient to bring about disorganized mitochondrial networks, reduced aerobic respiration and increased reactive oxygen species. We further show that mutant FUS associates with mitochondria and with mRNAs encoded by the mitochondrial genome. Importantly, similar results were also observed in fibroblasts derived from ALS patients with FUS mutations. Finally, we demonstrate that FUS loss of function does not underlie the observed mitochondrial dysfunction, and also provides a mechanism for the preferential sequestration of the respiratory chain complex mRNAs by FUS that does not involve sequence-specific binding. Together, our data reveal that respiratory chain complex mRNA sequestration underlies the mitochondrial defects characteristic of ALS/FTD and contributes to the FUS toxic gain of function linked to this disease spectrum.

CodonBert large language model for mRNA vaccines.

mRNA-based vaccines and therapeutics are gaining popularity and usage across a wide range of conditions. One of the critical issues when designing such mRNAs is sequence optimization. Even small proteins or peptides can be encoded by an enormously large number of mRNAs. The actual mRNA sequence can have a large impact on several properties including expression, stability, immunogenicity, and more. To enable the selection of an optimal sequence, we developed CodonBERT, a large language model (LLM) for mRNAs. Unlike prior models, CodonBERT uses codons as inputs which enables it to learn better representations. CodonBERT was trained using more than 10 million mRNA sequences from a diverse set of organisms. The resulting model captures important biological concepts. CodonBERT can also be extended to perform prediction tasks for various mRNA properties. CodonBERT outperforms previous mRNA prediction methods including on a new flu vaccine dataset.

The mRNA Export Receptor NXF1 Coordinates Transcriptional Dynamics, Alternative Polyadenylation, and mRNA Export.

Alternative polyadenylation (APA) produces mRNA isoforms with different 3' UTR lengths. Previous studies indicated that 3' end processing and mRNA export are intertwined in gene regulation. Here, we show that mRNA export factors generally facilitate usage of distal cleavage and polyadenylation sites (PASs), leading to long 3' UTR isoform expression. By focusing on the export receptor NXF1, which exhibits the most potent effect on APA in this study, we reveal several gene features that impact NXF1-dependent APA, including 3' UTR size, gene size, and AT content. Surprisingly, NXF1 downregulation results in RNA polymerase II (Pol II) accumulation at the 3' end of genes, correlating with its role in APA regulation. Moreover, NXF1 cooperates with CFI-68 to facilitate nuclear export of long 3' UTR isoform with UGUA motifs. Together, our work reveals important roles of NXF1 in coordinating transcriptional dynamics, 3' end processing, and nuclear export of long 3' UTR transcripts, implicating NXF1 as a nexus of gene regulation.

Alternative polyadenylation dysregulation contributes to the differentiation block of acute myeloid leukemia.

Posttranscriptional regulation has emerged as a driver for leukemia development and an avenue for therapeutic targeting. Among posttranscriptional processes, alternative polyadenylation (APA) is globally dysregulated across cancer types. However, limited studies have focused on the prevalence and role of APA in myeloid leukemia. Furthermore, it is poorly understood how altered poly(A) site usage of individual genes contributes to malignancy or whether targeting global APA patterns might alter oncogenic potential. In this study, we examined global APA dysregulation in patients with acute myeloid leukemia (AML) by performing 3' region extraction and deep sequencing (3'READS) on a subset of AML patient samples along with healthy hematopoietic stem and progenitor cells (HSPCs) and by analyzing publicly available data from a broad AML patient cohort. We show that patient cells exhibit global 3' untranslated region (UTR) shortening and coding sequence lengthening due to differences in poly(A) site (PAS) usage. Among APA regulators, expression of FIP1L1, one of the core cleavage and polyadenylation factors, correlated with the degree of APA dysregulation in our 3'READS data set. Targeting global APA by FIP1L1 knockdown reversed the global trends seen in patients. Importantly, FIP1L1 knockdown induced differentiation of t(8;21) cells by promoting 3'UTR lengthening and downregulation of the fusion oncoprotein AML1-ETO. In non-t(8;21) cells, FIP1L1 knockdown also promoted differentiation by attenuating mechanistic target of rapamycin complex 1 (mTORC1) signaling and reducing MYC protein levels. Our study provides mechanistic insights into the role of APA in AML pathogenesis and indicates that targeting global APA patterns can overcome the differentiation block in patients with AML.

The androgen receptor couples promoter recruitment of RNA processing factors to regulation of alternative polyadenylation at the 3' end of transcripts.

Prostate cancer (PC) relies on androgen receptor (AR) signaling. While hormonal therapy (HT) is efficacious, most patients evolve to an incurable castration-resistant stage (CRPC). To date, most proposed mechanisms of acquired resistance to HT have focused on AR transcriptional activity. Herein, we uncover a new role for the AR in alternative cleavage and polyadenylation (APA). Inhibition of the AR by Enzalutamide globally regulates APA in PC cells, with specific enrichment in genes related to transcription and DNA topology, suggesting their involvement in transcriptome reprogramming. AR inhibition selects promoter-distal polyadenylation sites (pAs) enriched in cis-elements recognized by the cleavage and polyadenylation specificity factor (CPSF) complex. Conversely, promoter-proximal intronic pAs relying on the cleavage stimulation factor (CSTF) complex are repressed. Mechanistically, Enzalutamide induces rearrangement of APA subcomplexes and impairs the interaction between CPSF and CSTF. AR inhibition also induces co-transcriptional CPSF recruitment to gene promoters, predisposing the selection of pAs depending on this complex. Importantly, the scaffold CPSF160 protein is up-regulated in CRPC cells and its depletion represses HT-induced APA patterns. These findings uncover an unexpected role for the AR in APA regulation and suggest that APA-mediated transcriptome reprogramming represents an adaptive response of PC cells to HT.

ALYREF links 3'-end processing to nuclear export of non-polyadenylated mRNAs.

The RNA-binding protein ALYREF plays key roles in nuclear export and also 3'-end processing of polyadenylated mRNAs, but whether such regulation also extends to non-polyadenylated RNAs is unknown. Replication-dependent (RD)-histone mRNAs are not polyadenylated, but instead end in a stem-loop (SL) structure. Here, we demonstrate that ALYREF prevalently binds a region next to the SL on RD-histone mRNAs. SL-binding protein (SLBP) directly interacts with ALYREF and promotes its recruitment. ALYREF promotes histone pre-mRNA 3'-end processing by facilitating U7-snRNP recruitment through physical interaction with the U7-snRNP-specific component Lsm11. Furthermore, ALYREF, together with other components of the TREX complex, enhances histone mRNA export. Moreover, we show that 3'-end processing promotes ALYREF recruitment and histone mRNA export. Together, our results point to an important role of ALYREF in coordinating 3'-end processing and nuclear export of non-polyadenylated mRNAs.

3'READS + RIP defines differential Staufen1 binding to alternative 3'UTR isoforms and reveals structures and sequence motifs influencing binding and polysome association.

Staufen1 (STAU1) is an RNA-binding protein (RBP) that interacts with double-stranded RNA structures and has been implicated in regulating different aspects of mRNA metabolism. Previous studies have indicated that STAU1 interacts extensively with RNA structures in coding regions (CDSs) and 3'-untranslated regions (3'UTRs). In particular, duplex structures formed within 3'UTRs by inverted-repeat Alu elements (IRAlus) interact with STAU1 through its double-stranded RNA-binding domains (dsRBDs). Using 3' region extraction and deep sequencing coupled to ribonucleoprotein immunoprecipitation (3'READS + RIP), together with reanalyzing previous STAU1 binding and RNA structure data, we delineate STAU1 interactions transcriptome-wide, including binding differences between alternative polyadenylation (APA) isoforms. Consistent with previous reports, RNA structures are dominant features for STAU1 binding to CDSs and 3'UTRs. Overall, relative to short 3'UTR counterparts, longer 3'UTR isoforms of genes have stronger STAU1 binding, most likely due to a higher frequency of RNA structures, including specific IRAlus sequences. Nevertheless, a sizable fraction of genes express transcripts showing the opposite trend, attributable to AU-rich sequences in their alternative 3'UTRs that may recruit antagonistic RBPs and/or destabilize RNA structures. Using STAU1-knockout cells, we show that strong STAU1 binding to mRNA 3'UTRs generally enhances polysome association. However, IRAlus generally have little impact on STAU1-mediated polysome association despite having strong interactions with the protein. Taken together, our work reveals complex interactions of STAU1 with its cognate RNA substrates. Our data also shed light on distinct post-transcriptional fates for the widespread APA isoforms in mammalian cells.

A compendium of conserved cleavage and polyadenylation events in mammalian genes.

Cleavage and polyadenylation is essential for 3' end processing of almost all eukaryotic mRNAs. Recent studies have shown widespread alternative cleavage and polyadenylation (APA) events leading to mRNA isoforms with different 3' UTRs and/or coding sequences. Here, we present a compendium of conserved cleavage and polyadenylation sites (PASs) in mammalian genes, based on approximately 1.2 billion 3' end sequencing reads from more than 360 human, mouse, and rat samples. We show that ∼80% of mammalian mRNA genes contain at least one conserved PAS, and ∼50% have conserved APA events. PAS conservation generally reduces promiscuous 3' end processing, stabilizing gene expression levels across species. Conservation of APA correlates with gene age, gene expression features, and gene functions. Genes with certain functions, such as cell morphology, cell proliferation, and mRNA metabolism, are particularly enriched with conserved APA events. Whereas tissue-specific genes typically have a low APA rate, brain-specific genes tend to evolve APA. In addition, we show enrichment of mRNA destabilizing motifs in alternative 3' UTR sequences, leading to substantial differences in mRNA stability between 3' UTR isoforms. Using conserved PASs, we reveal sequence motifs surrounding APA sites and a preference of adenosine at the cleavage site. Furthermore, we show that mutations of U-rich motifs around the PAS often accompany APA profile differences between species. Analysis of lncRNA PASs indicates a mechanism of PAS fixation through evolution of A-rich motifs. Taken together, our results present a comprehensive view of PAS evolution in mammals, and a phylogenic perspective on APA functions.

PolyA_DB 3 catalogs cleavage and polyadenylation sites identified by deep sequencing in multiple genomes.

PolyA_DB is a database cataloging cleavage and polyadenylation sites (PASs) in several genomes. Previous versions were based mainly on expressed sequence tags (ESTs), which had a limited amount and could lead to inaccurate PAS identification due to the presence of internal A-rich sequences in transcripts. Here, we present an updated version of the database based solely on deep sequencing data. First, PASs are mapped by the 3' region extraction and deep sequencing (3'READS) method, ensuring unequivocal PAS identification. Second, a large volume of data based on diverse biological samples increases PAS coverage by 3.5-fold over the EST-based version and provides PAS usage information. Third, strand-specific RNA-seq data are used to extend annotated 3' ends of genes to obtain more thorough annotations of alternative polyadenylation (APA) sites. Fourth, conservation information of PAS across mammals sheds light on significance of APA sites. The database (URL: http://www.polya-db.org/v3) currently holds PASs in human, mouse, rat and chicken, and has links to the UCSC genome browser for further visualization and for integration with other genomic data.

Novel mouse models of oculopharyngeal muscular dystrophy (OPMD) reveal early onset mitochondrial defects and suggest loss of PABPN1 may contribute to pathology.

Oculopharyngeal muscular dystrophy (OPMD) is a late onset disease caused by polyalanine expansion in the poly(A) binding protein nuclear 1 (PABPN1). Several mouse models have been generated to study OPMD; however, most of these models have employed transgenic overexpression of alanine-expanded PABPN1. These models do not recapitulate the OPMD patient genotype and PABPN1 overexpression could confound molecular phenotypes. We have developed a knock-in mouse model of OPMD (Pabpn1+/A17) that contains one alanine-expanded Pabpn1 allele under the control of the native promoter and one wild-type Pabpn1 allele. This mouse is the closest available genocopy of OPMD patients. We show that Pabpn1+/A17 mice have a mild myopathic phenotype in adult and aged animals. We examined early molecular and biochemical phenotypes associated with expressing native levels of A17-PABPN1 and detected shorter poly(A) tails, modest changes in poly(A) signal (PAS) usage, and evidence of mitochondrial damage in these mice. Recent studies have suggested that a loss of PABPN1 function could contribute to muscle pathology in OPMD. To investigate a loss of function model of pathology, we generated a heterozygous Pabpn1 knock-out mouse model (Pabpn1+/Δ). Like the Pabpn1+/A17 mice, Pabpn1+/Δ mice have mild histologic defects, shorter poly(A) tails, and evidence of mitochondrial damage. However, the phenotypes detected in Pabpn1+/Δ mice only partially overlap with those detected in Pabpn1+/A17 mice. These results suggest that loss of PABPN1 function could contribute to but may not completely explain the pathology detected in Pabpn1+/A17 mice.

Transcription elongation rate has a tissue-specific impact on alternative cleavage and polyadenylation in Drosophila melanogaster.

Alternative polyadenylation (APA) is a mechanism that generates multiple mRNA isoforms with different 3'UTRs and/or coding sequences from a single gene. Here, using 3' region extraction and deep sequencing (3'READS), we have systematically mapped cleavage and polyadenylation sites (PASs) in Drosophila melanogaster, expanding the total repertoire of PASs previously identified for the species, especially those located in A-rich genomic sequences. Cis-element analysis revealed distinct sequence motifs around fly PASs when compared to mammalian ones, including the greater enrichment of upstream UAUA elements and the less prominent presence of downstream UGUG elements. We found that over 75% of mRNA genes in Drosophila melanogaster undergo APA. The head tissue tends to use distal PASs when compared to the body, leading to preferential expression of APA isoforms with long 3'UTRs as well as with distal terminal exons. The distance between the APA sites and intron location of PAS are important parameters for APA difference between body and head, suggesting distinct PAS selection contexts. APA analysis of the RpII215C4 mutant strain, which harbors a mutant RNA polymerase II (RNAPII) with a slower elongation rate, revealed that a 50% decrease in transcriptional elongation rate leads to a mild trend of more usage of proximal, weaker PASs, both in 3'UTRs and in introns, consistent with the "first come, first served" model of APA regulation. However, this trend was not observed in the head, suggesting a different regulatory context in neuronal cells. Together, our data expand the PAS collection for Drosophila melanogaster and reveal a tissue-specific effect of APA regulation by RNAPII elongation rate.

Sizing up the poly(A) tail: insights from deep sequencing.

Global investigation of poly(A) tails has been hindered by technical challenges. In a recent advance, two groups developed deep sequencing methods to globally interrogate poly(A) tail length and sequence with high precision, opening new avenues for investigation of poly(A) tail functions in mRNA metabolism. Initial applications of these methods reveal insights into the relationship between poly(A) tail length and translational efficiency, and identify widespread uridylation and guanylation at the 3' ends of transcripts.

3'READS+, a sensitive and accurate method for 3' end sequencing of polyadenylated RNA.

Sequencing of the 3' end of poly(A)(+) RNA identifies cleavage and polyadenylation sites (pAs) and measures transcript expression. We previously developed a method, 3' region extraction and deep sequencing (3'READS), to address mispriming issues that often plague 3' end sequencing. Here we report a new version, named 3'READS+, which has vastly improved accuracy and sensitivity. Using a special locked nucleic acid oligo to capture poly(A)(+) RNA and to remove the bulk of the poly(A) tail, 3'READS+ generates RNA fragments with an optimal number of terminal A's that balance data quality and detection of genuine pAs. With improved RNA ligation steps for efficiency, the method shows much higher sensitivity (over two orders of magnitude) compared to the previous version. Using 3'READS+, we have uncovered a sizable fraction of previously overlooked pAs located next to or within a stretch of adenylate residues in human genes and more accurately assessed the frequency of alternative cleavage and polyadenylation (APA) in HeLa cells (∼50%). 3'READS+ will be a useful tool to accurately study APA and to analyze gene expression by 3' end counting, especially when the amount of input total RNA is limited.

Systematic profiling of poly(A)+ transcripts modulated by core 3' end processing and splicing factors reveals regulatory rules of alternative cleavage and polyadenylation.

Alternative cleavage and polyadenylation (APA) results in mRNA isoforms containing different 3' untranslated regions (3'UTRs) and/or coding sequences. How core cleavage/polyadenylation (C/P) factors regulate APA is not well understood. Using siRNA knockdown coupled with deep sequencing, we found that several C/P factors can play significant roles in 3'UTR-APA. Whereas Pcf11 and Fip1 enhance usage of proximal poly(A) sites (pAs), CFI-25/68, PABPN1 and PABPC1 promote usage of distal pAs. Strong cis element biases were found for pAs regulated by CFI-25/68 or Fip1, and the distance between pAs plays an important role in APA regulation. In addition, intronic pAs are substantially regulated by splicing factors, with U1 mostly inhibiting C/P events in introns near the 5' end of gene and U2 suppressing those in introns with features for efficient splicing. Furthermore, PABPN1 inhibits expression of transcripts with pAs near the transcription start site (TSS), a property possibly related to its role in RNA degradation. Finally, we found that groups of APA events regulated by C/P factors are also modulated in cell differentiation and development with distinct trends. Together, our results support an APA code where an APA event in a given cellular context is regulated by a number of parameters, including relative location to the TSS, splicing context, distance between competing pAs, surrounding cis elements and concentrations of core C/P factors.

Alternative cleavage and polyadenylation in spermatogenesis connects chromatin regulation with post-transcriptional control.

BACKGROUND: Most mammalian genes display alternative cleavage and polyadenylation (APA). Previous studies have indicated preferential expression of APA isoforms with short 3' untranslated regions (3'UTRs) in testes. RESULTS: By deep sequencing of the 3' end region of poly(A) + transcripts, we report widespread shortening of 3'UTR through APA during the first wave of spermatogenesis in mouse, with 3'UTR size being the shortest in spermatids. Using genes without APA as a control, we show that shortening of 3'UTR eliminates destabilizing elements, such as U-rich elements and transposable elements, which appear highly potent during spermatogenesis. We additionally found widespread regulation of APA events in introns and exons that can affect the coding sequence of transcripts and global activation of antisense transcripts upstream of the transcription start site, suggesting modulation of splicing and initiation of transcription during spermatogenesis. Importantly, genes that display significant 3'UTR shortening tend to have functions critical for further sperm maturation, and testis-specific genes display greater 3'UTR shortening than ubiquitously expressed ones, indicating functional relevance of APA to spermatogenesis. Interestingly, genes with shortened 3'UTRs tend to have higher RNA polymerase II and H3K4me3 levels in spermatids as compared to spermatocytes, features previously known to be associated with open chromatin state. CONCLUSIONS: Our data suggest that open chromatin may create a favorable cis environment for 3' end processing, leading to global shortening of 3'UTR during spermatogenesis. mRNAs with shortened 3'UTRs are relatively stable thanks to evasion of powerful mRNA degradation mechanisms acting on 3'UTR elements. Stable mRNAs generated in spermatids may be important for protein production at later stages of sperm maturation, when transcription is globally halted.

Analysis of alternative cleavage and polyadenylation by 3' region extraction and deep sequencing.

Alternative cleavage and polyadenylation (APA) generates diverse mRNA isoforms. We developed 3' region extraction and deep sequencing (3'READS) to address mispriming issues that commonly plague poly(A) site (pA) identification, and we used the method to comprehensively map pAs in the mouse genome. Thorough annotation of gene 3' ends revealed over 5,000 previously overlooked pAs (∼8% of total) flanked by A-rich sequences, underscoring the necessity of using an accurate tool for pA mapping. About 79% of mRNA genes and 66% of long noncoding RNA genes undergo APA, but these two gene types have distinct usage patterns for pAs in introns and upstream exons. Quantitative analysis of APA isoforms by 3'READS indicated that promoter-distal pAs, regardless of intron or exon locations, become more abundant during embryonic development and cell differentiation and that upregulated isoforms have stronger pAs, suggesting global modulation of the 3' end-processing activity in development and differentiation.

RNA-binding proteins in regulation of alternative cleavage and polyadenylation.

Almost all eukaryotic pre-mRNAs are processed at the 3' end by the cleavage and polyadenylation (C/P) reaction, which preludes termination of transcription and gives rise to the poly(A) tail of mature mRNA. Genomic studies in recent years have indicated that most eukaryotic mRNA genes have multiple cleavage and polyadenylation sites (pAs), leading to alternative cleavage and polyadenylation (APA) products. APA isoforms generally differ in their 3' untranslated regions (3' UTRs), but can also have different coding sequences (CDSs). APA expands the repertoire of transcripts expressed from the genome, and is highly regulated under various physiological and pathological conditions. Growing lines of evidence have shown that RNA-binding proteins (RBPs) play important roles in regulation of APA. Some RBPs are part of the machinery for C/P; others influence pA choice through binding to adjacent regions. In this chapter, we review cis elements and trans factors involved in C/P, the significance of APA, and increasingly elucidated roles of RBPs in APA regulation. We also discuss analysis of APA using transcriptome-wide techniques as well as molecular biology approaches.

Unraveling regulation and new components of human P-bodies through a protein interaction framework and experimental validation.

The cellular factors involved in mRNA degradation and translation repression can aggregate into cytoplasmic domains known as GW bodies or mRNA processing bodies (P-bodies). However, current understanding of P-bodies, especially the regulatory aspect, remains relatively fragmentary. To provide a framework for studying the mechanisms and regulation of P-body formation, maintenance, and disassembly, we compiled a list of P-body proteins found in various species and further grouped both reported and predicted human P-body proteins according to their functions. By analyzing protein-protein interactions of human P-body components, we found that many P-body proteins form complex interaction networks with each other and with other cellular proteins that are not recognized as P-body components. The observation suggests that these other cellular proteins may play important roles in regulating P-body dynamics and functions. We further used siRNA-mediated gene knockdown and immunofluorescence microscopy to demonstrate the validity of our in silico analyses. Our combined approach identifies new P-body components and suggests that protein ubiquitination and protein phosphorylation involving 14-3-3 proteins may play critical roles for post-translational modifications of P-body components in regulating P-body dynamics. Our analyses provide not only a global view of human P-body components and their physical interactions but also a wealth of hypotheses to help guide future research on the regulation and function of human P-bodies.

Neuronal Cells Display Distinct Stability Controls of Alternative Polyadenylation mRNA Isoforms, Long Non-Coding RNAs, and Mitochondrial RNAs.

RNA stability plays an important role in gene expression. Here, using 3' end sequencing of newly made and pre-existing poly(A)+ RNAs, we compare transcript stability in multiple human cell lines, including HEK293T, HepG2, and SH-SY5Y. We show that while mRNA stability is generally conserved across the cell lines, specific transcripts having a high GC content and possibly more stable secondary RNA structures are relatively more stable in SH-SY5Y cells compared to the other 2 cell lines. These features also differentiate stability levels of alternative polyadenylation (APA) 3'UTR isoforms in a cell type-specific manner. Using differentiation of a neural stem cell line as a model, we show that mRNA stability difference could contribute to gene expression changes in neurogenesis and confirm the neuronal identity of SH-SY5Y cells at both gene expression and APA levels. In addition, compared to transcripts using 3'-most exon cleavage/polyadenylation sites (PASs), those using intronic PASs are generally less stable, especially when the PAS-containing intron is large and has a strong 5' splice site, suggesting that intronic polyadenylation mostly plays a negative role in gene expression. Interestingly, the differential mRNA stability among APA isoforms appears to buffer PAS choice in these cell lines. Moreover, we found that several other poly(A)+ RNA species, including promoter-associated long noncoding RNAs and transcripts encoded by the mitochondrial genome, are more stable in SH-SY5Y cells than the other 2 cell lines, further highlighting distinct RNA metabolism in neuronal cells. Together, our results indicate that distinct RNA stability control in neuronal cells may contribute to the gene expression and APA programs that define their cell identity.

MAAPER: model-based analysis of alternative polyadenylation using 3' end-linked reads.

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.