Aberrant polyadenylation in HSPCs. (A) Heat-map showing genes (rows) with 3′ untranslated region (UTR) shortening (red) or lengthening (blue) in AA derived HSPCs. Bar plot shows the number of altered APA genes in each cell type (top). (B) Representative display of genes exhibiting 3′ UTR shortening (GATA2, JAK2) or lengthening (SETD2, TDP1) in HSC/MPP. Alternative usages of the 3′ UTR are highlighted by the dotted rectangles. (C) Top 5 enriched GO terms (ranked by logarithmic transformed adjusted P values, Benjamini-Hochberg correction) for altered APA genes in each HSPC subset. MBSP, multivesicular body sorting pathway; SDSA, synthesis-dependent strand annealing; TRBAN, transesterification reactions with bulged adenosine as nucleophile. (D) Enriched terms of differential APA genes between each patient with AA (n = 6) and 3 healthy controls (Ctrls) based on bulk RNA-seq. Dot sizes represent enrichment scores and dot colors represent logarithmic transformed adjusted P values (Benjamini-Hochberg correction). (E) Pie charts showing the number of ceRNA pairs shared between AA and Ctrl (blue). AA-specific ceRNA pairs (red) and control-specific ceRNA pairs (gray) have also been highlighted.