Single-cell multiomics reveals increased plasticity, resistant populations, and stem-cell-like blasts in KMT2A-rearranged leukemia

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

ABSTRACT

B-cell progenitor acute lymphoblastic leukemia (BCP-ALL) is the most common childhood malignancy and is 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+, ETV6::RUNX1+, 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. Although the BCR::ABL1+, ETV6::RUNX1+, and high hyperdiploidy cells mainly showed similarities to normal pro-B cells, 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 not only multilineage priming toward nonhematopoietic cells, myeloid, and T-cell lineages, but also an activation of phosphatidylinositol 3-kinase (PI3K)/AKT signaling that sensitized the cells to PI3K inhibition in vivo. Given the multilineage priming of DUX4-r blasts with aberrant expression of 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.

Immune escape of B-cell lymphoblastic leukemic cells through a lineage switch to acute myeloid leukemia

Acute leukemia (AL) with a lineage switch (LS) is associated with poor prognosis. The predisposing factors of LS are unknown, apart from KMT2A rearrangements that have been reported to be associated with LS. Herein, we present two cases and review all 104 published cases to identify risk factors for LS. Most of the patients (75.5%) experienced a switch from the lymphoid phenotype to the myeloid phenotype. Eighteen patients (17.0%) experienced a transformation from acute myelogenous leukemia (AML) to acute lymphoblastic leukemia (ALL). Forty-nine (46.2%) patients carried a KMT2A rearrangement. Most of the cases involved LS from B-cell ALL (B-ALL) to AML (59.4%), and 49 patients (46.2%) carried KMT2A-rearrangements. Forty patients (37.7%) received lineage-specific immunotherapy. Our findings suggest that the prevalence of KMT2A rearrangements together with the lineage-specific immunotherapy may trigger LS, which supports the thesis of the existence of leukemia stem cells that are capable of lymphoid or myeloid differentiation.

Treatment-induced stemness and lineage plasticity in driving prostate cancer therapy resistance

Most human cancers are heterogeneous consisting of cancer cells at different epigenetic and transcriptional states and with distinct phenotypes, functions, and drug sensitivities. This inherent cancer cell heterogeneity contributes to tumor resistance to clinical treatment, especially the molecularly targeted therapies such as tyrosine kinase inhibitors (TKIs) and androgen receptor signaling inhibitors (ARSIs). Therapeutic interventions, in turn, induce lineage plasticity (also called lineage infidelity) in cancer cells that also drives therapy resistance. In this Perspective, we focus our discussions on cancer cell lineage plasticity manifested as treatment-induced switching of epithelial cancer cells to basal/stem-like, mesenchymal, and neural lineages. We employ prostate cancer (PCa) as the prime example to highlight ARSI-induced lineage plasticity during and towards development of castration-resistant PCa (CRPC). We further discuss how the tumor microenvironment (TME) influences therapy-induced lineage plasticity. Finally, we offer an updated summary on the regulators and mechanisms driving cancer cell lineage infidelity, which should be therapeutically targeted to extend the therapeutic window and improve patients' survival.

KMT2A oncoproteins induce epigenetic resistance to targeted therapies

Chromosomal translocations involving the Lysine-Methyl-Transferase-2A ( KMT2A ) locus generate potent oncogenic fusion proteins (oncoproteins) that disrupt regulation of developmental gene expression. By profiling the oncoprotein-target sites of 36 broadly representative KMT2A -rearranged leukemia samples, including three samples that underwent a lymphoid-to-myeloid lineage-switching event in response to therapy, we find the genomic enrichment of the oncoprotein is highly variable between samples and subject to dynamic regulation. At high levels of expression, the oncoproteins preferentially activate either an acute lymphoblastic leukemia (ALL) program, enriched for pro-B-cell genes, or an acute myeloid leukemia (AML) program, enriched for hematopoietic-stem-cell genes. The fusion-partner-specific-binding patterns over these gene sets are highly correlated with the prevalence of each mutation in ALL versus AML. In lineage-switching samples the oncoprotein levels are reduced and the oncoproteins preferentially activate granulocyte-monocyte progenitor (GMP) genes. In a sample that lineage switched during treatment with the menin inhibitor revumenib, the oncoprotein and menin are reduced to undetectable levels, but ENL, a transcriptional cofactor of the oncoprotein, persists on numerous oncoprotein-target loci, including genes in the GMP-like lineage-switching program. We propose KMT2A oncoproteins promote lineage-switching events through dynamic chromatin binding and can induce epigenetic lesions, marked by ENL, that support resistance to targeted therapies.

scMaui: a widely applicable deep learning framework for single-cell multiomics integration in the presence of batch effects and missing data

The recent advances in high-throughput single-cell sequencing have created an urgent demand for computational models which can address the high complexity of single-cell multiomics data. Meticulous single-cell multiomics integration models are required to avoid biases towards a specific modality and overcome sparsity. Batch effects obfuscating biological signals must also be taken into account. Here, we introduce a new single-cell multiomics integration model, Single-cell Multiomics Autoencoder Integration (scMaui) based on variational product-of-experts autoencoders and adversarial learning. scMaui calculates a joint representation of multiple marginal distributions based on a product-of-experts approach which is especially effective for missing values in the modalities. Furthermore, it overcomes limitations seen in previous VAE-based integration methods with regard to batch effect correction and restricted applicable assays. It handles multiple batch effects independently accepting both discrete and continuous values, as well as provides varied reconstruction loss functions to cover all possible assays and preprocessing pipelines. We demonstrate that scMaui achieves superior performance in many tasks compared to other methods. Further downstream analyses also demonstrate its potential in identifying relations between assays and discovering hidden subpopulations.

