Studies of acute myeloid leukemia rely on DNA sequencing and immunophenotyping by flow cytometry as primary tools for disease characterization. However, leukemia tumor heterogeneity complicates integration of DNA variants and immunophenotypes from separate measurements. Here we introduce DAb-seq, a technology for simultaneous capture of DNA genotype and cell surface phenotype from single cells at high throughput, enabling direct profiling of proteogenomic states in tens of thousands of cells. To demonstrate the approach, we analyze the disease of three patients with leukemia over multiple treatment timepoints and disease recurrences. We observe complex genotype-phenotype dynamics that illustrate the subtlety of the disease process and the degree of incongruity between blast cell genotype and phenotype in different clinical scenarios. Our results highlight the importance of combined single-cell DNA and protein measurements to fully characterize the heterogeneity of leukemia.

Single-cell sequencing exposes previously unmeasurable complexity of tumor heterogeneity and clonal evolution on quizartinib. Single-cell sequencing reveals on- and off-target mechanisms of resistance to quizartinib, which can preexist therapy.

Mutations in the receptor tyrosine kinases (RTKs) FLT3 and KIT are frequent and associated with poor outcomes in acute myeloid leukemia (AML). Although selective FLT3 inhibitors (FLT3i) are clinically effective, remissions are short-lived due to secondary resistance characterized by acquired mutations constitutively activating the RAS/MAPK pathway. Hereby, we report the pre-clinical efficacy of co-targeting SHP2, a critical node in MAPK signaling, and BCL2 in RTK-driven AML. The allosteric SHP2 inhibitor RMC-4550 suppresses proliferation of AML cell lines with FLT3 and KIT mutations, including cell lines with acquired resistance to FLT3i. We demonstrate that pharmacologic SHP2 inhibition unveils an Achilles' heel of RTK-driven AML, increasing apoptotic dependency on BCL2 via MAPK-dependent mechanisms, including upregulation of BMF and downregulation of MCL1. Consequently, RMC-4550 and venetoclax are synergistically lethal in AML cell lines and in clinically relevant xenograft models. Our results provide mechanistic rationale and pre-clinical evidence for co-targeting SHP2 and BCL2 in RTK-driven AML.

Mutations in five canonical Ras pathway genes (NF1, NRAS, KRAS, PTPN11 and CBL) are detected in nearly 90% of patients with juvenile myelomonocytic leukemia (JMML), a frequently fatal malignant neoplasm of early childhood. In this report, we describe seven patients diagnosed with SH2B3-mutated JMML, including five patients who were found to have initiating, loss-of-function mutations in the gene. SH2B3 encodes the adaptor protein LNK, a negative regulator of normal hematopoiesis upstream of the Ras pathway. These mutations were identified to be germline, somatic or a combination of both. Loss of function of LNK, which has been observed in other myeloid malignancies, results in abnormal proliferation of hematopoietic cells due to cytokine hypersensitivity and activation of the JAK/STAT signaling pathway. In vitro studies of induced pluripotent stem cell-derived JMML-like hematopoietic progenitor cells also demonstrated sensitivity of SH2B3-mutated hematopoietic progenitor cells to JAK inhibition. Lastly, we describe two patients with JMML and SH2B3 mutations who were treated with the JAK1/2 inhibitor ruxolitinib. This report expands the spectrum of initiating mutations in JMML and raises the possibility of targeting the JAK/STAT pathway in patients with SH2B3 mutations.

Despite recent work linking mixed phenotype acute leukemia (MPAL) to certain genetic lesions, specific driver mutations remain undefined for a significant proportion of patients and no genetic subtype is predictive of clinical outcomes. Moreover, therapeutic strategy for MPAL remains unclear, and prognosis is overall poor. We performed multiomic single cell profiling of 14 newly diagnosed adult MPAL patients to characterize the inter- and intra-tumoral transcriptional, immunophenotypic, and genetic landscapes of MPAL. We show that neither genetic profile nor transcriptome reliably correlate with specific MPAL immunophenotypes. Despite this, we find that MPAL blasts express a shared stem cell-like transcriptional profile indicative of high differentiation potential. Patients with the highest differentiation potential demonstrate inferior survival in our dataset. A gene set score, MPAL95, derived from genes highly enriched in the most stem-like MPAL cells, is applicable to bulk RNA sequencing data and is predictive of survival in an independent patient cohort, suggesting a potential strategy for clinical risk stratification.

Current single-cell RNA-sequencing approaches have limitations that stem from the microfluidic devices or fluid handling steps required for sample processing. We develop a method that does not require specialized microfluidic devices, expertise or hardware. Our approach is based on particle-templated emulsification, which allows single-cell encapsulation and barcoding of cDNA in uniform droplet emulsions with only a vortexer. Particle-templated instant partition sequencing (PIP-seq) accommodates a wide range of emulsification formats, including microwell plates and large-volume conical tubes, enabling thousands of samples or millions of cells to be processed in minutes. We demonstrate that PIP-seq produces high-purity transcriptomes in mouse-human mixing studies, is compatible with multiomics measurements and can accurately characterize cell types in human breast tissue compared to a commercial microfluidic platform. Single-cell transcriptional profiling of mixed phenotype acute leukemia using PIP-seq reveals the emergence of heterogeneity within chemotherapy-resistant cell subsets that were hidden by standard immunophenotyping. PIP-seq is a simple, flexible and scalable next-generation workflow that extends single-cell sequencing to new applications.