Dynamic changes in the clonal structure of MDS and AML in response to epigenetic therapy.

Traditional response criteria in myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML) are based on bone marrow morphology and may not accurately reflect clonal tumor burden in patients treated with non-cytotoxic chemotherapy. We used next-generation sequencing of serial bone marrow samples to monitor MDS and AML tumor burden during treatment with epigenetic therapy (decitabine and panobinostat). Serial bone marrow samples (and skin as a source of normal DNA) from 25 MDS and AML patients were sequenced (exome or 285 gene panel). We observed that responders, including those in complete remission (CR), can have persistent measurable tumor burden (that is, mutations) for at least 1 year without disease progression. Using an ultrasensitive sequencing approach, we detected extremely rare mutations (equivalent to 1 heterozygous mutant cell in 2000 non-mutant cells) months to years before their expansion at disease relapse. While patients can live with persistent clonal hematopoiesis in a CR or stable disease, ultimately we find evidence that expansion of a rare subclone occurs at relapse or progression. Here we demonstrate that sequencing of serial samples provides an alternative measure of tumor burden in MDS or AML patients and augments traditional response criteria that rely on bone marrow blast percentage.

Epigenomic analysis of the HOX gene loci reveals mechanisms that may control canonical expression patterns in AML and normal hematopoietic cells.

HOX genes are highly expressed in many acute myeloid leukemia (AML) samples, but the patterns of expression and associated regulatory mechanisms are not clearly understood. We analyzed RNA sequencing data from 179 primary AML samples and normal hematopoietic cells to understand the range of expression patterns in normal versus leukemic cells. HOX expression in AML was restricted to specific genes in the HOXA or HOXB loci, and was highly correlated with recurrent cytogenetic abnormalities. However, the majority of samples expressed a canonical set of HOXA and HOXB genes that was nearly identical to the expression signature of normal hematopoietic stem/progenitor cells. Transcriptional profiles at the HOX loci were similar between normal cells and AML samples, and involved bidirectional transcription at the center of each gene cluster. Epigenetic analysis of a subset of AML samples also identified common regions of chromatin accessibility in AML samples and normal CD34(+) cells that displayed differences in methylation depending on HOX expression patterns. These data provide an integrated epigenetic view of the HOX gene loci in primary AML samples, and suggest that HOX expression in most AML samples represents a normal stem cell program that is controlled by epigenetic mechanisms at specific regulatory elements.

U2AF1 mutations alter sequence specificity of pre-mRNA binding and splicing.

We previously identified missense mutations in the U2AF1 splicing factor affecting codons S34 (S34F and S34Y) or Q157 (Q157R and Q157P) in 11% of the patients with de novo myelodysplastic syndrome (MDS). Although the role of U2AF1 as an accessory factor in the U2 snRNP is well established, it is not yet clear how these mutations affect splicing or contribute to MDS pathophysiology. We analyzed splice junctions in RNA-seq data generated from transfected CD34+ hematopoietic cells and found significant differences in the abundance of known and novel junctions in samples expressing mutant U2AF1 (S34F). For selected transcripts, splicing alterations detected by RNA-seq were confirmed by analysis of primary de novo MDS patient samples. These effects were not due to impaired U2AF1 (S34F) localization as it co-localized normally with U2AF2 within nuclear speckles. We further found evidence in the RNA-seq data for decreased affinity of U2AF1 (S34F) for uridine (relative to cytidine) at the e-3 position immediately upstream of the splice acceptor site and corroborated this finding using affinity-binding assays. These data suggest that the S34F mutation alters U2AF1 function in the context of specific RNA sequences, leading to aberrant alternative splicing of target genes, some of which may be relevant for MDS pathogenesis.

Clonal evolution revealed by whole genome sequencing in a case of primary myelofibrosis transformed to secondary acute myeloid leukemia.

Clonal architecture in myeloproliferative neoplasms (MPNs) is poorly understood. Here we report genomic analyses of a patient with primary myelofibrosis (PMF) transformed to secondary acute myeloid leukemia (sAML). Whole genome sequencing (WGS) was performed on PMF and sAML diagnosis samples, with skin included as a germline surrogate. Deep sequencing validation was performed on the WGS samples and an additional sample obtained during sAML remission/relapsed PMF. Clustering analysis of 649 validated somatic single-nucleotide variants revealed four distinct clonal groups, each including putative driver mutations. The first group (including JAK2 and U2AF1), representing the founding clone, included mutations with high frequency at all three disease stages. The second clonal group (including MYB) was present only in PMF, suggesting the presence of a clone that was dispensable for transformation. The third group (including ASXL1) contained mutations with low frequency in PMF and high frequency in subsequent samples, indicating evolution of the dominant clone with disease progression. The fourth clonal group (including IDH1 and RUNX1) was acquired at sAML transformation and was predominantly absent at sAML remission/relapsed PMF. Taken together, these findings illustrate the complex clonal dynamics associated with disease evolution in MPNs and sAML.

