Current classifications (World Health Organization-HAEM5/ICC) define up to 26 molecular B-cell precursor acute lymphoblastic leukemia (BCP-ALL) disease subtypes by genomic driver aberrations and corresponding gene expression signatures. Identification of driver aberrations by transcriptome sequencing (RNA-Seq) is well established, while systematic approaches for gene expression analysis are less advanced. Therefore, we developed ALLCatchR, a machine learning-based classifier using RNA-Seq gene expression data to allocate BCP-ALL samples to all 21 gene expression-defined molecular subtypes. Trained on n = 1869 transcriptome profiles with established subtype definitions (4 cohorts; 55% pediatric / 45% adult), ALLCatchR allowed subtype allocation in 3 independent hold-out cohorts (n = 1018; 75% pediatric / 25% adult) with 95.7% accuracy (averaged sensitivity across subtypes: 91.1% / specificity: 99.8%). High-confidence predictions were achieved in 83.7% of samples with 98.9% accuracy. Only 1.2% of samples remained unclassified. ALLCatchR outperformed existing tools and identified novel driver candidates in previously unassigned samples. Additional modules provided predictions of samples blast counts, patient's sex, and immunophenotype, allowing the imputation in cases where these information are missing. We established a novel RNA-Seq reference of human B-lymphopoiesis using 7 FACS-sorted progenitor stages from healthy bone marrow donors. Implementation in ALLCatchR enabled projection of BCP-ALL samples to this trajectory. This identified shared proximity patterns of BCP-ALL subtypes to normal lymphopoiesis stages, extending immunophenotypic classifications with a novel framework for developmental comparisons of BCP-ALL. ALLCatchR enables RNA-Seq routine application for BCP-ALL diagnostics with systematic gene expression analysis for accurate subtype allocation and novel insights into underlying developmental trajectories.

Childhood Leukaemia Investigation Prague 2nd Faculty of Medicine Charles University and University Hospital Motol Prague Czech Republic

Clinical Research Unit CATCH ALL funded by the Deutsche Forschungsgemeinschaft Bonn Germany

Department of Human Genetics Hannover Medical School Hannover Germany

Department of Medicine 2 Hematology Oncology Goethe University Hospital Frankfurt M Germany

Department of Pediatrics University Hospital Schleswig Holstein Kiel Germany

Institute of Experimental Medicine Research Group Medical Systems Biology Christian Albrechts University Kiel Germany

Medical Department 2 Hematology and Oncology University Hospital Schleswig Holstein Kiel Germany

MLL Munich Leukemia Laboratory Munich Germany

See more in PubMed

Alaggio R, Amador C, Anagnostopoulos I, et al. . The 5th edition of the World Health Organization classification of haematolymphoid tumours: lymphoid neoplasms. Leukemia. 2022;36:1720–1748. PubMed PMC

Arber DA, Orazi A, Hasserjian RP, et al. . International consensus classification of myeloid neoplasms and acute leukemia: integrating morphological, clinical, and genomic data. Blood. 2022;140:1200–1228. PubMed PMC

Bastian L, Schroeder MP, Eckert C, et al. . PAX5 biallelic genomic alterations define a novel subgroup of B-cell precursor acute lymphoblastic leukemia. Leukemia. 2019;33:1895–1909. PubMed

Bastian L, Hartmann AM, Beder T, et al. . UBTF::ATXN7L3 gene fusion defines novel B cell precursor ALL subtype with CDX2 expression and need for intensified treatment. Leukemia. 2022;36:1676–1680. PubMed PMC

Gu Z, Churchman ML, Roberts KG, et al. . PAX5-driven subtypes of B-progenitor acute lymphoblastic leukemia. Nat Genet. 2019;51:296–307. PubMed PMC

Zaliova M, Stuchly J, Winkowska L, et al. . Genomic landscape of pediatric B-other acute lymphoblastic leukemia in a consecutive European cohort. Haematologica. 2019;104:1396–1406. PubMed PMC

Chouvarine P, Antić Z, Lentes J, et al. . Transcriptional and mutational profiling of B-other acute lymphoblastic leukemia for improved diagnostics. Cancers. 2021;13:5653. PubMed PMC

Walter W, Shahswar R, Stengel A, et al. . Clinical application of whole transcriptome sequencing for the classification of patients with acute lymphoblastic leukemia. BMC Cancer. 2021;21:886. PubMed PMC

Nicorici D, Satalan M, Edgren H, et al. . FusionCatcher – a tool for finding somatic fusion genes in paired-end RNA-sequencing data. bioRxiv. 2014:011650.

Rehn J, Mayoh C, Heatley SL, et al. . RaScALL: Rapid (Ra) screening (Sc) of RNA-seq data for prognostically significant genomic alterations in acute lymphoblastic leukaemia (ALL). PLoS Genet. 2022;18:e1010300. PubMed PMC

Bařinka J, Hu Z, Wang L, et al. . RNAseqCNV: analysis of large-scale copy number variations from RNA-seq data. Leukemia. 2022;36:1492–1498. PubMed PMC

Schmidt B, Brown LM, Ryland GL, et al. . ALLSorts: an RNA-Seq subtype classifier for B-cell acute lymphoblastic leukemia. Blood Adv. 2022;6:4093–4097. PubMed PMC

Mäkinen V-P, Rehn J, Breen J, et al. . Multi-cohort transcriptomic subtyping of B-cell acute lymphoblastic leukemia. Int J Mol Sci. 2022;23:4574. PubMed PMC

