Single-cell genomics details the maturation block in BCP-ALL and identifies therapeutic vulnerabilities in DUX4-r cases.

B cell progenitor acute lymphoblastic leukemia (BCP-ALL) is the most common childhood malignancy, driven by multiple genetic alterations that cause maturation arrest and accumulation of abnormal progenitor B cells. Current treatment protocols with chemotherapy have led to favorable outcomes but are associated with significant toxicity and risk of side effects, highlighting the necessity for highly effective, less toxic, targeted drugs, even in subtypes with a favorable outcome. Here, we used multimodal single-cell sequencing to delineate the transcriptional, epigenetic, and immunophenotypic characteristics of 23 childhood BCP-ALLs, belonging to the BCR::ABL1-positive, ETV6::RUNX1-positive, high hyperdiploid, and recently discovered DUX4-rearranged (DUX4-r) subtypes. Projection of the ALL cells along the normal hematopoietic differentiation axis revealed a diversity in the maturation pattern between the different BCP-ALL subtypes. Whereas the BCR::ABL1-, ETV6::RUNX1-positive, and high hyperdiploidy cells mainly showed similarities to normal pro-B cells, the DUX4-r ALL cells also displayed transcriptional signatures resembling mature B cells. Focusing on the DUX4-r subtype, we found that the blast population displayed multilineage priming toward non-hematopoietic cells, myeloid, and T cell lineages, but also an activation of PI3K/AKT signaling that sensitized the cells to PI3K inhibition in vivo. Given the multilineage priming of the DUX4-r blasts with aberrant expression of the myeloid marker CD371 (CLL-1), we generated chimeric antigen receptor T cells, which effectively eliminated DUX4-r ALL cells in vivo. These results provide a detailed characterization of BCP-ALL at the single-cell level and reveal therapeutic vulnerabilities in the DUX4-r subtype with implications for the understanding of ALL biology and new therapeutic strategies.

Distinct Infiltration of T Cell Populations in Bladder Cancer Molecular Subtypes.

Bladder cancer is a heterogenous disease, and molecular subtyping is a promising method to capture this variability. Currently, the immune compartment in relation to subtypes is poorly characterized. Here, we analyzed the immune compartment in bladder tumors and normal bladder urothelium with a focus on T cell subpopulations using flow cytometry and RNA sequencing. The results were investigated in relation to tumor invasiveness (NMIBC/MIBC) and molecular subtypes according to the Lund Taxonomy system. Whereas the NMIBC/MIBC differed in the overall immune infiltration only, the molecular subtypes differed both in terms of immune infiltration and immune compartment compositions. The Basal/Squamous (Ba/Sq) and genomically unstable (GU) tumors displayed increased immune infiltration compared to urothelial-like (Uro) tumors. Additionally, the GU tumors had a higher proportion of regulatory T cells within the immune compartment compared to Uro tumors. Furthermore, sequencing showed higher levels of exhaustion in CD8+ T cells from GU tumors compared to both Uro tumors and the control. Although no such difference was detected at the transcriptomic level in Uro tumors compared to the controls, CD8+ T cells in Uro tumors showed higher expression of several exhaustion markers at the protein level. Taken together, our findings indicate that depending on the molecular subtype, different immunotherapeutic interventions might be warranted.

Constitutional and acquired genetic variants in ARID5B in pediatric B-cell precursor acute lymphoblastic leukemia.

Constitutional polymorphisms in ARID5B are associated with an increased risk of developing high hyperdiploid (HeH; 51-67 chromosomes) pediatric B-cell precursor acute lymphoblastic leukemia (BCP ALL). Here, we investigated constitutional and somatic ARID5B variants in 1335 BCP ALL cases from five different cohorts, with a particular focus on HeH cases. In 353 HeH ALL that were heterozygous for risk alleles and trisomic for chromosome 10, where ARID5B is located, a significantly higher proportion of risk allele duplication was seen for the SNPs rs7090445 (p = 0.009), rs7089424 (p = 0.005), rs7073837 (p = 0.03), and rs10740055 (p = 0.04). Somatic ARID5B deletions were seen in 16/1335 cases (1.2%), being more common in HeH than in other genetic subtypes (2.2% vs. 0.4%; p = 0.002). The expression of ARID5B in HeH cases with genomic deletions was reduced, consistent with a functional role in leukemogenesis. Whole-genome sequencing and RNA-sequencing in HeH revealed additional somatic events involving ARID5B, resulting in a total frequency of 3.6% of HeH cases displaying a somatic ARID5B aberration. Overall, our results show that both constitutional and somatic events in ARID5B are involved in the leukemogenesis of pediatric BCP ALL, particularly in the HeH subtype.