Sharing Data from the Human Tumor Atlas Network through Standards, Infrastructure, and Community Engagement

The Data Coordinating Center (DCC) of the Human Tumor Atlas Network (HTAN) has played a crucial role in enabling the broad sharing and effective utilization of HTAN data within the scientific community. Data from the first phase of HTAN are now available publicly. We describe the diverse datasets and modalities shared, multiple access routes to HTAN assay data and metadata, data standards, technical infrastructure and governance approaches, as well as our approach to sustained community engagement. HTAN data can be accessed via the HTAN Portal, explored in visualization tools-including CellxGene, Minerva, and cBioPortal-and analyzed in the cloud through the NCI Cancer Research Data Commons nodes. We have developed a streamlined infrastructure to ingest and disseminate data by leveraging the Synapse platform. Taken together, the HTAN DCC's approach demonstrates a successful model for coordinating, standardizing, and disseminating complex cancer research data via multiple resources in the cancer data ecosystem, offering valuable insights for similar consortia, and researchers looking to leverage HTAN data.

Concepts in B cell acute lymphoblastic leukemia pathogenesis

B cell acute lymphoblastic leukemia (B-ALL) arises from genetic alterations impacting B cell progenitors, ultimately leading to clinically overt disease. Extensive collaborative efforts in basic and clinical research have significantly improved patient prognoses. Nevertheless, a subset of patients demonstrate resistance to conventional chemotherapeutic approaches and emerging immunotherapeutic interventions. This review highlights the mechanistic underpinnings governing B-ALL transformation. Beginning with exploring normative B cell lymphopoiesis, we delineate the influence of recurrent germline and somatic genetic aberrations on the perturbation of B cell progenitor differentiation and protumorigenic signaling, thereby facilitating the neoplastic transformation underlying B-ALL progression. Additionally, we highlight recent advances in the multifaceted landscape of B-ALL, encompassing metabolic reprogramming, microbiome influences, inflammation, and the discernible impact of socioeconomic and racial disparities on B-ALL transformation and patient survival.

Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging

Non-hematopoietic cells are essential contributors to hematopoiesis. However, heterogeneity and spatial organization of these cells in human bone marrow remain largely uncharacterized. We used single-cell RNA sequencing (scRNA-seq) to profile 29,325 non-hematopoietic cells and discovered nine transcriptionally distinct subtypes. We simultaneously profiled 53,417 hematopoietic cells and predicted their interactions with non-hematopoietic subsets. We employed co-detection by indexing (CODEX) to spatially profile over 1.2 million cells. We integrated scRNA-seq and CODEX data to link predicted cellular signaling with spatial proximity. Our analysis revealed a hyperoxygenated arterio-endosteal neighborhood for early myelopoiesis, and an adipocytic localization for early hematopoietic stem and progenitor cells (HSPCs). We used our CODEX atlas to annotate new images and uncovered mesenchymal stromal cell (MSC) expansion and spatial neighborhoods co-enriched for leukemic blasts and MSCs in acute myeloid leukemia (AML) patient samples. This spatially resolved, multiomic atlas of human bone marrow provides a reference for investigation of cellular interactions that drive hematopoiesis.

Sequential drug treatment targeting cell cycle and cell fate regulatory programs blocks non-genetic cancer evolution in acute lymphoblastic leukemia

BACKGROUND

Targeted therapies exploiting vulnerabilities of cancer cells hold promise for improving patient outcome and reducing side-effects of chemotherapy. However, efficacy of precision therapies is limited in part because of tumor cell heterogeneity. A better mechanistic understanding of how drug effect is linked to cancer cell state diversity is crucial for identifying effective combination therapies that can prevent disease recurrence.

RESULTS

Here, we characterize the effect of G2/M checkpoint inhibition in acute lymphoblastic leukemia (ALL) and demonstrate that WEE1 targeted therapy impinges on cell fate decision regulatory circuits. We find the highest inhibition of recovery of proliferation in ALL cells with KMT2A-rearrangements. Single-cell RNA-seq and ATAC-seq of RS4;11 cells harboring KMT2A::AFF1, treated with the WEE1 inhibitor AZD1775, reveal diversification of cell states, with a fraction of cells exhibiting strong activation of p53-driven processes linked to apoptosis and senescence, and disruption of a core KMT2A-RUNX1-MYC regulatory network. In this cell state diversification induced by WEE1 inhibition, a subpopulation transitions to a drug tolerant cell state characterized by activation of transcription factors regulating pre-B cell fate, lipid metabolism, and pre-BCR signaling in a reversible manner. Sequential treatment with BCR-signaling inhibitors dasatinib, ibrutinib, or perturbing metabolism by fatostatin or AZD2014 effectively counteracts drug tolerance by inducing cell death and repressing stemness markers.

CONCLUSIONS

Collectively, our findings provide new insights into the tight connectivity of gene regulatory programs associated with cell cycle and cell fate regulation, and a rationale for sequential administration of WEE1 inhibitors with low toxicity inhibitors of pre-BCR signaling or metabolism.

Single-Cell Omics for Transcriptome CHaracterization (SCOTCH): isoform-level characterization of gene expression through long-read single-cell RNA sequencing