Clonal diversity of recurrently mutated genes in myelodysplastic syndromes.

Recent studies suggest that most cases of myelodysplastic syndrome (MDS) are clonally heterogeneous, with a founding clone and multiple subclones. It is not known whether specific gene mutations typically occur in founding clones or subclones. We screened a panel of 94 candidate genes in a cohort of 157 patients with MDS or secondary acute myeloid leukemia (sAML). This included 150 cases with samples obtained at MDS diagnosis and 15 cases with samples obtained at sAML transformation (8 were also analyzed at the MDS stage). We performed whole-genome sequencing (WGS) to define the clonal architecture in eight sAML genomes and identified the range of variant allele frequencies (VAFs) for founding clone mutations. At least one mutation or cytogenetic abnormality was detected in 83% of the 150 MDS patients and 17 genes were significantly mutated (false discovery rate ≤0.05). Individual genes and patient samples displayed a wide range of VAFs for recurrently mutated genes, indicating that no single gene is exclusively mutated in the founding clone. The VAFs of recurrently mutated genes did not fully recapitulate the clonal architecture defined by WGS, suggesting that comprehensive sequencing may be required to accurately assess the clonal status of recurrently mutated genes in MDS.

Recurrent DNMT3A mutations in patients with myelodysplastic syndromes.

Alterations in DNA methylation have been implicated in the pathogenesis of myelodysplastic syndromes (MDS), although the underlying mechanism remains largely unknown. Methylation of CpG dinucleotides is mediated by DNA methyltransferases, including DNMT1, DNMT3A and DNMT3B. DNMT3A mutations have recently been reported in patients with de novo acute myeloid leukemia (AML), providing a rationale for examining the status of DNMT3A in MDS samples. In this study, we report the frequency of DNMT3A mutations in patients with de novo MDS, and their association with secondary AML. We sequenced all coding exons of DNMT3A using DNA from bone marrow and paired normal cells from 150 patients with MDS and identified 13 heterozygous mutations with predicted translational consequences in 12/150 patients (8.0%). Amino acid R882, located in the methyltransferase domain of DNMT3A, was the most common mutation site, accounting for 4/13 mutations. DNMT3A mutations were expressed in the majority of cells in all tested mutant samples regardless of myeloblast counts, suggesting that DNMT3A mutations occur early in the course of MDS. Patients with DNMT3A mutations had worse overall survival compared with patients without DNMT3A mutations (P=0.005) and more rapid progression to AML (P=0.007), suggesting that DNMT3A mutation status may have prognostic value in de novo MDS.

POU4F1 is associated with t(8;21) acute myeloid leukemia and contributes directly to its unique transcriptional signature.

The t(8;21)(q22;q22) translocation, present in approximately 5% of adult acute myeloid leukemia (AML) cases, produces the AML1/ETO (AE) fusion protein. Dysregulation of the Pit/Oct/Unc (POU) domain-containing transcription factor POU4F1 is a recurring abnormality in t(8;21) AML. In this study, we showed that POU4F1 overexpression is highly correlated with, but not caused by, AE. We observed that AE markedly increases the self-renewal capacity of myeloid progenitors from murine bone marrow or fetal liver and drives the expansion of these cells in liquid culture. POU4F1 is neither necessary nor sufficient for these AE-dependent properties, suggesting that it contributes to leukemia through novel mechanisms. To identify targets of POU4F1, we performed gene expression profiling in primary mouse cells with genetically defined levels of POU4F1 and identified 140 differentially expressed genes. This expression signature was significantly enriched in human t(8;21) AML samples and was sufficient to cluster t(8;21) AML samples in an unsupervised hierarchical analysis. Among the most highly differentially expressed genes, half are known AML1/ETO targets, implying that the unique transcriptional signature of t(8;21) AML is, in part, attributable to POU4F1 and not AML1/ETO itself. These genes provide novel candidates for understanding the biology and developing therapeutic approaches for t(8;21) AML.