Single-Cell Proteogenomics Analysis of AML Samples Treated with Decitabine and Venetoclax Reveals Divergent Treatment Escape Mechanisms

Background: The combination therapy of hypomethylating agents and venetoclax (HMA+VEN) has emerged as a standard of care for acute myeloid leukemia (AML) patients unfit for intensive chemotherapy. Nonetheless, significant clinical challenges persist, including relapse and treatment resistance. Previous studies have revealed genetic ( RAS/RTK pathway mutations and TP53 mutations) and phenotypic (monocytic differentiation) resistance mechanisms in the context of HMA+VEN therapy. Despite these advancements, our comprehensive understanding of the HMA+VEN resistance mechanisms remains incomplete. To gain further insights into the resistance mechanism in HMA+VEN therapy, we implemented single-cell proteogenomics analysis on longitudinally collected samples from AML patients treated with decitabine+venetoclax (DAC+VEN) trial. Methods: We performed a high-throughput, single-cell targeted DNA-antibody sequencing (DAb-seq) using Tapestri platform from Mission Bio Inc. Sixty-three AML samples sequentially collected from 33 unique patients who participated in DAC+VEN clinical trial (NCT03404193) were analyzed. Our custom DNA panel covered 37 recurrently mutated genes in AML, whereas the custom antibody panel included either 15 cell surface proteins (N = 20 patients) or 48 proteins (N = 13 patients). Cellular phenotypes were inferred based on the patterns of cell surface protein expression. Median 2,767 cells were sequenced per sample, with total 196,719 cells sequenced for this study. Results: The median age of the cohort was 70 years (range 21-84). Twenty-four (73%) patients had primary AML and 9 (27%) had secondary or therapy-related AML. Based on the FAB classification, 13 patients had M0/M1/M2, 6 had M4/M5, one had M7, and the rest had unknown subtypes. The most frequent gene mutation at baseline was DNMT3A (45%), followed by NPM1 (39%), TET2 (30%), and TP53 (24%). Among the 33 patients, 32 had evaluable treatment response; 27 (82%) were defined as responders by achieving CR or CRi, whereas 5 (15%) were non-responders. Six patients experienced relapse after responding to the therapy, and their baseline-relapse pairs were analyzed by the single-cell platform. Using the composition of cellular phenotype data, AML samples were classified into either stem-like (N = 5), progenitor-like (N = 1), predominantly monocytic (N = 4), predominantly erythroid (N = 2) and mixed population AML (mixed population of stem, progenitor, monocyte or erythroid-like population, N = 21).