The advent of long-read single-cell transcriptome sequencing (lr-scRNA-Seq) represents a significant leap forward in single-cell genomics. With the recent introduction of R10 flowcells by Oxford Nanopore, we propose that previous computational methods designed to handle high sequencing error rates are no longer relevant, and that the prevailing approach using short reads to compile "barcode space" (candidate barcode list) to de-multiplex long reads are no longer necessary. Instead, computational methods should now shift focus on harnessing the unique benefits of long reads to analyze transcriptome complexity. In this context, we introduce a comprehensive suite of computational methods named Single-Cell Omics for Transcriptome CHaracterization (SCOTCH). Our method is compatible with the single-cell library preparation platform from both 10X Genomics and Parse Biosciences, facilitating the analysis of special cell populations, such as neurons, hepatocytes and developing cardiomyocytes. We specifically re-formulated the transcript mapping problem with a compatibility matrix and addressed the multiple-mapping issue using probabilistic inference, which allows the discovery of novel isoforms as well as the detection of differential isoform usage between cell populations. We evaluated SCOTCH through analysis of real data across different combinations of single-cell libraries and sequencing technologies (10X + Illumina, Parse + Illumina, 10X + Nanopore_R9, 10X + Nanopore_R10, Parse + Nanopore_R10), and showed its ability to infer novel biological insights on cell type-specific isoform expression. These datasets enhance the availability of publicly available data for continued development of computational approaches. In summary, SCOTCH allows extraction of more biological insights from the new advancements in single-cell library construction and sequencing technologies, facilitating the examination of transcriptome complexity at the single-cell level.

Aberrant stem cell and developmental programs in pediatric leukemia

Hematopoiesis is a finely orchestrated process, whereby hematopoietic stem cells give rise to all mature blood cells. Crucially, they maintain the ability to self-renew and/or differentiate to replenish downstream progeny. This process starts at an embryonic stage and continues throughout the human lifespan. Blood cancers such as leukemia occur when normal hematopoiesis is disrupted, leading to uncontrolled proliferation and a block in differentiation of progenitors of a particular lineage (myeloid or lymphoid). Although normal stem cell programs are crucial for tissue homeostasis, these can be co-opted in many cancers, including leukemia. Myeloid or lymphoid leukemias often display stem cell-like properties that not only allow proliferation and survival of leukemic blasts but also enable them to escape treatments currently employed to treat patients. In addition, some leukemias, especially in children, have a fetal stem cell profile, which may reflect the developmental origins of the disease. Aberrant fetal stem cell programs necessary for leukemia maintenance are particularly attractive therapeutic targets. Understanding how hijacked stem cell programs lead to aberrant gene expression in place and time, and drive the biology of leukemia, will help us develop the best treatment strategies for patients.

A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma

Pediatric high-grade glioma (pHGG) is an incurable central nervous system malignancy that is a leading cause of pediatric cancer death. While pHGG shares many similarities to adult glioma, it is increasingly recognized as a molecularly distinct, yet highly heterogeneous disease. In this study, we longitudinally profiled a molecularly diverse cohort of 16 pHGG patients before and after standard therapy through single-nucleus RNA and ATAC sequencing, whole-genome sequencing, and CODEX spatial proteomics to capture the evolution of the tumor microenvironment during progression following treatment. We found that the canonical neoplastic cell phenotypes of adult glioblastoma are insufficient to capture the range of tumor cell states in a pediatric cohort and observed differential tumor-myeloid interactions between malignant cell states. We identified key transcriptional regulators of pHGG cell states and did not observe the marked proneural to mesenchymal shift characteristic of adult glioblastoma. We showed that essential neuromodulators and the interferon response are upregulated post-therapy along with an increase in non-neoplastic oligodendrocytes. Through in vitro pharmacological perturbation, we demonstrated novel malignant cell-intrinsic targets. This multiomic atlas of longitudinal pHGG captures the key features of therapy response that support distinction from its adult counterpart and suggests therapeutic strategies which are targeted to pediatric gliomas.

The dark side of stemness - the role of hematopoietic stem cells in development of blood malignancies

Hematopoietic stem cells (HSCs) produce all blood cells throughout the life of the organism. However, the high self-renewal and longevity of HSCs predispose them to accumulate mutations. The acquired mutations drive preleukemic clonal hematopoiesis, which is frequent among elderly people. The preleukemic state, although often asymptomatic, increases the risk of blood cancers. Nevertheless, the direct role of preleukemic HSCs is well-evidenced in adult myeloid leukemia (AML), while their contribution to other hematopoietic malignancies remains less understood. Here, we review the evidence supporting the role of preleukemic HSCs in different types of blood cancers, as well as present the alternative models of malignant evolution. Finally, we discuss the clinical importance of preleukemic HSCs in choosing the therapeutic strategies and provide the perspective on further studies on biology of preleukemic HSCs.

Scalable single-cell profiling of chromatin modifications with sciCUT&Tag

