Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia.

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

Allele-specific expression of GATA2 due to epigenetic dysregulation in CEBPA double-mutant AML.

Transcriptional deregulation is a central event in the development of acute myeloid leukemia (AML). To identify potential disturbances in gene regulation, we conducted an unbiased screen of allele-specific expression (ASE) in 209 AML cases. The gene encoding GATA binding protein 2 (GATA2) displayed ASE more often than any other myeloid- or cancer-related gene. GATA2 ASE was strongly associated with CEBPA double mutations (DMs), with 95% of cases presenting GATA2 ASE. In CEBPA DM AML with GATA2 mutations, the mutated allele was preferentially expressed. We found that GATA2 ASE was a somatic event lost in complete remission, supporting the notion that it plays a role in CEBPA DM AML. Acquisition of GATA2 ASE involved silencing of 1 allele via promoter methylation and concurrent overactivation of the other allele, thereby preserving expression levels. Notably, promoter methylation was also lost in remission along with GATA2 ASE. In summary, we propose that GATA2 ASE is acquired by epigenetic mechanisms and is a prerequisite for the development of AML with CEBPA DMs. This finding constitutes a novel example of an epigenetic hit cooperating with a genetic hit in the pathogenesis of AML.

BCL6 maintains survival and self-renewal of primary human acute myeloid leukemia cells.

B-cell lymphoma 6 (BCL6) is a transcription repressor and proto-oncogene that plays a crucial role in the innate and adaptive immune system and lymphoid neoplasms. However, its role in myeloid malignancies remains unclear. Here, we explored the role of BCL6 in acute myeloid leukemia (AML). BCL6 was expressed at variable and often high levels in AML cell lines and primary AML samples. AMLs with higher levels of BCL6 were generally sensitive to treatment with BCL6 inhibitors, with the exception of those with monocytic differentiation. Gene expression profiling of AML cells treated with a BCL6 inhibitor revealed induction of BCL6-repressed target genes and transcriptional programs linked to DNA damage checkpoints and downregulation of stem cell genes. Ex vivo treatment of primary AML cells with BCL6 inhibitors induced apoptosis and decreased colony-forming capacity, which correlated with the levels of BCL6 expression. Importantly, inhibition or knockdown of BCL6 in primary AML cells resulted in a significant reduction of leukemia-initiating capacity in mice, suggesting ablation of leukemia repopulating cell functionality. In contrast, BCL6 knockout or inhibition did not suppress the function of normal hematopoietic stem cells. Treatment with cytarabine further induced BCL6 expression, and the levels of BCL6 induction were correlated with resistance to cytarabine. Treatment of AML patient-derived xenografts with BCL6 inhibitor plus cytarabine suggested enhanced antileukemia activity with this combination. Hence, pharmacologic inhibition of BCL6 might provide a novel therapeutic strategy for ablation of leukemia-repopulating cells and increased responsiveness to chemotherapy.

Therapeutic Targeting of RNA Splicing Catalysis through Inhibition of Protein Arginine Methylation.

Cancer-associated mutations in genes encoding RNA splicing factors (SFs) commonly occur in leukemias, as well as in a variety of solid tumors, and confer dependence on wild-type splicing. These observations have led to clinical efforts to directly inhibit the spliceosome in patients with refractory leukemias. Here, we identify that inhibiting symmetric or asymmetric dimethylation of arginine, mediated by PRMT5 and type I protein arginine methyltransferases (PRMTs), respectively, reduces splicing fidelity and results in preferential killing of SF-mutant leukemias over wild-type counterparts. These data identify genetic subsets of cancer most likely to respond to PRMT inhibition, synergistic effects of combined PRMT5 and type I PRMT inhibition, and a mechanistic basis for the therapeutic efficacy of PRMT inhibition in cancer.

Minimal Residual Disease Monitoring of Acute Myeloid Leukemia by Massively Multiplex Digital PCR in Patients with NPM1 Mutations.

