Single-cell RNA-seq reveals developmental plasticity with coexisting oncogenic states and immune evasion programs in ETP-ALL.

Lineage plasticity and stemness have been invoked as causes of therapy resistance in cancer, because these flexible states allow cancer cells to dedifferentiate and alter their dependencies. We investigated such resistance mechanisms in relapsed/refractory early T-cell progenitor acute lymphoblastic leukemia (ETP-ALL) carrying activating NOTCH1 mutations via full-length single-cell RNA sequencing (scRNA-seq) of malignant and microenvironmental cells. We identified 2 highly distinct stem-like states that critically differed with regard to cell cycle and oncogenic signaling. Fast-cycling stem-like leukemia cells demonstrated Notch activation and were effectively eliminated in patients by Notch inhibition, whereas slow-cycling stem-like cells were Notch independent and rather relied on PI3K signaling, likely explaining the poor efficacy of Notch inhibition in this disease. Remarkably, we found that both stem-like states could differentiate into a more mature leukemia state with prominent immunomodulatory functions, including high expression of the LGALS9 checkpoint molecule. These cells promoted an immunosuppressive leukemia ecosystem with clonal accumulation of dysfunctional CD8+ T cells that expressed HAVCR2, the cognate receptor for LGALS9. Our study identified complex interactions between signaling programs, cellular plasticity, and immune programs that characterize ETP-ALL, illustrating the multidimensionality of tumor heterogeneity. In this scenario, combination therapies targeting diverse oncogenic states and the immune ecosystem seem most promising to successfully eliminate tumor cells that escape treatment through coexisting transcriptional programs.

Implementation of a prostate cancer-specific targeted sequencing panel for credentialing of patient-derived cell lines and genomic characterization of patient samples.

BACKGROUND: Primary and metastatic prostate cancers have low mutation rates and recurrent alterations in a small set of genes, enabling targeted sequencing of prostate cancer-associated genes as an efficient approach to characterizing patient samples (compared to whole-exome and whole-genome sequencing). For example, targeted sequencing provides a flexible, rapid, and cost-effective method for genomic assessment of patient-derived cell lines to evaluate fidelity to initial patient tumor samples. METHODS: We developed a prostate cancer-specific targeted next-generation sequencing (NGS) panel to detect alterations in 62 prostate cancer-associated genes as well as recurring gene fusions with ETS family members, representing the majority of common alterations in prostate cancer. We tested this panel on primary prostate cancer tissues and blood biopsies from patients with metastatic prostate cancer. We generated patient-derived cell lines from primary prostate cancers using conditional reprogramming methods and applied targeted sequencing to evaluate the fidelity of these cell lines to the original patient tumors. RESULTS: The prostate cancer-specific panel identified biologically and clinically relevant alterations, including point mutations in driver oncogenes and ETS family fusion genes, in tumor tissues from 29 radical prostatectomy samples. The targeted panel also identified genomic alterations in cell-free DNA and circulating tumor cells (CTCs) from patients with metastatic prostate cancer, and in standard prostate cancer cell lines. We used the targeted panel to sequence our set of patient-derived cell lines; however, no prostate cancer-specific mutations were identified in the tumor-derived cell lines, suggesting preferential outgrowth of normal prostate epithelial cells. CONCLUSIONS: We evaluated a prostate cancer-specific targeted NGS panel to detect common and clinically relevant alterations (including ETS family gene fusions) in prostate cancer. The panel detected driver mutations in a diverse set of clinical samples of prostate cancer, including fresh-frozen tumors, cell-free DNA, CTCs, and cell lines. Targeted sequencing of patient-derived cell lines highlights the challenge of deriving cell lines from primary prostate cancers and the importance of genomic characterization to credential candidate cell lines. Our study supports that a prostate cancer-specific targeted sequencing panel provides an efficient, clinically feasible approach to identify genetic alterations across a spectrum of prostate cancer samples and cell lines.

BMT CTN Myeloma Intergroup Workshop on Minimal Residual Disease and Immune Profiling: Summary and Recommendations from the Organizing Committee.

The Blood and Marrow Transplant Clinical Trials Network Myeloma Intergroup Workshop on Minimal Residual Disease and Immune Profiling was convened on December 1, 2016 at the American Society of Hematology meeting to discuss the emerging data and technologies for minimal residual disease assessment and immune profiling in myeloma. Particular emphasis was placed on developing strategies to incorporate these techniques into clinical trial design. This document reviews the literature, summarizes the topics discussed in the workshop, and provides recommendations for integration of these techniques into future clinical trial design.