Sub-Nanomolar Detection of Oligonucleotides Using Molecular Beacons Immobilized on Lightguiding Nanowires.

The detection of oligonucleotides is a central step in many biomedical investigations. The most commonly used methods for detecting oligonucleotides often require concentration and amplification before detection. Therefore, developing detection methods with a direct read-out would be beneficial. Although commonly used for the detection of amplified oligonucleotides, fluorescent molecular beacons have been proposed for such direct detection. However, the reported limits of detection using molecular beacons are relatively high, ranging from 100 nM to a few µM, primarily limited by the beacon fluorescence background. In this study, we enhanced the relative signal contrast between hybridized and non-hybridized states of the beacons by immobilizing them on lightguiding nanowires. Upon hybridization to a complementary oligonucleotide, the fluorescence from the surface-bound beacon becomes coupled in the lightguiding nanowire core and is re-emitted at the nanowire tip in a narrower cone of light compared with the standard 4π emission. Prior knowledge of the nanowire positions allows for the continuous monitoring of fluorescence signals from each nanowire, which effectively facilitates the discrimination of signals arising from hybridization events against background signals. This resulted in improved signal-to-background and signal-to-noise ratios, which allowed for the direct detection of oligonucleotides at a concentration as low as 0.1 nM.

Multimodal classification of molecular subtypes in pediatric acute lymphoblastic leukemia.

Genomic analyses have redefined the molecular subgrouping of pediatric acute lymphoblastic leukemia (ALL). Molecular subgroups guide risk-stratification and targeted therapies, but outcomes of recently identified subtypes are often unclear, owing to limited cases with comprehensive profiling and cross-protocol studies. We developed a machine learning tool (ALLIUM) for the molecular subclassification of ALL in retrospective cohorts as well as for up-front diagnostics. ALLIUM uses DNA methylation and gene expression data from 1131 Nordic ALL patients to predict 17 ALL subtypes with high accuracy. ALLIUM was used to revise and verify the molecular subtype of 281 B-cell precursor ALL (BCP-ALL) cases with previously undefined molecular phenotype, resulting in a single revised subtype for 81.5% of these cases. Our study shows the power of combining DNA methylation and gene expression data for resolving ALL subtypes and provides a comprehensive population-based retrospective cohort study of molecular subtype frequencies in the Nordic countries.

Feasibility to use whole-genome sequencing as a sole diagnostic method to detect genomic aberrations in pediatric B-cell acute lymphoblastic leukemia.

Introduction: The suitability of whole-genome sequencing (WGS) as the sole method to detect clinically relevant genomic aberrations in B-cell acute lymphoblastic leukemia (ALL) was investigated with the aim of replacing current diagnostic methods. Methods: For this purpose, we assessed the analytical performance of 150 bp paired-end WGS (90x leukemia/30x germline). A set of 88 retrospective B-cell ALL samples were selected to represent established ALL subgroups as well as ALL lacking stratifying markers by standard-of-care (SoC), so-called B-other ALL. Results: Both the analysis of paired leukemia/germline (L/N)(n=64) as well as leukemia-only (L-only)(n=88) detected all types of aberrations mandatory in the current ALLTogether trial protocol, i.e., aneuploidies, structural variants, and focal copy-number aberrations. Moreover, comparison to SoC revealed 100% concordance and that all patients had been assigned to the correct genetic subgroup using both approaches. Notably, WGS could allocate 35 out of 39 B-other ALL samples to one of the emerging genetic subgroups considered in the most recent classifications of ALL. We further investigated the impact of high (90x; n=58) vs low (30x; n=30) coverage on the diagnostic yield and observed an equally perfect concordance with SoC; low coverage detected all relevant lesions. Discussion: The filtration of the WGS findings with a short list of genes recurrently rearranged in ALL was instrumental to extract the clinically relevant information efficiently. Nonetheless, the detection of DUX4 rearrangements required an additional customized analysis, due to multiple copies of this gene embedded in the highly repetitive D4Z4 region. We conclude that the diagnostic performance of WGS as the standalone method was remarkable and allowed detection of all clinically relevant genomic events in the diagnostic setting of B-cell ALL.

Combined GLUT1 and OXPHOS inhibition eliminates acute myeloid leukemia cells by restraining their metabolic plasticity.