The presence of minimal residual disease (MRD) is widely recognized as a powerful predictor of therapeutic outcome in acute myeloid leukemia (AML), but methods of measurement and quantification of MRD in AML are not yet standardized in clinical practice. There is an urgent, unmet need for robust and sensitive assays that can be readily adopted as real-time tools for disease monitoring. NPM1 frameshift mutations are an established MRD marker present in half of patients with cytogenetically normal AML. However, detection is complicated by the existence of hundreds of potential frameshift insertions, clonal heterogeneity, and absence of sequence information when the NPM1 mutation is identified using capillary electrophoresis. Thus, some patients are ineligible for NPM1 MRD monitoring. Furthermore, a subset of patients with NPM1-mutated AML will have false-negative MRD results because of clonal evolution. To simplify and improve MRD testing for NPM1, we present a novel digital PCR technique composed of massively multiplex pools of insertion-specific primers that selectively detect mutated but not wild-type NPM1. By measuring reaction end points using digital PCR technology, the resulting single assay enables sensitive and specific quantification of most NPM1 exon 12 mutations in a manner that is robust to clonal heterogeneity, does not require NPM1 sequence information, and obviates the need for maintenance of hundreds of type-specific assays and associated plasmid standards.

Epigenetic Identity in AML Depends on Disruption of Nonpromoter Regulatory Elements and Is Affected by Antagonistic Effects of Mutations in Epigenetic Modifiers.

We performed cytosine methylation sequencing on genetically diverse patients with acute myeloid leukemia (AML) and found leukemic DNA methylation patterning is primarily driven by nonpromoter regulatory elements and CpG shores. Enhancers displayed stronger differential methylation than promoters, consisting predominantly of hypomethylation. AMLs with dominant hypermethylation featured greater epigenetic disruption of promoters, whereas those with dominant hypomethylation displayed greater disruption of distal and intronic regions. Mutations in IDH and DNMT3A had opposing and mutually exclusive effects on the epigenome. Notably, co-occurrence of both mutations resulted in epigenetic antagonism, with most CpGs affected by either mutation alone no longer affected in double-mutant AMLs. Importantly, this epigenetic antagonism precedes malignant transformation and can be observed in preleukemic LSK cells from Idh2R140Q or Dnmt3aR882H single-mutant and Idh2R140Q/Dnmt3aR882H double-mutant mice. Notably, IDH/DNMT3A double-mutant AMLs manifested upregulation of a RAS signaling signature and displayed unique sensitivity to MEK inhibition ex vivo as compared with AMLs with either single mutation.Significance: AML is biologically heterogeneous with subtypes characterized by specific genetic and epigenetic abnormalities. Comprehensive DNA methylation profiling revealed that differential methylation of nonpromoter regulatory elements is a driver of epigenetic identity, that gene mutations can be context-dependent, and that co-occurrence of mutations in epigenetic modifiers can result in epigenetic antagonism. Cancer Discov; 7(8); 868-83. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 783.

Epigenetics and approaches to targeted epigenetic therapy in acute myeloid leukemia.

Acute myeloid leukemia (AML) is the most common type of acute leukemia in adults. AML is a heterogeneous malignancy characterized by distinct genetic abnormalities. Recent discoveries have highlighted an additional important role of dysregulated epigenetic mechanisms in the pathogenesis of the disease. In contrast to genetic changes, epigenetic modifications are frequently reversible, which provides opportunities for targeted treatment using specific inhibitors. In this review, we will provide an overview of the current state of epigenetics and epigenetic therapy in AML and will describe perspectives on how to identify promising new approaches for epigenetic targeted treatment.

Risk stratification of intermediate-risk acute myeloid leukemia: integrative analysis of a multitude of gene mutation and gene expression markers.

Numerous molecular markers have been recently discovered as potential prognostic factors in acute myeloid leukemia (AML). It has become of critical importance to thoroughly evaluate their interrelationships and relative prognostic importance. Gene expression profiling was conducted in a well-characterized cohort of 439 AML patients (age < 60 years) to determine expression levels of EVI1, WT1, BCL2, ABCB1, BAALC, FLT3, CD34, INDO, ERG and MN1. A variety of AML-specific mutations were evaluated, that is, FLT3, NPM1, N-RAS, K-RAS, IDH1, IDH2, and CEBPA(DM/SM) (double/single). Univariable survival analysis shows that (1) patients with FLT3(ITD) mutations have inferior overall survival (OS) and event-free survival (EFS), whereas CEBPA(DM) and NPM1 mutations indicate favorable OS and EFS in intermediate-risk AML, and (2) high transcript levels of BAALC, CD34, MN1, EVl1, and ERG predict inferior OS and EFS. In multivariable survival analysis, CD34, ERG, and CEBPA(DM) remain significant. Using survival tree and regression methodologies, we show that CEBPA(DM), CD34, and IDH2 mutations are capable of separating the intermediate group into 2 AML subgroups with highly distinctive survival characteristics (OS at 60 months: 51.9% vs 14.9%). The integrated statistical approach demonstrates that from the multitude of biomarkers a greatly condensed subset can be selected for improved stratification of intermediate-risk AML.

