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.