Discovery and prioritization of somatic mutations in diffuse large B-cell lymphoma (DLBCL) by whole-exome sequencing.

To gain insight into the genomic basis of diffuse large B-cell lymphoma (DLBCL), we performed massively parallel whole-exome sequencing of 55 primary tumor samples from patients with DLBCL and matched normal tissue. We identified recurrent mutations in genes that are well known to be functionally relevant in DLBCL, including MYD88, CARD11, EZH2, and CREBBP. We also identified somatic mutations in genes for which a functional role in DLBCL has not been previously suspected. These genes include MEF2B, MLL2, BTG1, GNA13, ACTB, P2RY8, PCLO, and TNFRSF14. Further, we show that BCL2 mutations commonly occur in patients with BCL2/IgH rearrangements as a result of somatic hypermutation normally occurring at the IgH locus. The BCL2 point mutations are primarily synonymous, and likely caused by activation-induced cytidine deaminase-mediated somatic hypermutation, as shown by comprehensive analysis of enrichment of mutations in WRCY target motifs. Those nonsynonymous mutations that are observed tend to be found outside of the functionally important BH domains of the protein, suggesting that strong negative selection against BCL2 loss-of-function mutations is at play. Last, by using an algorithm designed to identify likely functionally relevant but infrequent mutations, we identify KRAS, BRAF, and NOTCH1 as likely drivers of DLBCL pathogenesis in some patients. Our data provide an unbiased view of the landscape of mutations in DLBCL, and this in turn may point toward new therapeutic strategies for the disease.

Dual role of signaling pathways in myeloma requires cell type-specific targeting of ligand-receptor interactions.

ABSTRACT: Although most patients with multiple myeloma respond to treatment initially, therapy resistance develops almost invariably, and only a subset of patients show durable responses to immunomodulatory therapies. Although the immune microenvironment has been extensively studied in patients with myeloma, its composition is currently not used as prognostic markers in clinical routine. We hypothesized that the outcome of immune signaling pathway engagement can be highly variable, depending on which 2 cellular populations participate in this interaction. This would have important prognostic and therapeutic implications, suggesting that it is crucial for immune pathways to be targeted in a specific cellular context. To test this hypothesis, we investigated a cohort of 25 patients with newly diagnosed multiple myeloma. We examined the complex regulatory networks within the immune compartment and their impact on disease progression. Analysis of immune cell composition and expression profiles revealed significant differences in the B-cell compartment associated with treatment response. Transcriptional states in patients with short time to progression demonstrated an enrichment of pathways promoting B-cell differentiation and inflammatory responses, which may indicate immune dysfunction. Importantly, the analysis of molecular interactions within the immune microenvironment highlights the dual role of signaling pathways, which can either be associated with good or poor prognosis depending on the cell types involved. Our findings therefore argue that therapeutic strategies targeting ligand-receptor interactions should take into consideration the composition of the microenvironment and the specific cell types involved in molecular interactions.

Genomics of lymphoid malignancies reveal major activation pathways in lymphocytes.

Breakdown of tolerance leads to autoimmunity due to emergence of autoreactive T or B cell clones. Autoimmune diseases predispose to lymphoid malignancies and lymphoid malignancies, conversely, can manifest as autoimmune diseases. While it has been clear for a long time that a competitive advantage and uncontrolled growth of lymphocytes contribute to the pathogenesis of both lymphoid malignancies as well as autoimmune diseases, the overlap of the underlying mechanisms has been less well described. Next generation sequencing has led to massive expansion of the available genomic data in many diseases over the last five years. These data allow for comparison of the molecular pathogenesis between autoimmune diseases and lymphoid malignancies. Here, we review the similarities between autoimmune diseases and lymphoid malignancies: 1) Both, autoimmune diseases and lymphoid malignancies are characterized by activation of the same T and B cell signaling pathways, and dysregulation of these pathways can occur through genetic or epigenetic events. 2) In both scenarios, clonal and subclonal evolution of lymphocytes contribute to disease. 3) Development of both diseases not only depends on T or B cell intrinsic factors, such as germline or somatic mutations, but also on environmental factors. These include infections, the presence of other immune cells in the microenvironment, and the cytokine milieu. A better mechanistic understanding of the parallels between lymphomagenesis and autoimmunity may help the development of precision treatment strategies with rationally designed therapeutic agents.