Cleavage under targets and tagmentation (CUT&Tag) is an antibody-directed in situ chromatin profiling strategy that is rapidly replacing immune precipitation-based methods, such as chromatin immunoprecipitation-sequencing. The efficiency of the method enables chromatin profiling in single cells but is limited by the numbers of cells that can be profiled. Here, we describe a combinatorial barcoding strategy for CUT&Tag that harnesses a nanowell dispenser for simple, high-resolution, high-throughput, single-cell chromatin profiling. In this single-cell combinatorial indexing CUT&Tag (sciCUT&Tag) protocol, lightly cross-linked nuclei are bound to magnetic beads and incubated with primary and secondary antibodies in bulk and then arrayed in a 96-well plate for a first round of cellular indexing by antibody-directed Tn5 tagmentation. The sample is then repooled, mixed and arrayed across 5,184 nanowells at a density of 12-24 nuclei per well for a second round of cellular indexing during PCR amplification of the sequencing-ready library. This protocol can be completed in 1.5 days by a research technician, and we illustrate the optimized protocol by profiling histone modifications associated with developmental gene repression (H3K27me3) as well as transcriptional activation (H3K4me1-2-3) in human peripheral blood mononuclear cells and use single-nucleotide polymorphisms to facilitate collision removal. We have also used sciCUT&Tag for simultaneous profiling of multiple chromatin epitopes in single cells. The reduced cost, improved resolution and scalability of sciCUT&Tag make it an attractive platform to profile chromatin features in single cells.

INSPIRED Symposium Part 4A: Access to CAR T Cell Therapy in Unique Populations with B Cell Acute Lymphoblastic Leukemia

The approval of tisagenlecleucel (tisa-cel) for use in children with B cell acute lymphoblastic leukemia (B-ALL) was based on the phase 2 ELIANA trial, a global registration study. However, the ELIANA trial excluded specific subsets of patients facing unique challenges and did not include a sufficient number of patients to adequately evaluate outcomes in rare subpopulations. Since the commercialization of tisa-cel, data have become available that support therapeutic indications beyond the specific cohorts previously eligible for chimeric antigen receptor (CAR) T cells targeted to CD19 (CD19 CAR-T) therapy on the registration clinical trial. Substantial real-world data and aggregate clinical trial data have addressed gaps in our understanding of response rates, longer-term efficacy, and toxicities associated with CD19 CAR-T in special populations and rare clinical scenarios. These include patients with central nervous system relapsed disease, who were excluded from ELIANA and other early CAR-T trials owing to concerns about risk of neurotoxicity that have not been born out. There is also interest in the use of CD19 CAR-T for very-high-risk patients earlier in the course of therapy, such as patients with persistent minimal residual disease after 2 cycles of upfront chemotherapy and patients with first relapse of B-ALL. However, these indications are not specified on the label for tisa-cel and historically were not included in eligibility criteria for most clinical trials; data addressing these populations are needed. Populations at high risk of relapse, including patients with high-risk cytogenetic lesions, infants with B-ALL, patients with trisomy 21, and young adults with B-ALL, also may benefit from earlier treatment with CD19 CAR-T. It is important to prospectively study patient-reported outcomes given the differential toxicity expected between CD19 CAR-T and the historic standard of care, hematopoietic cell transplantation. Now that CD19 CAR-T therapy is commercially available, studies evaluating potential access disparities created by this very expensive novel therapy are increasingly pressing.

A longitudinal single-cell atlas of treatment response in pediatric AML

Pediatric acute myeloid leukemia (pAML) is characterized by heterogeneous cellular composition, driver alterations and prognosis. Characterization of this heterogeneity and how it affects treatment response remains understudied in pediatric patients. We used single-cell RNA sequencing and single-cell ATAC sequencing to profile 28 patients representing different pAML subtypes at diagnosis, remission and relapse. At diagnosis, cellular composition differed between genetic subgroups. Upon relapse, cellular hierarchies transitioned toward a more primitive state regardless of subtype. Primitive cells in the relapsed tumor were distinct compared to cells at diagnosis, with under-representation of myeloid transcriptional programs and over-representation of other lineage programs. In some patients, this was accompanied by the appearance of a B-lymphoid-like hierarchy. Our data thus reveal the emergence of apparent subtype-specific plasticity upon treatment and inform on potentially targetable processes.

Genetic alterations and MRD refine risk assessment for KMT2A-rearranged B-cell precursor ALL in adults: a GRAALL study

KMT2A-rearranged (KMT2A-r) B-cell precursor acute lymphoblastic leukemia (BCP-ALL) is widely recognized as a high-risk leukemia in both children and adults. However, there is a paucity of data on adults treated in recent protocols, and the optimal treatment strategy for these patients is still a matter of debate. In this study, we set out to refine the prognosis of adult KMT2A-r BCP-ALL treated with modern chemotherapy regimen and investigate the prognostic impact of comutations and minimal residual disease (MRD). Of 1091 adult patients with Philadelphia-negative BCP-ALL enrolled in 3 consecutive trials from the Group for Research on Adult Acute Lymphoblastic Leukemia (GRAALL), 141 (12.9%) had KMT2A-r, with 5-year cumulative incidence of relapse (CIR) and overall survival (OS) rates of 40.7% and 53.3%, respectively. Molecular profiling highlighted a low mutational burden in this subtype, reminiscent of infant BCP-ALL. However, the presence of TP53 and/or IKZF1 alterations defined a subset of patients with significantly poorer CIR (69.3% vs 36.2%; P = .001) and OS (28.1% vs 60.7%; P = .006) rates. Next, we analyzed the prognostic implication of MRD measured after induction and first consolidation, using both immunoglobulin (IG) or T-cell receptor (TR) gene rearrangements and KMT2A genomic fusion as markers. In approximately one-third of patients, IG/TR rearrangements were absent or displayed clonal evolution during the disease course, compromising MRD monitoring. In contrast, KMT2A-based MRD was highly reliable and strongly associated with outcome, with early good responders having an excellent outcome (3-year CIR, 7.1%; OS, 92.9%). Altogether, our study reveals striking heterogeneity in outcomes within adults with KMT2A-r BCP-ALL and provides new biomarkers to guide risk-based therapeutic stratification.

