Alternative splicing analysis of pediatric B-cell acute lymphoblastic leukemia compared to normal bone marrow derived pro-B cells

Aberrant splicing is a hallmark of leukemias with mutations in splicing factor (SF)-encoding genes. Here we investigated its prevalence in pediatric B-cell acute lymphoblastic leukemias (B-ALL), where SFs are not mutated. By comparing them to normal pro-B cells, we found thousands of aberrant local splice variations (LSVs) per sample, with 279 LSVs in 241 genes present in every comparison. These genes were enriched in RNA processing pathways and encoded ~100 SFs, e.g. hnRNPA1. hnRNPA1 3’UTR was most pervasively misspliced, yielding the transcript subject to nonsense-mediated decay. Thus, we knocked it down in B-lymphoblastoid cells, identified 213 hnRNPA1-dependent splicing events, and defined the hnRNPA1 splicing signature in pediatric leukemias. One of its elements was DICER1, a known tumor suppressor gene; its LSVs involved the 5’ UTR, suggestive of splicing as a mechanism of translational deregulation. Additionally, we searched for LSVs in other leukemia and lymphoma drivers and discovered 81 LSVs in 41 genes. 77 LSVs were confirmed using two large independent B-ALL RNA-seq datasets. In fact, the twenty most common B-ALL drivers showed higher prevalence of aberrant splicing than of somatic mutations. Thus, post-transcriptional deregulation of SF can drive widespread changes in B-ALL splicing and likely contribute to disease pathogenesis.