Loss of metabolic fitness drives tumor resistance after CAR-NK cell therapy and can be overcome by cytokine engineering.

Chimeric antigen receptor (CAR) engineering of natural killer (NK) cells is promising, with early-phase clinical studies showing encouraging responses. However, the transcriptional signatures that control the fate of CAR-NK cells after infusion and factors that influence tumor control remain poorly understood. We performed single-cell RNA sequencing and mass cytometry to study the heterogeneity of CAR-NK cells and their in vivo evolution after adoptive transfer, from the phase of tumor control to relapse. Using a preclinical model of noncurative lymphoma and samples from a responder and a nonresponder patient treated with CAR19/IL-15 NK cells, we observed the emergence of NK cell clusters with distinct patterns of activation, function, and metabolic signature associated with different phases of in vivo evolution and tumor control. Interaction with the highly metabolically active tumor resulted in loss of metabolic fitness in NK cells that could be partly overcome by incorporation of IL-15 in the CAR construct.

Regulatory Programs of B-cell Activation and Germinal Center Reaction Allow B-ALL Escape from CD19 CAR T-cell Therapy.

Chimeric antigen receptor (CAR) T-cell therapy has led to tremendous successes in the treatment of B-cell malignancies. However, a large fraction of treated patients relapse, often with disease expressing reduced levels of the target antigen. Here, we report that exposing CD19+ B-cell acute lymphoblastic leukemia (B-ALL) cells to CD19 CAR T cells reduced CD19 expression within hours. Initially, CD19 CAR T cells caused clustering of CD19 at the T cell-leukemia cell interface followed by CD19 internalization and decreased CD19 surface expression on the B-ALL cells. CD19 expression was then repressed by transcriptional rewiring. Using single-cell RNA sequencing and single-cell assay for transposase-accessible chromatin using sequencing, we demonstrated that a subset of refractory CD19low cells sustained decreased CD19 expression through transcriptional programs of physiologic B-cell activation and germinal center reaction. Inhibiting B-cell activation programs with the Bruton's tyrosine kinase inhibitor ibrutinib increased the cytotoxicity of CD19 CAR T cells without affecting CAR T-cell viability. These results demonstrate transcriptional plasticity as an underlying mechanism of escape from CAR T cells and highlight the importance of combining CAR T-cell therapy with targeted therapies that aim to overcome this plasticity. See related Spotlight by Zhao and Melenhorst, p. 1040.

Cell-free DNA for the detection of emerging treatment failure in relapsed/ refractory multiple myeloma.

Interrogation of cell-free DNA (cfDNA) represents an emerging approach to non-invasively estimate disease burden in multiple myeloma (MM). Here, we examined low-pass whole genome sequencing (LPWGS) of cfDNA for its predictive value in relapsed/ refractory MM (RRMM). We observed that cfDNA positivity, defined as ≥10% tumor fraction by LPWGS, was associated with significantly shorter progression-free survival (PFS) in an exploratory test cohort of 16 patients who were actively treated on diverse regimens. We prospectively determined the predictive value of cfDNA in 86 samples from 45 RRMM patients treated with elotuzumab, pomalidomide, bortezomib, and dexamethasone in a phase II clinical trial (NCT02718833). PFS in patients with tumor-positive and -negative cfDNA after two cycles of treatment was 1.6 and 17.6 months, respectively (HR 7.6, P < 0.0001). Multivariate hazard modelling confirmed cfDNA as independent risk factor (HR 96.6, P = 6.92e-05). While correlating with serum-free light chains and bone marrow, cfDNA additionally discriminated patients with poor PFS among those with the same response by IMWG criteria. In summary, detectability of MM-derived cfDNA, as a measure of substantial tumor burden with therapy, independently predicts poor PFS and may provide refinement for standard-of-care response parameters to identify patients with poor response to treatment earlier than is currently feasible.

A phase II trial of all-trans retinoic acid (ATRA) in advanced adenoid cystic carcinoma.