Modelling acquired resistance to DOT1L inhibition exhibits the adaptive potential of KMT2A-rearranged acute lymphoblastic leukemia

In KMT2A-rearranged acute lymphoblastic leukemia (ALL), an aggressive malignancy, oncogenic KMT2A-fusion proteins inappropriately recruit DOT1L to promote leukemogenesis, highlighting DOT1L as an attractive therapeutic target. Unfortunately, treatment with the first-in-class DOT1L inhibitor pinometostat eventually leads to non-responsiveness. To understand this we established acquired pinometostat resistance in pediatric KMT2A::AFF1+ B-ALL cells. Interestingly, these cells became mostly independent of DOT1L-mediated H3K79 methylation, but still relied on the physical presence of DOT1L, HOXA9 and the KMT2A::AFF1 fusion. Moreover, these cells selectively lost the epigenetic regulation and expression of various KMT2A-fusion target genes such as PROM1/CD133, while other KMT2A::AFF1 target genes, including HOXA9 and CDK6 remained unaffected. Concomitantly, these pinometostat-resistant cells showed upregulation of several myeloid-associated genes, including CD33 and LILRB4/CD85k. Taken together, this model comprehensively shows the adaptive potential of KMT2A-rearranged ALL cells upon losing dependency on one of its main oncogenic properties.

The Gene Expression Classifier ALLCatchR Identifies B-cell Precursor ALL Subtypes and Underlying Developmental Trajectories Across Age

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.

Clinical and translational relevance of intratumor heterogeneity

Intratumor heterogeneity (ITH) is a driver of tumor evolution and a main cause of therapeutic resistance. Despite its importance, measures of ITH are still not incorporated into clinical practice. Consequently, standard treatment is frequently ineffective for patients with heterogeneous tumors as changes to treatment regimens are made only after recurrence and disease progression. More effective combination therapies require a mechanistic understanding of ITH and ways to assess it in clinical samples. The growth of technologies enabling the spatially intact analysis of tumors at the single-cell level and the development of sophisticated preclinical models give us hope that ITH will not simply be used as a predictor of a poor outcome but will guide treatment decisions from diagnosis through treatment.

Single-cell chromatin accessibility profiling of acute myeloid leukemia reveals heterogeneous lineage composition upon therapy-resistance

Acute myeloid leukemia (AML) is a heterogeneous disease characterized by high rate of therapy resistance. Since the cell of origin can impact response to therapy, it is crucial to understand the lineage composition of AML cells at time of therapy resistance. Here we leverage single-cell chromatin accessibility profiling of 22 AML bone marrow aspirates from eight patients at time of therapy resistance and following subsequent therapy to characterize their lineage landscape. Our findings reveal a complex lineage architecture of therapy-resistant AML cells that are primed for stem and progenitor lineages and spanning quiescent, activated and late stem cell/progenitor states. Remarkably, therapy-resistant AML cells are also composed of cells primed for differentiated myeloid, erythroid and even lymphoid lineages. The heterogeneous lineage composition persists following subsequent therapy, with early progenitor-driven features marking unfavorable prognosis in The Cancer Genome Atlas AML cohort. Pseudotime analysis further confirms the vast degree of heterogeneity driven by the dynamic changes in chromatin accessibility. Our findings suggest that therapy-resistant AML cells are characterized not only by stem and progenitor states, but also by a continuum of differentiated cellular lineages. The heterogeneity in lineages likely contributes to their therapy resistance by harboring different degrees of lineage-specific susceptibilities to therapy.

Layered immunity and layered leukemogenicity: Developmentally restricted mechanisms of pediatric leukemia initiation

Hematopoietic stem cells (HSCs) and multipotent progenitor cells (MPPs) arise in successive waves during ontogeny, and their properties change significantly throughout life. Ontological changes in HSCs/MPPs underlie corresponding changes in mechanisms of pediatric leukemia initiation. As HSCs and MPPs progress from fetal to neonatal, juvenile and adult stages of life, they undergo transcriptional and epigenetic reprogramming that modifies immune output to meet age-specific pathogenic challenges. Some immune cells arise exclusively from fetal HSCs/MPPs. We propose that this layered immunity instructs cell fates that underlie a parallel layered leukemogenicity. Indeed, some pediatric leukemias, such as juvenile myelomonocytic leukemia, myeloid leukemia of Down syndrome, and infant pre-B-cell acute lymphoblastic leukemia, are age-restricted. They only present during infancy or early childhood. These leukemias likely arise from fetal progenitors that lose competence for transformation as they age. Other childhood leukemias, such as non-infant pre-B-cell acute lymphoblastic leukemia and acute myeloid leukemia, have mutation profiles that are common in childhood but rare in morphologically similar adult leukemias. These differences could reflect temporal changes in mechanisms of mutagenesis or changes in how progenitors respond to a given mutation at different ages. Interactions between leukemogenic mutations and normal developmental switches offer potential targets for therapy.

The EGR3 regulome of infant KMT2A-r acute lymphoblastic leukemia identifies differential expression of B-lineage genes predictive for outcome