Acute myeloid leukemia (AML) is initiated and propagated by leukemia stem cells (LSCs), a self-renewing population of leukemia cells responsible for therapy resistance. Hence, there is an urgent need to identify new therapeutic opportunities targeting LSCs. Here, we performed an in vivo CRISPR knockout screen to identify potential therapeutic targets by interrogating cell surface dependencies of LSCs. The facilitated glucose transporter type 1 (GLUT1) emerged as a critical in vivo metabolic dependency for LSCs in a murine MLL::AF9-driven model of AML. GLUT1 disruption by genetic ablation or pharmacological inhibition led to suppression of leukemia progression and improved survival of mice that received transplantation with LSCs. Metabolic profiling revealed that Glut1 inhibition suppressed glycolysis, decreased levels of tricarboxylic acid cycle intermediates and increased the levels of amino acids. This metabolic reprogramming was accompanied by an increase in autophagic activity and apoptosis. Moreover, Glut1 disruption caused transcriptional, morphological, and immunophenotypic changes, consistent with differentiation of AML cells. Notably, dual inhibition of GLUT1 and oxidative phosphorylation (OXPHOS) exhibited synergistic antileukemic effects in the majority of tested primary AML patient samples through restraining of their metabolic plasticity. In particular, RUNX1-mutated primary leukemia cells displayed striking sensitivity to the combination treatment compared with normal CD34+ bone marrow and cord blood cells. Collectively, our study reveals a GLUT1 dependency of murine LSCs in the bone marrow microenvironment and demonstrates that dual inhibition of GLUT1 and OXPHOS is a promising therapeutic approach for AML.

The complement receptor C3AR constitutes a novel therapeutic target in NPM1-mutated AML.

Mutated nucleophosmin 1 (NPM1) is the most common genetic alteration in acute myeloid leukemia (AML), found in ∼30% of cases. Although mutations in this gene are considered favorable according to current risk stratification guidelines, a large fraction of patients will experience relapse, demonstrating the urgent need for new treatment options. Therefore, we aimed to identify cell surface proteins specifically expressed on NPM1-mutated AML cells, allowing for potential targeting with antibody-based therapies. Herein, we report on an arrayed flow cytometry-based screen directed to 362 cell surface markers. In comparing the cell surface expression on NPM1-mutated AML cells with primitive (CD34+ CD38-) normal bone marrow cells, we identified the complement receptor C3AR as being specifically expressed in NPM1-mutated AML. By flow cytometry and single-cell RNA sequencing, we further show that normal hematopoietic stem and progenitor cells lack detectable C3AR gene and protein expression, making it particularly suitable as a target for antibody therapy. We also demonstrate that C3AR in combination with GPR56 distinguishes the leukemic stem cells (LSCs) in NPM1-mutated AML from the normal hematopoietic stem cells, defining the LSC population, as shown by transplantation into immunodeficient mice. Mechanistically, the stimulation of C3AR-expressing cells with C3a, the ligand of C3AR, leads to the activation of ERK1/2 and increased survival of AML cells, suggesting that this is an important signaling axis in this subtype of AML. Finally, we show that antibodies directed against C3AR efficiently elicit natural killer cell-mediated killing of primary AML cells ex vivo, highlighting C3AR as a candidate therapeutic target in NPM1-mutated AML.

High CD34 surface expression in BCP-ALL predicts poor induction therapy response and is associated with altered expression of genes related to cell migration and adhesion.

Minimal residual disease (MRD) constitutes the most important prognostic factor in B-cell precursor acute lymphoblastic leukemia (BCP-ALL). Flow cytometry is widely used in MRD assessment, yet little is known regarding the effect of different immunophenotypic subsets on outcome. In this study of 200 BCP-ALL patients, we found that a CD34-positive, CD38 dim-positive, nTdT dim-positive immunophenotype on the leukemic blasts was associated with poor induction therapy response and predicted an MRD level at the end of induction therapy (EOI) of ≥ 0.001. CD34 expression was strongly and positively associated with EOI MRD, whereas CD34-negative patients had a low relapse risk. Further, CD34 expression increased from diagnosis to relapse. CD34 is a stemness-associated cell-surface molecule, possibly involved in cell adhesion/migration or survival. Accordingly, genes associated with stemness were overrepresented among the most upregulated genes in CD34-positive leukemias, and protein-protein interaction networks showed an overrepresentation of genes associated with cell migration, cell adhesion, and negative regulation of apoptosis. The present work is the first to demonstrate a CD34-negative immunophenotype as a good prognostic factor in ALL, whereas high CD34 expression is associated with poor therapy response and an altered gene expression profile reminiscent of migrating cancer stem-like cells.