BACKGROUND: Effective therapies are lacking for recurrent, metastatic adenoid cystic carcinoma (R/M ACC) and preclinical models suggest retinoic acid agonists inhibit ACC growth. This phase II trial evaluated all-trans retinoic acid (ATRA) as a novel therapy for ACC. METHODS: Patients with R/M ACC (any site) with clinical and/or radiographic progression ≤12 months prior to study entry were eligible. Cohort 1 (CH1) received ATRA 45 mg/m2 split oral daily dosing on days 1-14 of a 28-day cycle; Cohort 2 (CH2) received the same dosing continuously. Primary endpoint was best overall response rate (CR + PR) (RECIST v1.1). Secondary endpoints: safety and progression-free survival (PFS). Exploratory analyses: ATRA impact on MYB expression and genomic predictors of response. RESULTS: Eighteen patients enrolled. There were no responses, but 61% (11/18) had stable disease (SD) and 28% (5/18) progression as best response; 11% (2/18) unevaluable. Median duration of stability: 3.7 months (95%CI, 1.9-3.9). One patient (CH1) remains on drug with SD approaching 1 year. Half of those who received prior VEGFR therapy achieved SD (4/8). At median follow up of 7.9 months, median PFS was 3.2 months (95%CI, 1.8-3.9). N = 1 required dose adjustment; N = 1 came off drug for toxicity. There were no grade 3-4 adverse events. NOTCH1 and PI3K pathway alterations were most frequent. Low MYB protein expression was associated with longer duration of stability on ATRA (P < 0.01). CONCLUSION(S): While the trial did not meet its prespecified response endpoint, ATRA alone or in combination may be a low toxicity treatment for disease growth stabilization in R/M ACC.

Dynamic transcriptional reprogramming leads to immunotherapeutic vulnerabilities in myeloma.

While there is extensive evidence for genetic variation as a basis for treatment resistance, other sources of variation result from cellular plasticity. Using multiple myeloma as an example of an incurable lymphoid malignancy, we show how cancer cells modulate lineage restriction, adapt their enhancer usage and employ cell-intrinsic diversity for survival and treatment escape. By using single-cell transcriptome and chromatin accessibility profiling, we show that distinct transcriptional states co-exist in individual cancer cells and that differential transcriptional regulon usage and enhancer rewiring underlie these alternative transcriptional states. We demonstrate that exposure to standard treatment further promotes transcriptional reprogramming and differential enhancer recruitment while simultaneously reducing developmental potential. Importantly, treatment generates a distinct complement of actionable immunotherapy targets, such as CXCR4, which can be exploited to overcome treatment resistance. Our studies therefore delineate how to transform the cellular plasticity that underlies drug resistance into immuno-oncologic therapeutic opportunities.

Single-Cell Profiling Reveals Metabolic Reprogramming as a Resistance Mechanism in BRAF-Mutated Multiple Myeloma.

PURPOSE: Although remarkably effective in some patients, precision medicine typically induces only transient responses despite initial absence of resistance-conferring mutations. Using BRAF-mutated myeloma as a model for resistance to precision medicine we investigated if BRAF-mutated cancer cells have the ability to ensure their survival by rapidly adapting to BRAF inhibitor treatment. EXPERIMENTAL DESIGN: Full-length single-cell RNA (scRNA) sequencing (scRNA-seq) was conducted on 3 patients with BRAF-mutated myeloma and 1 healthy donor. We sequenced 1,495 cells before, after 1 week, and at clinical relapse to BRAF/MEK inhibitor treatment. We developed an in vitro model of dabrafenib resistance using genetically homogeneous single-cell clones from two cell lines with established BRAF mutations (U266, DP6). Transcriptional and epigenetic adaptation in resistant cells were defined by RNA-seq and H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq). Mitochondrial metabolism was characterized by metabolic flux analysis. RESULTS: Profiling by scRNA-seq revealed rapid cellular state changes in response to BRAF/MEK inhibition in patients with myeloma and cell lines. Transcriptional adaptation preceded detectable outgrowth of genetically discernible drug-resistant clones and was associated with widespread enhancer remodeling. As a dominant vulnerability, dependency on oxidative phosphorylation (OxPhos) was induced. In treated individuals, OxPhos was activated at the time of relapse and showed inverse correlation to MAPK activation. Metabolic flux analysis confirmed OxPhos as a preferential energetic resource of drug-persistent myeloma cells. CONCLUSIONS: This study demonstrates that cancer cells have the ability to rapidly adapt to precision treatments through transcriptional state changes, epigenetic adaptation, and metabolic rewiring, thus facilitating the development of refractory disease while simultaneously exposing novel vulnerabilities.

Tracking myeloma tumor DNA in peripheral blood.