KMT2A-rearranged acute lymphoblastic infant leukemia (KMT2A-r iALL) is associated with outsize risk of relapse and relapse mortality. We previously reported strong upregulation of the immediate early gene EGR3 in KMT2A::AFF1 iALL at relapse; now we provide analyses of the EGR3 regulome, which we assessed through binding and expression target analysis of an EGR3-overexpressing t(4;11) cell culture model. Our data identify EGR3 as a regulator of early B-lineage commitment. Principal component analysis of 50 KMT2A-r iALL patients at diagnosis and 18 at relapse provided strictly dichotomous separation of patients based on the expression of four B-lineage genes. Absence of B-lineage gene expression translates to more than two-fold poorer long-term event-free survival. In conclusion, our study presents four B-lineage genes with prognostic significance, suitable for gene expression-based risk stratification of KMT2A-r iALL patients.

Phenotypic and genotypic infidelity in B-lineage neoplasms, including transdifferentiation following targeted therapy: Report from the 2021 SH/EAHP Workshop

OBJECTIVES

Session 2 of the 2021 Society for Hematopathology and European Association for Haematopathology Workshop collected examples of lineage infidelity and transdifferentiation in B-lineage neoplasms, including after targeted therapy.

METHODS

Twenty cases were submitted. Whole-exome sequencing and genome-wide RNA expression analysis were available on a limited subsample.

RESULTS

A diagnosis of B-cell acute lymphoblastic leukemia (B-ALL) was rendered on at least 1 biopsy from 13 patients. There was 1 case of acute myeloid leukemia (AML); the remaining 6 cases were mature B-cell neoplasms. Targeted therapy was administered in 7 cases of B-ALL and 4 cases of mature B-cell neoplasms. Six cases of B-ALL underwent lineage switch to AML or mixed-phenotype acute leukemia at relapse, 5 of which had rearranged KMT2A. Changes in maturational state without lineage switch were observed in 2 cases. Examples of de novo aberrant T-cell antigen expression (n = 2) were seen among the mature B-cell lymphoma cohort, and their presence correlated with alterations in tumor cell gene expression patterns.

CONCLUSIONS

This cohort of cases enabled us to illustrate, discuss, and review current concepts of lineage switch and aberrant antigen expression in a variety of B-cell neoplasms and draw attention to the role targeted therapies may have in predisposing neoplasms to transdifferentiation as well as other, less expected changes in maturational status.

Single-cell analysis of acute lymphoblastic and lineage-ambiguous leukemia: approaches and molecular insights

Despite recent progress in identifying the genetic drivers of acute lymphoblastic leukemia (ALL), prognosis remains poor for those individuals who experience disease recurrence. Moreover, acute leukemias of ambiguous lineage lack a biologically informed framework to guide classification and therapy. These needs have driven the adoption of multiple complementary single-cell sequencing approaches to explore key issues in the biology of these leukemias, including cell of origin, developmental hierarchy and ontogeny, and the molecular heterogeneity driving pathogenesis, progression, and therapeutic responsiveness. There are multiple single-cell techniques for profiling a specific modality, including RNA, DNA, chromatin accessibility and methylation; and an expanding range of approaches for simultaneous analysis of multiple modalities. Single-cell sequencing approaches have also enabled characterization of cell-intrinsic and -extrinsic features of ALL biology. In this review we describe these approaches and highlight the extensive heterogeneity that underpins ALL gene expression, cellular differentiation, and clonal architecture throughout disease pathogenesis and treatment resistance. In addition, we discuss the importance of the dynamic interactions that occur between leukemia cells and the nonleukemia microenvironment. We discuss potential opportunities and limitations of single-cell sequencing for the study of ALL biology and treatment responsiveness.

Preclinical efficacy of azacitidine and venetoclax for infant KMT2A-rearranged acute lymphoblastic leukemia reveals a new therapeutic strategy

Infants with KMT2A-rearranged B-cell acute lymphoblastic leukemia (ALL) have a dismal prognosis. Survival outcomes have remained static in recent decades despite treatment intensification and novel therapies are urgently required. KMT2A-rearranged infant ALL cells are characterized by an abundance of promoter hypermethylation and exhibit high BCL-2 expression, highlighting potential for therapeutic targeting. Here, we show that hypomethylating agents exhibit in vitro additivity when combined with most conventional chemotherapeutic agents. However, in a subset of samples an antagonistic effect was seen between several agents. This was most evident when hypomethylating agents were combined with methotrexate, with upregulation of ATP-binding cassette transporters identified as a potential mechanism. Single agent treatment with azacitidine and decitabine significantly prolonged in vivo survival in KMT2A-rearranged infant ALL xenografts. Treatment of KMT2A-rearranged infant ALL cell lines with azacitidine and decitabine led to differential genome-wide DNA methylation, changes in gene expression and thermal proteome profiling revealed the target protein-binding landscape of these agents. The selective BCL-2 inhibitor, venetoclax, exhibited in vitro additivity in combination with hypomethylating or conventional chemotherapeutic agents. The addition of venetoclax to azacitidine resulted in a significant in vivo survival advantage indicating the therapeutic potential of this combination to improve outcome for infants with KMT2A-rearranged ALL.

Dissecting the tumor microenvironment of epigenetically driven gliomas: Opportunities for single-cell and spatial multiomics