Krali O, Marincevic-Zuniga Y, Arvidsson G, et al. . Multimodal classification of molecular subtypes in pediatric acute lymphoblastic leukemia. Hematology. 2023:2023.03.24.23287613. PubMed PMC

Boer JM, Marchante JRM, Evans WE, et al. . BCR-ABL1-like cases in pediatric acute lymphoblastic leukemia: a comparison between DCOG/Erasmus MC and COG/St. Jude signatures. Haematologica. 2015;100:e354–e357. PubMed PMC

Tibshirani R. The lasso method for variable selection in the cox model. Stat Med. 1997;16:385–395. PubMed

Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33:1–22. PubMed PMC

Kursa MB, Rudnicki WR. Feature selection with the boruta package. J Stat Softw. 2010;36:1–13.

Kuhn M. Building predictive models in R using the caret package. J Stat Softw. 2008;28:1–26. PubMed

Foroutan M, Bhuva DD, Lyu R, et al. . Single sample scoring of molecular phenotypes. BMC Bioinf. 2018;19:404. PubMed PMC

Harrison CJ, Moorman AV, Schwab C, et al. . An international study of intrachromosomal amplification of chromosome 21 (iAMP21): cytogenetic characterization and outcome. Leukemia. 2014;28:1015–1021. PubMed PMC

Schwab CJ, Chilton L, Morrison H, et al. . Genes commonly deleted in childhood B-cell precursor acute lymphoblastic leukemia: association with cytogenetics and clinical features. Haematologica. 2013;98:1081–1088. PubMed PMC

Bene MC, Castoldi G, Knapp W, et al. . Proposals for the immunological classification of acute leukemias. European Group for the Immunological Characterization of Leukemias (EGIL). Leukemia. 1995;9:1783–1786. PubMed

van Zelm MC, van der Burg M, de Ridder D, et al. . Ig gene rearrangement steps are initiated in early human precursor B cell subsets and correlate with specific transcription factor expression. J Immunol Baltim Md 1950. 2005;175:5912–5922. PubMed

Zoutman WH, Nell RJ, Versluis M, et al. . A novel digital PCR-based method to quantify (switched) B cells reveals the extent of allelic involvement in different recombination processes in the IGH locus. Mol Immunol. 2022;145:109–123. PubMed

Wagener R, López C, Kleinheinz K, et al. . IG-MYC+ neoplasms with precursor B-cell phenotype are molecularly distinct from Burkitt lymphomas. Blood. 2018;132:2280–2285. PubMed PMC

Chen C, Yu W, Alikarami F, et al. . Single-cell multiomics reveals increased plasticity, resistant populations, and stem-cell-like blasts in KMT2A-rearranged leukemia. Blood. 2022;139:2198–2211. PubMed PMC

Jia Z, Hu Z, Damirchi B, et al. . Characterization of PAX5 Mutations in B Progenitor Acute Lymphoblastic Leukemia. Blood. 2022;140:1001–1002.

Leek JT, Johnson WE, Parker HS, et al. . The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28:882–883. PubMed PMC

Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. PubMed PMC

Sinclair PB, Ryan S, Bashton M, et al. . SH2B3 inactivation through CN-LOH 12q is uniquely associated with B-cell precursor ALL with iAMP21 or other chromosome 21 gain. Leukemia. 2019;33:1881–1894. PubMed PMC

McClure BJ, Heatley SL, Kok CH, et al. . Pre-B acute lymphoblastic leukaemia recurrent fusion, EP300-ZNF384, is associated with a distinct gene expression. Br J Cancer. 2018;118:1000–1004. PubMed PMC

Yasuda T, Sanada M, Kawazu M, et al. . Two novel high-risk adult B-cell acute lymphoblastic leukemia subtypes with high expression of CDX2 and IDH1/2 mutations. Blood. 2022;139:1850–1862. PubMed

Passet M, Kim R, Gachet S, et al. . Concurrent CDX2 cis-deregulation and UBTF-ATXN7L3 fusion define a novel high-risk subtype of B-cell ALL. Blood. 2022;139:3505–3518. PubMed PMC

Kimura S, Montefiori L, Iacobucci I, et al. . Enhancer retargeting of CDX2 and UBTF::ATXN7L3 define a subtype of high-risk B-progenitor acute lymphoblastic leukemia. Blood. 2022;139:3519–3531. PubMed PMC

Li J-F, Dai Y-T, Lilljebjörn H, et al. . Transcriptional landscape of B cell precursor acute lymphoblastic leukemia based on an international study of 1,223 cases. Proc Natl Acad Sci U S A. 2018;115:E11711–E11720. PubMed PMC

Akasaka T, Balasas T, Russell LJ, et al. . Five members of the CEBP transcription factor family are targeted by recurrent IGH translocations in B-cell precursor acute lymphoblastic leukemia (BCP-ALL). Blood. 2007;109:3451–3461. PubMed

Zaliova M, Potuckova E, Lukes J, et al. . Frequency and prognostic impact of ZEB2 H1038 and 1072 mutations in childhood B-other acute lymphoblastic leukemia. Haematologica. 2020;106:886–890. PubMed PMC

Khabirova E, Jardine L, Coorens THH, et al. . Single-cell transcriptomics reveals a distinct developmental state of KMT2A-rearranged infant B-cell acute lymphoblastic leukemia. Nat Med. 2022;28:743–751. PubMed PMC

Uhrig S, Ellermann J, Walther T, et al. . Accurate and efficient detection of gene fusions from RNA sequencing data. Genome Res. 2021;31:448–460. PubMed PMC