Prognostic impact, concurrent genetic mutations, and gene expression features of AML with CEBPA mutations in a cohort of 1182 cytogenetically normal AML patients: further evidence for CEBPA double mutant AML as a distinctive disease entity.

We evaluated concurrent gene mutations, clinical outcome, and gene expression signatures of CCAAT/enhancer binding protein alpha (CEBPA) double mutations (CEBPA(dm)) versus single mutations (CEBPA(sm)) in 1182 cytogenetically normal acute myeloid leukemia (AML) patients (16-60 years of age). We identified 151 (12.8%) patients with CEBPA mutations (91 CEBPA(dm) and 60 CEBPA(sm)). The incidence of germline mutations was 7% (5 of 71), including 3 C-terminal mutations. CEBPA(dm) patients had a lower frequency of concurrent mutations than CEBPA(sm) patients (P < .0001). Both, groups were associated with a favorable outcome compared with CEBPA(wt) (5-year overall survival [OS] 63% and 56% vs 39%; P < .0001 and P = .05, respectively). However, in multivariable analysis only CEBPA(dm) was a prognostic factor for favorable OS outcome (hazard ratio [HR] 0.36, P < .0001; event-free survival, HR 0.41, P < .0001; relapse-free survival, HR 0.55, P = .001). Outcome in CEBPA(sm) is dominated by concurrent NPM1 and/or FLT3 internal tandem duplication mutations. Unsupervised and supervised GEP analyses showed that CEBPA(dm) AML (n = 42), but not CEBPA(sm) AML (n = 18), expressed a unique gene signature. A 25-probe set prediction signature for CEBPA(dm) AML showed 100% sensitivity and specificity. Based on these findings, we propose that CEBPA(dm) should be clearly defined from CEBPA(sm) AML and considered as a separate entity in the classification of AML.

HAT: hypergeometric analysis of tiling-arrays with application to promoter-GeneChip data.

BACKGROUND: Tiling-arrays are applicable to multiple types of biological research questions. Due to its advantages (high sensitivity, resolution, unbiased), the technology is often employed in genome-wide investigations. A major challenge in the analysis of tiling-array data is to define regions-of-interest, i.e., contiguous probes with increased signal intensity (as a result of hybridization of labeled DNA) in a region. Currently, no standard criteria are available to define these regions-of-interest as there is no single probe intensity cut-off level, different regions-of-interest can contain various numbers of probes, and can vary in genomic width. Furthermore, the chromosomal distance between neighboring probes can vary across the genome among different arrays. RESULTS: We have developed Hypergeometric Analysis of Tiling-arrays (HAT), and first evaluated its performance for tiling-array datasets from a Chromatin Immunoprecipitation study on chip (ChIP-on-chip) for the identification of genome-wide DNA binding profiles of transcription factor Cebpa (used for method comparison). Using this assay, we can refine the detection of regions-of-interest by illustrating that regions detected by HAT are more highly enriched for expected motifs in comparison with an alternative detection method (MAT). Subsequently, data from a retroviral insertional mutagenesis screen were used to examine the performance of HAT among different applications of tiling-array datasets. In both studies, detected regions-of-interest have been validated with (q)PCR. CONCLUSIONS: We demonstrate that HAT has increased specificity for analysis of tiling-array data in comparison with the alternative method, and that it accurately detects regions-of-interest in two different applications of tiling-arrays. HAT has several advantages over previous methods: i) as there is no single cut-off level for probe-intensity, HAT can detect regions-of-interest at various thresholds, ii) it can detect regions-of-interest of any size, iii) it is independent of probe-resolution across the genome, and across tiling-array platforms and iv) it employs a single user defined parameter: the significance level. Regions-of-interest are detected by computing the hypergeometric-probability, while controlling the Family Wise Error. Furthermore, the method does not require experimental replicates, common regions-of-interest are indicated, a sequence-of-interest can be examined for every detected region-of-interest, and flanking genes can be reported.

Gene expression profiling of minimally differentiated acute myeloid leukemia: M0 is a distinct entity subdivided by RUNX1 mutation status.