Malignant gliomas are incurable brain neoplasms with dismal prognoses and near-universal fatality, with minimal therapeutic progress despite billions of dollars invested in research and clinical trials over the last 2 decades. Many glioma studies have utilized disparate histologic and genomic platforms to characterize the stunning genomic, transcriptomic, and immunologic heterogeneity found in gliomas. Single-cell and spatial omics technologies enable unprecedented characterization of heterogeneity in solid malignancies and provide a granular annotation of transcriptional, epigenetic, and microenvironmental states with limited resected tissue. Heterogeneity in gliomas may be defined, at the broadest levels, by tumors ostensibly driven by epigenetic alterations (IDH- and histone-mutant) versus non-epigenetic tumors (IDH-wild type). Epigenetically driven tumors are defined by remarkable transcriptional programs, immunologically distinct microenvironments, and incompletely understood topography (unique cellular neighborhoods and cell-cell interactions). Thus, these tumors are the ideal substrate for single-cell multiomic technologies to disentangle the complex intra-tumoral features, including differentiation trajectories, tumor-immune cell interactions, and chromatin dysregulation. The current review summarizes the applications of single-cell multiomics to existing datasets of epigenetically driven glioma. More importantly, we discuss future capabilities and applications of novel multiomic strategies to answer outstanding questions, enable the development of potent therapeutic strategies, and improve personalized diagnostics and treatment via digital pathology.

Updates in infant acute lymphoblastic leukemia and the potential for targeted therapy

Outcomes for infants diagnosed under 1 year of age with KMT2A-rearranged acute lymphoblastic leukemia (ALL) have remained stagnant over the past 20 years. Successive treatment protocols have previously focused on intensification of conventional chemotherapy, but increased treatment-related toxicity and chemoresistance have led to a plateau in survival. We have now entered an era of immunotherapy with integration of agents, such as blinatumomab or chimeric antigen receptor T-cell therapy, into the standard chemotherapy backbone, showing significant promise for improving the dismal outcomes for this disease. There remains much optimism for the future as a wealth of preclinical studies have identified additional novel targeted agents, such as venetoclax or menin inhibitors, ready for incorporation into treatment, providing further ammunition to combat this aggressive disease. In contrast, infants with KMT2A-germline ALL have demonstrated excellent survival outcomes with current therapy, but there remains a high burden of treatment-related morbidity. Greater understanding of the underlying blast genetics for infants with KMT2A-germline ALL and incorporation of immunotherapeutic approaches may enable a reduction in the intensity of chemotherapy while maintaining the excellent outcomes.

Single-cell multiomics sequencing reveals the reprogramming defects in embryos generated by round spermatid injection

Round spermatid injection (ROSI) technique holds great promise for clinical treatment of a proportion of infertile men. However, the compromised developmental potential of ROSI embryos largely limits the clinical application, and the mechanisms are not fully understood. Here, we describe the transcriptome, chromatin accessibility, and DNA methylation landscapes of mouse ROSI embryos derived from early-stage round spermatids using a single-cell multiomics sequencing approach. By interrogating these data, we identify the reprogramming defects in ROSI embryos at the pronuclear stages, which are mainly associated with the misexpression of a cohort of minor zygotic genome activation genes. We screen a small compound, A366, that can significantly increase the developmental potential of ROSI embryos, in which A366 can partially overcome the reprogramming defects by amending the epigenetic and transcriptomic states. Collectively, our study uncovers the reprogramming defects in ROSI embryos for understanding the mechanisms underlying compromised developmental potential and offers an avenue for ROSI technique optimization.

Poor treatment responses were related to poor outcomes in pediatric B cell acute lymphoblastic leukemia with KMT2A rearrangements

Background

The KMT2A gene, formerly named the MLL gene, is rearranged (KMT2Ar) in 70–75% of infants, 5–6% of children and 10–15% of adult patients with B cell acute lymphoblastic leukemia (B-ALL). The outcome after chemotherapy of pediatric cases remains poor, and only a few studies have investigated the clinical and laboratory features, treatment response and prognosis in Chinese populations.

Methods

A total of 48 B-ALL children with KMT2Ar were enrolled in the study, and clinical and laboratory data were collected and analyzed by age group. The relationship between prognosis and traditional risk factors and treatment response was investigated for these patients who received chemotherapy.

Results

The 48 enrolled patients included 28 males and 20 females; 18 (37.50%) or 30 (62.50%) patients were an age of < 12 m (infant B-ALL) or of > 12 m at onset. An initial WBC count of 300 × 10⁹/L was detected in 7 (14.58%) patients; testicular leukemia (TL) or central nervous system involvement was found in 5 (10.41%) or 3 (6.25%) patients, respectively. Statistical differences were not found in the age groups of sex or initial WBC count, whereas TL was more common in the infant group (P < 0.05). 11q23 was detected in 18 patients; KMT2Ar was detected in 46 (95.83%) or 45 (93.75%) patients by FISH or multiplex RT–PCR technology, respectively; RNA-seq data were obtained for 18 patients, and 3 patients with uncommon KMT2Ar were identified. KMT2A-AFF1, KMT2A-MLLT3 and KMT2A-MLLT1 were the most common transcripts. Statistical differences were not found in treatment response by age groups, including dexamethasone induction, bone marrow (BM) smear status and minimal residual disease (MRD) level at different time points (TP), treatment-related mortality (TRM), or complete remission (CR) rate (P > 0.05); MRD levels monitored by FCM or PCR were unequal at the same TP. Four patients died of treatment, and TRM was 8.33%; 40 patients achieved CR, and the CR rate for the cohort was 83.33%. Seven patients quit, 15 patients relapsed, and the 5 yr cumulative relapse rate was 59.16 ± 9.16%; the 5 yr prospective EFS (pEFS) for patients who were included or excluded from the TRM group was 36.86 ± 8.48% or 40.84 ± 9.16%, respectively. Multivariate analysis for prognosis and hazard ratio was performed for 37 patients without TRM and revealed that an initial WBC count of > 300 × 10⁹/L and a positive level of FCM-MRD were strongly related to a poor outcome for B-ALL patients with KMT2Ar (P < 0.05).

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09804-w.