Transcriptional profiling demonstrates altered characteristics of CD8+ cytotoxic T-cells and regulatory T-cells in TP53-mutated acute myeloid leukemia.

BACKGROUND: Acute myeloid leukemia (AML) patients have limited effect from T-cell-based therapies, such as PD-1 and CTLA-4 blockade. However, recent data indicate that AML patients with TP53 mutation have higher immune infiltration and other immunomodulatory therapies could thus potentially be effective. Here, we performed the transcriptional analysis of distinct T-cell subpopulations from TP53-mutated AML to identify gene expression signatures suggestive of altered functional properties. METHODS: CD8+ cytotoxic T lymphocytes (CTLs), conventional helper T cells (Th), and regulatory T cells (Tregs) were sorted from peripheral blood of AML patients with TP53 mutation (n = 5) and healthy donors (n = 3), using FACS, and the different subpopulations were subsequently subjected to RNA-sequencing. Differentially expressed genes were identified and gene set enrichment analysis (GSEA) was performed to outline altered pathways and exhaustion status. Also, expression levels for a set of genes encoding established and emerging immuno-oncological targets were defined. RESULTS: The results showed altered transcriptional profiles for each of the T-cell subpopulations from TP53-mutated AML as compared to control subjects. IFN-α and IFN-γ signaling were stronger in TP53-mutated AML for both CTLs and Tregs. Furthermore, in TP53-mutated AML as compared to healthy controls, Tregs showed gene expression signatures suggestive of metabolic adaptation to their environment, whereas CTLs exhibited features of exhaustion/dysfunction with a stronger expression of TIM3 as well as enrichment of a gene set related to exhaustion. CONCLUSIONS: The results provide insights on mechanisms underlying the inadequate immune response to leukemic cells in TP53-mutated AML and open up for further exploration toward novel treatment regimens for these patients.

N-terminus DUX4-immunohistochemistry is a reliable methodology for the diagnosis of DUX4-fused B-lymphoblastic leukemia/lymphoma (N-terminus DUX4 IHC for DUX4-fused B-ALL).

B-lymphoblastic leukemia/lymphoma (B-ALL) is the most common pediatric malignancy and the most commonly diagnosed adult lymphoblastic leukemia. Recent advances have broadened the spectrum of B-ALL, with DUX4 gene fusions implicated in a subclass occurring in adolescents and young adults and harboring a favorable prognosis. DUX4 fusions have been challenging to identify. We aimed to determine whether expression of the DUX4 oncoprotein, as detected by targeted immunohistochemistry, might serve as a surrogate for molecular detection of DUX4 fusions in B-ALL. A cohort of investigational B-ALLs was generated with enrichment for DUX4 fusions by the inclusion of cases with characteristic demographic features and immunophenotypic properties. B-ALLs with mutually exclusive cytogenetics were collected. Immunohistochemical staining by a monoclonal antibody raised against the N-terminus of the DUX4 protein was performed. N-DUX4 immunohistochemistry demonstrated strong, crisp nuclear staining in blasts of seven investigational cases, six of which had nucleic acid material available for molecular evaluation. Five of these cases demonstrated RNA-seq DUX4-fusion positivity. One N-DUX4 immunohistochemistry positive case lacked a definitive DUX4-fusion by RNA-seq, though demonstrated a gene expression profile characteristic of DUX4-rearranged B-ALLs, a CD2+ immunophenotype, and a lack of staining by C-terminus DUX4 antibody immunohistochemistry. At least 83.3% [5/6] positive predictive value. N-DUX4 immunohistochemistry was negative in blasts of three RNA-seq DUX4-fusion-negative cases (3/3; 100% negative predictive value). B-ALLs with mutually exclusive cytogenetic profiles were all N-DUX4 negative (0/10, specificity 100%). N-DUX4 immunohistochemistry is reliable for the distinction of DUX4-rearranged B-ALLs from other B-ALLs. We recommend its use for subclassification of B-ALLs in adolescents and young adults and in B-ALLs that remain "not otherwise specified."

Transcriptomics paving the way for improved diagnostics and precision medicine of acute leukemia.