Minimally differentiated acute myeloid leukemia (AML-M0) is defined by immature morphology and expression of early hematologic markers. By gene expression profiling (GEP) and subsequent unsupervised analysis of 35 AML-M0 samples and 253 previously reported AML cases, we demonstrate that AML-M0 cases express a unique signature that is largely separated from other molecular subtypes. Hematologic transcription regulators such as CEBPA, CEBPD, and ETV6, and the differentiation associated gene MPO appeared strongly down-regulated, in line with the primitive state of this leukemia. AML-M0 frequently carries loss-of-function RUNX1 mutation. Unsupervised analyses revealed a subdivision between AML-M0 cases with and without RUNX1 mutations. RUNX1 mutant AML-M0 samples showed a distinct up-regulation of B cell-related genes such as members of the B-cell receptor complex, transcription regulators RUNX3, ETS2, IRF8, or PRDM1, and major histocompatibility complex class II genes. Importantly, prediction with high accuracy of the AML-M0 subtype and prediction of patients carrying RUNX1 mutation within this subtype were possible based on the expression level of only a few transcripts. We propose that RUNX1 mutations in this AML subgroup cause lineage infidelity, leading to aberrant coexpression of myeloid and B-lymphoid genes. Furthermore, our results imply that AML-M0, although originally determined by morphology, constitutes a leukemia subgroup.

Genome-wide epigenetic analysis delineates a biologically distinct immature acute leukemia with myeloid/T-lymphoid features.

Acute myeloid leukemia is a heterogeneous disease from the molecular and biologic standpoints, and even patients with a specific gene expression profile may present clinical and molecular heterogeneity. We studied the epigenetic profiles of a cohort of patients who shared a common gene expression profile but differed in that only half of them harbored mutations of the CEBPA locus, whereas the rest presented with silencing of this gene and coexpression of certain T-cell markers. DNA methylation studies revealed that these 2 groups of patients could be readily segregated in an unsupervised fashion based on their DNA methylation profiles alone. Furthermore, CEBPA silencing was associated with the presence of an aberrant DNA hypermethylation signature, which was not present in the CEBPA mutant group. This aberrant hypermethylation occurred more frequently at sites within CpG islands. CEBPA-silenced leukemias also displayed marked hypermethylation compared with normal CD34(+) hematopoietic cells, whereas CEBPA mutant cases showed only mild changes in DNA methylation compared with these normal progenitors. Biologically, CEBPA-silenced leukemias presented with a decreased response to myeloid growth factors in vitro.

Double CEBPA mutations, but not single CEBPA mutations, define a subgroup of acute myeloid leukemia with a distinctive gene expression profile that is uniquely associated with a favorable outcome.

Mutations in CCAAT/enhancer binding protein alpha (CEBPA) are seen in 5% to 14% of acute myeloid leukemia (AML) and have been associated with a favorable clinical outcome. Most AMLs with CEBPA mutations simultaneously carry 2 mutations (CEBPA(double-mut)), usually biallelic, whereas single heterozygous mutations (CEBPA(single-mut)) are less frequently seen. Using denaturing high-performance liquid chromatography and nucleotide sequencing, we identified among a cohort of 598 newly diagnosed AMLs a subset of 41 CEBPA mutant cases (28 CEBPA(double-mut) and 13 CEBPA(single-mut) cases). CEBPA(double-mut) associated with a unique gene expression profile as well as favorable overall and event-free survival, retained in multivariable analysis that included cytogenetic risk, FLT3-ITD and NPM1 mutation, white blood cell count, and age. In contrast, CEBPA(single-mut) AMLs did not express a discriminating signature and could not be distinguished from wild-type cases as regards clinical outcome. These results demonstrate significant underlying heterogeneity within CEBPA mutation-positive AML with prognostic relevance.

Prediction of molecular subtypes in acute myeloid leukemia based on gene expression profiling.