Elevated Mutational Age in Blood of Children Treated for Cancer Contributes to Therapy-Related Myeloid Neoplasms

Childhood cancer survivors are confronted with various chronic health conditions like therapy-related malignancies. However, it is unclear how exposure to chemotherapy contributes to the mutation burden and clonal composition of healthy tissues early in life. Here, we studied mutation accumulation in hematopoietic stem and progenitor cells (HSPC) before and after cancer treatment of 24 children. Of these children, 19 developed therapy-related myeloid neoplasms (t-MN). Posttreatment HSPCs had an average mutation burden increase comparable to what treatment-naïve cells accumulate during 16 years of life, with excesses up to 80 years. In most children, these additional mutations were induced by clock-like processes, which are also active during healthy aging. Other patients harbored mutations that could be directly attributed to treatments like platinum-based drugs and thiopurines. Using phylogenetic inference, we demonstrate that most t-MN in children originate after the start of treatment and that leukemic clones become dominant during or directly after chemotherapy exposure.

SIGNIFICANCE

Our study shows that chemotherapy increases the mutation burden of normal blood cells in cancer survivors. Only few drugs damage the DNA directly, whereas in most patients, chemotherapy-induced mutations are caused by processes similar to those present during normal aging. This article is highlighted in the In This Issue feature, p. 1825.

To B- or not to B-: A review of lineage switched acute leukemia

Acute leukemia is a heterogeneous disorder of hematologic malignancies composed primarily of hematopoietic precursors that have acquired unregulated self-renewal and proliferation. Hematology classification systems typically divide these neoplasms into lymphoid (B- or T-) and myeloid-lineage subtypes, with therapy dependent upon this distinction. Infrequently, certain acute leukemias may undergo a complete lineage switch at relapse, subsequently complicating the diagnosis and treatment of these recurrent diseases. Transformation from B-lineage to myeloid lineage is the most common switch observed, and is frequently associated with a balanced 11q23 translocation, involving KMT2A. The mechanisms involved in the lineage-switch are unclear, but modern therapies targeting the B-cell-specific marker, CD19, have proven to promote this conversion as one means of treatment escape. Broadly speaking, therapy-mediated selection of alternate lineage-committed subclones derived from the same initial pluripotent progenitors, clonal evolution and reprogramming of lineage-committed blasts, and de novo clonally unrelated leukemias may account for the clinical impression of lineage switched acute leukemia during treatment. This review will explore the phenomenon and potential mechanisms of lineage transformation during the treatment of acute leukemia.

Neonatal congenital leukemia caused by several missense mutations and AFF1‑KMT2A fusion: A case report

Neonatal leukemia, a congenital form of leukemia, is a rare and fatal disease occurring in the neonatal period. Its etiology and pathogenesis have remained to be fully elucidated and the clinical manifestations differ due to age variability. Acute myeloid leukemia (AML) occurring after birth indicates genetic abnormalities and possibly intrauterine exposure to radiation, drugs or other toxins. The present report described the case of a premature neonate without phenotypic signs of Down syndrome, but with an elevated white blood cell count, mainly pertaining to the monocytes of peripheral blood. At 31 weeks of gestation, delivery by Caesarean section was performed due to fetal distress; however, the infant died three days after birth. Further laboratory examination indicated pediatric myeloid leukemia. The present case report described a case of fetal AML. According to the results of peripheral blood smear and targeted-panel sequencing, 5 missense mutations with clinical significance and a novel AFF1-KMT2A fusion gene were detected, which may be the main causes of AML and death.

Hypoxic, glycolytic metabolism is a vulnerability of B-acute lymphoblastic leukemia-initiating cells

High-risk forms of B-acute lymphoblastic leukemia (B-ALL) remain a therapeutic challenge. Leukemia-initiating cells (LICs) self-renew and spark relapse and therefore have been the subject of intensive investigation; however, the properties of LICs in high-risk B-ALL are not well understood. Here, we use single-cell transcriptomics and quantitative xenotransplantation to understand LICs in MLL-rearranged (MLL-r) B-ALL. Compared with reported LIC frequencies in acute myeloid leukemia (AML), engraftable LICs in MLL-r B-ALL are abundant. Although we find that multipotent, self-renewing LICs are enriched among phenotypically undifferentiated B-ALL cells, LICs with the capacity to replenish the leukemic cellular diversity can emerge from more mature fractions. While inhibiting oxidative phosphorylation blunts blast proliferation, this intervention promotes LIC emergence. Conversely, inhibiting hypoxia and glycolysis impairs MLL-r B-ALL LICs, providing a therapeutic benefit in xenotransplantation systems. These findings provide insight into the aggressive nature of MLL-r B-ALL and provide a rationale for therapeutic targeting of hypoxia and glycolysis.