Over the past years, the emergence of liquid biopsy technologies has dramatically expanded our ability to assess multiple myeloma without the need for invasive sampling. Interrogation of cell-free DNA from the peripheral blood recapitulates the mutational landscape at excellent concordance with matching bone marrow aspirates. It can quantify disease burden and identify previously undetected resistance mechanisms which may inform clinical management in real-time. The convenience of sample acquisition and storage provides strong procedural benefits over currently available testing. Further investigations will have to define the role of cell-free DNA as a diagnostic measure by determining clinically relevant tumor thresholds in comparison to existing routine parameters. This review presents an overview of currently available assays and discusses the clinical value, potential and limitations of cell-free DNA technologies for the assessment of this challenging disease.

Comprehensive characterization of circulating and bone marrow-derived multiple myeloma cells at minimal residual disease.

The presence or absence of minimal residual disease (MRD) in patients with multiple myeloma (MM) has emerged as a useful marker to determine the depth of remission. MRD negativity as an endpoint has been shown to be associated with improved progression-free survival in many studies. MRD detection is therefore part of numerous clinical trial protocols for MM. At the present time, two methodologies are most widely accepted for MRD detection: (1) multicolor flow cytometry and (2) next-generation sequencing-based clonotype detection. While both of those methodologies enable accurate quantification of MRD in the bone marrow (BM), with sensitivity as low as 10-5 to 10-6, there are several limitations to these methods. First, these approaches reveal the presence or absence of MRD but provide limited molecular information about MM. More comprehensive characterization of MM cells at the MRD stage may identify molecular mechanisms of drug resistance. Second, MRD detection in the BM is typically performed at one time point only, but more frequent detection may define the duration of the MRD status and thus refine its prognostic value. Third, less-invasive approaches that avoid the discomfort and risk associated with BM biopsy would be highly desirable, especially in elderly or frail patients. "Liquid biopsy" for the detection and characterization of circulating MM cells may address these issues. Although MRD detection in the peripheral blood at the same sensitivity as in the BM may be challenging, the identification of patients who do not achieve MRD negativity might reduce the need for BM biopsies. Here, we give an overview of approaches that have been described to detect and characterize MM cells when they occur at very low frequencies in the peripheral blood or in the BM, emphasizing recently described next-generation sequencing approaches for more comprehensive characterization of circulating MM cells.

The Mutational Landscape of Circulating Tumor Cells in Multiple Myeloma.

The development of sensitive and non-invasive "liquid biopsies" presents new opportunities for longitudinal monitoring of tumor dissemination and clonal evolution. The number of circulating tumor cells (CTCs) is prognostic in multiple myeloma (MM), but there is little information on their genetic features. Here, we have analyzed the genomic landscape of CTCs from 29 MM patients, including eight cases with matched/paired bone marrow (BM) tumor cells. Our results show that 100% of clonal mutations in patient BM were detected in CTCs and that 99% of clonal mutations in CTCs were present in BM MM. These include typical driver mutations in MM such as in KRAS, NRAS, or BRAF. These data suggest that BM and CTC samples have similar clonal structures, as discordances between the two were restricted to subclonal mutations. Accordingly, our results pave the way for potentially less invasive mutation screening of MM patients through characterization of CTCs.

The Role of Minimal Residual Disease Testing in Myeloma Treatment Selection and Drug Development: Current Value and Future Applications.

Treatment of myeloma has benefited from the introduction of more effective and better tolerated agents, improvements in supportive care, better understanding of disease biology, revision of diagnostic criteria, and new sensitive and specific tools for disease prognostication and management. Assessment of minimal residual disease (MRD) in response to therapy is one of these tools, as longer progression-free survival (PFS) is seen consistently among patients who have achieved MRD negativity. Current therapies lead to unprecedented frequency and depth of response, and next-generation flow and sequencing methods to measure MRD in bone marrow are in use and being developed with sensitivities in the range of 10-5 to 10-6 cells. These technologies may be combined with functional imaging to detect MRD outside of bone marrow. Moreover, immune profiling methods are being developed to better understand the immune environment in myeloma and response to immunomodulatory agents while methods for molecular profiling of myeloma cells and circulating DNA in blood are also emerging. With the continued development and standardization of these methodologies, MRD has high potential for use in gaining new drug approvals in myeloma. The FDA has outlined two pathways by which MRD could be qualified as a surrogate endpoint for clinical studies directed at obtaining accelerated approval for new myeloma drugs. Most importantly, better understanding of MRD should also contribute to better treatment monitoring. Potentially, MRD status could be used as a prognostic factor for making treatment decisions and for informing timing of therapeutic interventions. Clin Cancer Res; 23(15); 3980-93. ©2017 AACR.

Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors.