Transcriptional profiling of acute leukemia, specifically by RNA-sequencing or whole transcriptome sequencing (WTS), has provided fundamental insights into its underlying disease biology and allows unbiased detection of oncogenic gene fusions, as well as of gene expression signatures that can be used for improved disease classification. While used as a research tool for many years, RNA-sequencing is becoming increasingly used in clinical diagnostics. Here, we highlight key transcriptomic studies of acute lymphoblastic leukemia (ALL) and acute myeloid leukemia (AML) that have improved our biological understanding of these heterogeneous malignant disorders and have paved the way for translation into clinical diagnostics. Recent single-cell transcriptomic studies of ALL and AML, which provide new insights into the cellular ecosystem of acute leukemia and point to future clinical utility, are also reviewed. Finally, we discuss current challenges that need to be overcome for a more wide-spread adoption of RNA-sequencing in clinical diagnostics and how this technology significantly can aid the identification of genetic alterations in current guidelines and of newly emerging disease entities, some of which are critical to identify because of the availability of targeted therapies, thereby paving the way for improved precision medicine of acute leukemia.

Clinical and genomic characterization of patients diagnosed with the provisional entity acute myeloid leukemia with BCR-ABL1, a Swedish population-based study.

Acute myeloid leukemia (AML) with t(9;22)(q34;q11), also known as AML with BCR-ABL1, is a rare, provisional entity in the WHO 2016 classification and is considered a high-risk disease according to the European LeukemiaNet 2017 risk stratification. We here present a retrospective, population-based study of this disease entity from the Swedish Acute Leukemia Registry. By strict clinical inclusion criteria we aimed to identify genetic markers further distinguishing AML with t(9;22) as a separate entity. Twenty-five patients were identified and next-generation sequencing using a 54-gene panel was performed in 21 cases. Interestingly, no mutations were found in NPM1, FLT3, or DNMT3A, three frequently mutated genes in AML. Instead, RUNX1 was the most commonly mutated gene, with aberrations present in 38% of the cases compared to around 10% in de novo AML. Additional mutations were identified in genes involved in RNA splicing (SRSF2, SF3B1) and chromatin regulation (ASXL1, STAG2, BCOR, BCORL1). Less frequently, mutations were found in IDH2, NRAS, TET2, and TP53. The mutational landscape exhibited a similar pattern as recently described in patients with chronic myeloid leukemia (CML) in myeloid blast crisis (BC). Despite the concomitant presence of BCR-ABL1 and RUNX1 mutations in our cohort, both features of high-risk AML, the RUNX1-mutated cases showed a superior overall survival compared to RUNX1 wildtype cases. Our results suggest that the molecular characteristics of AML with t(9;22)/BCR-ABL1 and CML in myeloid BC are similar and do not support a distinction of the two disease entities based on their underlying molecular alterations.

Clonal competition within complex evolutionary hierarchies shapes AML over time.

Clonal heterogeneity and evolution has major implications for disease progression and relapse in acute myeloid leukemia (AML). To model clonal dynamics in vivo, we serially transplanted 23 AML cases to immunodeficient mice and followed clonal composition for up to 15 months by whole-exome sequencing of 84 xenografts across two generations. We demonstrate vast changes in clonality that both progress and reverse over time, and define five patterns of clonal dynamics: Monoclonal, Stable, Loss, Expansion and Burst. We also show that subclonal expansion in vivo correlates with a more adverse prognosis. Furthermore, clonal expansion enabled detection of very rare clones with AML driver mutations that were undetectable by sequencing at diagnosis, demonstrating that the vast majority of AML cases harbor multiple clones already at diagnosis. Finally, the rise and fall of related clones enabled deconstruction of the complex evolutionary hierarchies of the clones that compete to shape AML over time.

A high-content cytokine screen identifies myostatin propeptide as a positive regulator of primitive chronic myeloid leukemia cells.

Aberrantly expressed cytokines in the bone marrow (BM) niche are increasingly recognized as critical mediators of survival and expansion of leukemic stem cells. To identify regulators of primitive chronic myeloid leukemia (CML) cells, we performed a high-content cytokine screen using primary CD34+ CD38low chronic phase CML cells. Out of the 313 unique human cytokines evaluated, 11 were found to expand cell numbers ≥2-fold in a 7-day culture. Focusing on novel positive regulators of primitive CML cells, the myostatin antagonist myostatin propeptide gave the largest increase in cell expansion and was chosen for further studies. Herein, we demonstrate that myostatin propeptide expands primitive CML and normal BM cells, as shown by increased colony-forming capacity. For primary CML samples, retention of CD34-expression was also seen after culture. Furthermore, we show expression of MSTN by CML mesenchymal stromal cells, and that myostatin propeptide has a direct and instant effect on CML cells, independent of myostatin, by demonstrating binding of myostatin propeptide to the cell surface and increased phosphorylation of STAT5 and SMAD2/3. In summary, we identify myostatin propeptide as a novel positive regulator of primitive CML cells and corresponding normal hematopoietic cells.