We examined the gene expression profiles of two independent cohorts of patients with acute myeloid leukemia [n=247 and n=214 (younger than or equal to 60 years)] to study the applicability of gene expression profiling as a single assay in prediction of acute myeloid leukemia-specific molecular subtypes. The favorable cytogenetic acute myeloid leukemia subtypes, i.e., acute myeloid leukemia with t(8;21), t(15;17) or inv(16), were predicted with maximum accuracy (positive and negative predictive value: 100%). Mutations in NPM1 and CEBPA were predicted less accurately (positive predictive value: 66% and 100%, and negative predictive value: 99% and 97% respectively). Various other characteristic molecular acute myeloid leukemia subtypes, i.e., mutant FLT3 and RAS, abnormalities involving 11q23, -5/5q-, -7/7q-, abnormalities involving 3q (abn3q) and t(9;22), could not be correctly predicted using gene expression profiling. In conclusion, gene expression profiling allows accurate prediction of certain acute myeloid leukemia subtypes, e.g. those characterized by expression of chimeric transcription factors. However, detection of mutations affecting signaling molecules and numerical abnormalities still requires alternative molecular methods.

A decade of genome-wide gene expression profiling in acute myeloid leukemia: flashback and prospects.

The past decade has shown a marked increase in the use of high-throughput assays in clinical research into human cancer, including acute myeloid leukemia (AML). In particular, genome-wide gene expression profiling (GEP) using DNA microarrays has been extensively used for improved understanding of the diagnosis, prognosis, and pathobiology of this heterogeneous disease. This review discusses the progress that has been made, places the technologic limitations in perspective, and highlights promising future avenues.

Gene expression profiling for improved dissection of acute leukemia: a recently identified immature myeloid/T-lymphoid subgroup as an example.

In this concise overview, we discuss recent findings concerning a distinct subgroup of acute myeloid/T-lymphoid leukemia. We describe how we identified these leukemias in multiple cohorts of acute myeloid leukemia (AML) using a combination of gene expression profiling and additional analytic approaches, and how we obtained insight in possible mechanisms leading to their phenotype.

Distinct gene expression profiles of acute myeloid/T-lymphoid leukemia with silenced CEBPA and mutations in NOTCH1.

Gene expression profiling of acute myeloid leukemia (AML) allows the discovery of previously unrecognized molecular entities. Here, we identified a specific subgroup of AML, defined by an expression profile resembling that of AMLs with mutations in the myeloid transcription factor CCAAT/enhancer-binding protein alpha (C/EBPalpha), while lacking such mutations. We found that in these leukemias, the CEBPA gene was silenced, which was associated with frequent promoter hypermethylation. The leukemias phenotypically showed aberrant expression of T-cell genes, of which CD7 was most consistent. We identified 2 mechanisms that may contribute to this phenotype. First, absence of Cebpa led to up-regulation of specific T-cell transcripts (ie, Cd7 and Lck) in hematopoietic stem cells isolated from conditional Cebpa knockout mice. Second, the enhanced expression of TRIB2, which we identify here as a direct target of the T-cell commitment factor NOTCH1, suggested aberrantly activated Notch signaling. Putatively activating NOTCH1 mutations were found in several specimens of the newly identified subgroup, while a large set of control AMLs was mutation negative. A gene expression prediction signature allowed the detection of similar cases of leukemia in independent series of AML.

Text-derived concept profiles support assessment of DNA microarray data for acute myeloid leukemia and for androgen receptor stimulation.

BACKGROUND: High-throughput experiments, such as with DNA microarrays, typically result in hundreds of genes potentially relevant to the process under study, rendering the interpretation of these experiments problematic. Here, we propose and evaluate an approach to find functional associations between large numbers of genes and other biomedical concepts from free-text literature. For each gene, a profile of related concepts is constructed that summarizes the context in which the gene is mentioned in literature. We assign a weight to each concept in the profile based on a likelihood ratio measure. Gene concept profiles can then be clustered to find related genes and other concepts. RESULTS: The experimental validation was done in two steps. We first applied our method on a controlled test set. After this proved to be successful the datasets from two DNA microarray experiments were analyzed in the same way and the results were evaluated by domain experts. The first dataset was a gene-expression profile that characterizes the cancer cells of a group of acute myeloid leukemia patients. For this group of patients the biological background of the cancer cells is largely unknown. Using our methodology we found an association of these cells to monocytes, which agreed with other experimental evidence. The second data set consisted of differentially expressed genes following androgen receptor stimulation in a prostate cancer cell line. Based on the analysis we put forward a hypothesis about the biological processes induced in these studied cells: secretory lysosomes are involved in the production of prostatic fluid and their development and/or secretion are androgen-regulated processes. CONCLUSION: Our method can be used to analyze DNA microarray datasets based on information explicitly and implicitly available in the literature. We provide a publicly available tool, dubbed Anni, for this purpose.