Abstract:Introduction: Pediatric acute myeloid leukemia (AML) is a heterogeneous hematological malignancy characterized by various chromosomal abnormalities and somatic mutations. Investigating the origin of tumorigenesis and deciphering the complex interplay between AML cells and their surrounding tumor immune microenvironment (TIME) is crucial to understand the underlying mechanisms driving disease progression and response to therapy. Recent advancements in single-cell multi-omics technologies enabled us to characterize the transcriptional (GEX) and epigenetic (ATAC) landscapes of both tumors and TIMEs. Leveraging our in-house algorithm, scMINER (single-cell Mutual Information-based Network Engineering Ranger, Ding et al., 2023), this study aims to build a comprehensive cell type atlas in diverse pediatric AML bone marrow samples and construct intra- and inter-cellular networks, enabling the identification of hidden drivers that shape subtype-specific AML-TIME interactions. Methods: Using CD45 and CD33 as selective markers by flow cytometry, we have successfully isolated blasts and enriched for T cell and B cell populations for 10X single-cell multiome profiling across AML molecular subtypes, including UBTF-TD (n=1), NPM1 (n=1), and KMT2A rearrangement (n=1). We interrogated the GEX data using scMINER, including two key steps: (i) mutual information-based clustering analysis (MICA) and (ii) mutual information-based network inference engine (MINIE). MICA characterizes the intrinsic nonlinear similarity of gene expression distributions among cells, enabling a high resolution of clusters. MINIE reverse-engineers cluster-specific intracellular gene networks and infers protein activity of transcription factors (TF) or signaling factors (SIG) drivers based on the expression of its predicted regulon targets in the corresponding cluster. The activity profiles overcome the dropout effects of single-cell RNA-seq data and reflect protein activities across all cells. ArchR is used to analyze ATAC data. CellChat is used to construct and evaluate intercellular networks. Results: AML cells from three samples were divided into four major clusters, and we labeled them as NPM1, KMT2A, UBTF_one, and UBTF_two based on sample subtypes (Fig.1A). MICA managed to obtain distinct AML cell sub-clusters. Comparison with single-cell profiles of normal hematopoietic cells and annotation with marker genes revealed their correspondence with different myeloid developmental stages. Notably, the KMT2A cluster shared a high similarity in cell type proportion with the UBTF_two cluster, and the NPM1 cluster is closer to the UBTF_one cluster. This result demonstrated significant tumor heterogeneity within and across AML molecular subtypes. With the help of intracellular gene networks generated by MINIE, we identified multiple subtype-specific hidden drivers by differential activity analysis. For example, KMT2A and UBTF_two clusters exhibited elevated SAMHD1 levels, consistent with previous GSE data (Zhang et al., 2022). Chromatin accessibility analysis of ATAC data showed increased peaks near SAMHD1 in KMT2A and UBTF_two clusters, further supporting this observation. IRF8 showed a similar expression pattern to SAMHD1, suggesting its context-dependent role across AML subtypes. Accurate annotation of myeloid and immune cells greatly enriches the preciseness of the AML and TIME interaction patterns, which vary across AML molecular subtypes (Fig.1B). Interestingly, cells originated from KMT2A cluster exhibited similar interaction patterns to UBTF_two, while NPM1 and UBTF_one were grouped. This result indicated that different AML molecular subtypes may partially utilize the same mechanisms of tumor-TIME interaction, resulting in similar tumor phenotypes. Conclusions: Accurate clustering and well-defined annotation of myeloid and immune cells enabled us to build a comprehensive cell type atlas in AML samples. The subtype-specific network and the following analysis identified candidate hidden drivers that possibly contribute to the heterogeneity of tumor populations. AML subtypes sharing the same hidden drivers have similar cell type proportions and interaction patterns with TIMEs. This observation implies that AML subtypes share cell type-specific dependencies, which can be therapeutic targets to overcome the refractoriness posed by tumor heterogeneity.

Single-cell epigenetic and clonal analysis decodes disease progression in pediatric acute myeloid leukemia

Transcriptome-wide subtyping of pediatric and adult T cell acute lymphoblastic leukemia in an international study of 707 cases

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

Single-cell analysis reveals altered tumor microenvironments of relapse- and remission-associated pediatric acute myeloid leukemia

Single-cell analysis of pediatric acute myeloid leukemia samples uncovers treatment-resistant stem and mast cells

Subtype-specific 3D Genome Alteration in Acute Myeloid Leukaemia

Epigenomic machinery regulating pediatric AML: Clonal expansion mechanisms, therapies, and future perspectives

Single-cell transcriptomic analysis of the immune microenvironment in pediatric acute leukemia

Revealing the Intratumoral Heterogeneity of Non-Ds Acute Megakaryoblastic Leukemia in Single-Cell Resolution

Single-cell transcriptional mapping reveals genetic and non-genetic determinants of aberrant differentiation in AML

Abstract 6234: Single-cell analysis of state-transition from health to leukemia and the leukemic cell heterogeneities in inv(16) acute myeloid leukemia (AML)

A single-cell survey of cellular hierarchy in acute myeloid leukemia

Profiling chromatin accessibility in pediatric acute lymphoblastic leukemia identifies subtype-specific chromatin landscapes and gene regulatory networks

Clonal evolution of acute myeloid leukemia revealed by high-throughput single-cell genomics

Single-cell RNA sequencing distinctly characterizes the wide heterogeneity in pediatric mixed phenotype acute leukemia

Next-generation sequencing reveals gene mutations landscape and clonal evolution in patients with acute myeloid leukemia

Genomic Profiling of Adult and Pediatric B-cell Acute Lymphoblastic Leukemia.

Epigenetic Modifications and Targeted Therapy in Pediatric Acute Myeloid Leukemia

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

Comprehensive genomic analysis reveals molecular heterogeneity in pediatric ALK-positive anaplastic large cell lymphoma

Dissecting Subtype-Specific Tumor-Time Interactions and Underlying Hidden Drivers in Pediatric Acute Myeloid Leukemia Via Single-Cell Multi-Omics