PAX5 is part of a functional transcription factor network targeted in lymphoid leukemia.

One of the most frequently mutated proteins in human B-lineage leukemia is the transcription factor PAX5. These mutations often result in partial rather than complete loss of function of the transcription factor. While the functional dose of PAX5 has a clear connection to human malignancy, there is limited evidence for that heterozygote loss of PAX5 have a dramatic effect on the development and function of B-cell progenitors. One possible explanation comes from the finding that PAX5 mutated B-ALL often display complex karyotypes and additional mutations. Thus, PAX5 might be one component of a larger transcription factor network targeted in B-ALL. To investigate the functional network associated with PAX5 we used BioID technology to isolate proteins associated with this transcription factor in the living cell. This identified 239 proteins out of which several could be found mutated in human B-ALL. Most prominently we identified the commonly mutated IKZF1 and RUNX1, involved in the formation of ETV6-AML1 fusion protein, among the interaction partners. ChIP- as well as PLAC-seq analysis supported the idea that these factors share a multitude of target genes in human B-ALL cells. Gene expression analysis of mouse models and primary human leukemia suggested that reduced function of PAX5 increased the ability of an oncogenic form of IKZF1 or ETV6-AML to modulate gene expression. Our data reveals that PAX5 belong to a regulatory network frequently targeted by multiple mutations in B-ALL shedding light on the molecular interplay in leukemia cells.

Proteogenomics and Hi-C reveal transcriptional dysregulation in high hyperdiploid childhood acute lymphoblastic leukemia.

Hyperdiploidy, i.e. gain of whole chromosomes, is one of the most common genetic features of childhood acute lymphoblastic leukemia (ALL), but its pathogenetic impact is poorly understood. Here, we report a proteogenomic analysis on matched datasets from genomic profiling, RNA-sequencing, and mass spectrometry-based analysis of >8,000 genes and proteins as well as Hi-C of primary patient samples from hyperdiploid and ETV6/RUNX1-positive pediatric ALL. We show that CTCF and cohesin, which are master regulators of chromatin architecture, display low expression in hyperdiploid ALL. In line with this, a general genome-wide dysregulation of gene expression in relation to topologically associating domain (TAD) borders were seen in the hyperdiploid group. Furthermore, Hi-C of a limited number of hyperdiploid childhood ALL cases revealed that 2/4 cases displayed a clear loss of TAD boundary strength and 3/4 showed reduced insulation at TAD borders, with putative leukemogenic effects.

Arrayed molecular barcoding identifies TNFSF13 as a positive regulator of acute myeloid leukemia-initiating cells.

Dysregulation of cytokines in the bone marrow (BM) microenvironment promotes acute myeloid leukemia (AML) cell growth. Due to the complexity and low throughput of in vivo stem-cell based assays, studying the role of cytokines in the BM niche in a screening setting is challenging. Here, we developed an ex vivo cytokine screen using 11 arrayed molecular barcodes, allowing for a competitive in vivo readout of leukemia-initiating capacity. With this approach, we assessed the effect of 114 murine cytokines on MLL-AF9 AML mouse cells and identified the tumor necrosis factor ligand superfamily member 13 (TNFSF13) as a positive regulator of leukemia-initiating cells. By using Tnfsf13-/- recipient mice, we confirmed that TNFSF13 supports leukemia initiation also under physiological conditions. TNFSF13 was secreted by normal myeloid cells but not by leukemia mouse cells, suggesting that mature myeloid BM cells support leukemia cells by secreting TNFSF13. TNFSF13 supported leukemia cell proliferation in an NF-κB-dependent manner by binding TNFRSF17 and suppressed apoptosis. Moreover, TNFSF13 supported the growth and survival of several human myeloid leukemia cell lines, demonstrating that our findings translate to human disease. Taken together, using arrayed molecular barcoding, we identified a previously unrecognized role of TNFSF13 as a positive regulator of AML-initiating cells. The arrayed barcoded screening methodology is not limited to cytokines and leukemia, but can be extended to other types of ex vivo screens, where a multiplexed in vivo read-out of stem cell functionality is needed.