Cancer cell heterogeneity is a major driver of therapy resistance. To characterize resistant cells and their vulnerabilities, we studied the PLZF-RARA variant of acute promyelocytic leukemia, resistant to retinoic acid (RA), using single-cell multiomics. We uncovered transcriptional and chromatin heterogeneity in leukemia cells. We identified a subset of cells resistant to RA with proliferation, DNA replication, and repair signatures that depend on a fine-tuned E2F transcriptional network targeting the epigenetic regulator enhancer of zeste homolog 2 (EZH2). Epigenomic and functional analyses validated the driver role of EZH2 in RA resistance. Targeting pan-EZH2 activities (canonical/noncanonical) was necessary to eliminate leukemia relapse-initiating cells, which underlies a dependency of resistant cells on an EZH2 noncanonical activity and the necessity to degrade EZH2 to overcome resistance. Our study provides critical insights into the mechanisms of RA resistance that allow us to eliminate treatment-resistant leukemia cells by targeting EZH2, thus highlighting a potential targeted therapy approach. Beyond RA resistance and acute promyelocytic leukemia context, our study also demonstrates the power of single-cell multiomics to identify, characterize, and clear therapy-resistant cells.