Whole-exome sequencing of cell-free DNA (cfDNA) could enable comprehensive profiling of tumors from blood but the genome-wide concordance between cfDNA and tumor biopsies is uncertain. Here we report ichorCNA, software that quantifies tumor content in cfDNA from 0.1× coverage whole-genome sequencing data without prior knowledge of tumor mutations. We apply ichorCNA to 1439 blood samples from 520 patients with metastatic prostate or breast cancers. In the earliest tested sample for each patient, 34% of patients have ≥10% tumor-derived cfDNA, sufficient for standard coverage whole-exome sequencing. Using whole-exome sequencing, we validate the concordance of clonal somatic mutations (88%), copy number alterations (80%), mutational signatures, and neoantigens between cfDNA and matched tumor biopsies from 41 patients with ≥10% cfDNA tumor content. In summary, we provide methods to identify patients eligible for comprehensive cfDNA profiling, revealing its applicability to many patients, and demonstrate high concordance of cfDNA and metastatic tumor whole-exome sequencing.

Genetic interrogation of circulating multiple myeloma cells at single-cell resolution.

Multiple myeloma (MM) remains an incurable disease, with a treatment-refractory state eventually developing in all patients. Constant clonal evolution and genetic heterogeneity of MM are a likely explanation for the emergence of drug-resistant disease. Monitoring of MM genomic evolution on therapy by serial bone marrow biopsy is unfortunately impractical because it involves an invasive and painful procedure. We describe how noninvasive and highly sensitive isolation and characterization of circulating tumor cells (CTCs) from peripheral blood at single-cell resolution recapitulate MM in the bone marrow. We demonstrate that CTCs provide the same genetic information as bone marrow MM cells and even reveal mutations with greater sensitivity than bone marrow biopsies in some cases. Single CTC RNA sequencing enables classification of MM and quantitative assessment of genes that are relevant for prognosis. We propose that the genomic characterization of CTCs should be included in clinical trials to follow the emergence of resistant subclones after MM therapy.

Widespread genetic heterogeneity in multiple myeloma: implications for targeted therapy.

We performed massively parallel sequencing of paired tumor/normal samples from 203 multiple myeloma (MM) patients and identified significantly mutated genes and copy number alterations and discovered putative tumor suppressor genes by determining homozygous deletions and loss of heterozygosity. We observed frequent mutations in KRAS (particularly in previously treated patients), NRAS, BRAF, FAM46C, TP53, and DIS3 (particularly in nonhyperdiploid MM). Mutations were often present in subclonal populations, and multiple mutations within the same pathway (e.g., KRAS, NRAS, and BRAF) were observed in the same patient. In vitro modeling predicts only partial treatment efficacy of targeting subclonal mutations, and even growth promotion of nonmutated subclones in some cases. These results emphasize the importance of heterogeneity analysis for treatment decisions.

An epigenetic mechanism of resistance to targeted therapy in T cell acute lymphoblastic leukemia.

The identification of activating NOTCH1 mutations in T cell acute lymphoblastic leukemia (T-ALL) led to clinical testing of γ-secretase inhibitors (GSIs) that prevent NOTCH1 activation. However, responses to these inhibitors have been transient, suggesting that resistance limits their clinical efficacy. Here we modeled T-ALL resistance, identifying GSI-tolerant 'persister' cells that expand in the absence of NOTCH1 signaling. Rare persisters are already present in naive T-ALL populations, and the reversibility of their phenotype suggests an epigenetic mechanism. Relative to GSI-sensitive cells, persister cells activate distinct signaling and transcriptional programs and exhibit chromatin compaction. A knockdown screen identified chromatin regulators essential for persister viability, including BRD4. BRD4 binds enhancers near critical T-ALL genes, including MYC and BCL2. The BRD4 inhibitor JQ1 downregulates expression of these targets and induces growth arrest and apoptosis in persister cells, at doses well tolerated by GSI-sensitive cells. Consistently, the GSI-JQ1 combination was found to be effective against primary human leukemias in vivo. Our findings establish a role for epigenetic heterogeneity in leukemia resistance that may be addressed by incorporating epigenetic modulators in combination therapy.