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.

Granzyme B can cause mitochondrial depolarization and cell death in the absence of BID, BAX, and BAK.

Granzyme B (GzmB) is a serine protease that is used by activated cytotoxic T lymphocytes to induce target cell apoptosis. Although GzmB directly cleaves the Bcl2 family member BID on target cell entry, Bid-deficient (and Bax, Bak doubly deficient) cells are susceptible to GzmB-induced death, even though they fail to release cytochrome c from mitochondria. GzmB still induces mitochondrial depolarization in Bax, Bak double knockout cells without cytochrome c release or opening of the permeability transition pore. Because GzmB cannot directly cause depolarization of isolated mitochondria, novel intracellular factor(s) may be required for GzmB to depolarize mitochondria in situ. GzmB therefore utilizes two distinct mitochondrial pathways to amplify the proapoptotic signal that it delivers to target cells.

Mouse models of acute promyelocytic leukemia.

Translocations involving a variety of fusion partners, such as promyelocytic leukemia gene, promyelocytic leukemia zinc finger, nucleophosmin, nuclear matrix protein, and signal transducer and activator of transcription protein 5B, with the retinoic acid receptor alpha gene are commonly associated with development of acute promyelocytic leukemia. Through the development of transgenic mouse models, some retinoic acid receptor alpha translocation fusion proteins have been shown to be capable of initiating acute promyelocytic leukemia development, and dictate the leukemias' responsiveness to retinoic acid. Transgenic mouse models also have identified the influence of reciprocal translocation fusion proteins on acute promyelocytic leukemia development, and have demonstrated that additional mutations can contribute to the development of acute promyelocytic leukemia. In this review, the authors summarize current mouse models of acute promyelocytic leukemia and describe current knowledge about additional genetic alterations that occur during development of acute promyelocytic leukemia in the mouse.

Sudden death among patients with acute promyelocytic leukemia treated with arsenic trioxide.

Arsenic trioxide has been shown to be effective in treating acute promyelocytic leukemia (APL), with minimal overall toxicity reported to date. A phase I/II study was initiated in June 1998 using arsenic trioxide for relapsed APL to determine the maximum tolerated or minimal effective dose and to determine the efficacy of treatment at that dose. Ten patients received 1 to 4 monthly cycles of treatment with 0.1 mg/kg per day intravenous arsenic trioxide. Six of 7 patients evaluable for response achieved cytogenetic or molecular complete remission. However, 3 patients died suddenly during the first cycle of treatment. Autopsies obtained on 2 of these failed to identify a cause of sudden death, despite evidence of pulmonary hemorrhage in one. A third patient, for whom an autopsy was not performed, became asystolic and died while on continuous cardiac telemetry. These observations suggest that arsenic trioxide may be significantly or even fatally toxic at doses currently used and that caution is warranted in its use.

Lack of neighborhood effects from a transcriptionally active phosphoglycerate kinase-neo cassette located between the murine beta-major and beta-minor globin genes.

For the treatment of beta-globin gene defects, a homologous recombination-mediated gene correction approach would provide advantages over random integration-based gene therapy strategies. However, "neighborhood effects" from retained selectable marker genes in the targeted locus are among the key issues that must be taken into consideration for any attempt to use this strategy for gene correction. An Ala-to-Ile mutation was created in the beta6 position of the mouse beta-major globin gene (beta(6I)) as a step toward the development of a murine model system that could serve as a platform for therapeutic gene correction studies. The marked beta-major gene can be tracked at the level of DNA, RNA, and protein, allowing investigation of the impact of a retained phosphoglycerate kinase (PGK)-neo cassette located between the mutant beta-major and beta-minor globin genes on expression of these 2 neighboring genes. Although the PGK-neo cassette was expressed at high levels in adult erythroid cells, the abundance of the beta(6I) mRNA was indistinguishable from that of the wild-type counterpart in bone marrow cells. Similarly, the output from the beta-minor globin gene was also normal. Therefore, in this specific location, the retained, transcriptionally active PGK-neo cassette does not disrupt the regulated expression of the adult beta-globin genes. (Blood. 2001;98:65-73)

Dipeptidyl peptidase I is essential for activation of mast cell chymases, but not tryptases, in mice.

Dipeptidyl peptidase I (DPPI) is the sole activator in vivo of several granule-associated serine proteases of cytotoxic lymphocytes. In vitro, DPPI also activates mast cell chymases and tryptases. To determine whether DPPI is essential for their activation in vivo, we used enzyme histochemical and immunohistochemical approaches and solution-based activity assays to study these enzymes in tissues and bone marrow-derived mast cells (BMMCs) from DPPI +/+ and DPPI -/- mice. We find that DPPI -/- mast cells contain normal amounts of immunoreactive chymases but no chymase activity, indicating that DPPI is essential for chymase activation and suggesting that DPPI -/- mice are functional chymase knockouts. The absence of DPPI and chymase activity does not affect the growth, granularity, or staining characteristics of BMMCs and, despite prior predictions, does not alter IgE-mediated exocytosis of histamine. In contrast, the level of active tryptase (mMCP-6) in DPPI -/- BMMCs is 25% that of DPPI +/- BMMCs. These findings indicate that DPPI is not essential for mMCP-6 activation but does influence the total amount of active mMCP-6 in mast cells and therefore may be an important, but not exclusive mechanism for tryptase activation.

Acquired, nonrandom chromosomal abnormalities associated with the development of acute promyelocytic leukemia in transgenic mice.

We previously generated a transgenic mouse model for acute promyelocytic leukemia (APL) by expressing the promyelocytic leukemia (PML)-retinoic acid receptor (RARalpha) cDNA in early myeloid cells. This fusion protein causes a myeloproliferative disease in 100% of animals, but only 15-20% of the animals develop acute leukemia after a long latency period (6-13 months). PML-RARalpha is therefore necessary, but not sufficient, for APL development. The coexpression of a reciprocal form of the fusion, RARalpha-PML, increased the likelihood of APL development (55-60%), but did not shorten latency. Together, these results suggested that additional genetic events are required for the development of APL. We therefore evaluated the splenic tumor cells from 18 transgenic mice with APL for evidence of secondary genetic events, by using spectral karyotyping analysis. Interstitial or terminal deletions of the distal region of one copy of chromosome 2 [del(2)] were found in 1/5 tumors expressing PML-RARalpha, but in 11/13 tumors expressing both PML-RARalpha and RARalpha-PML (P < 0.05). Leukemic cells that contained a deletion on chromosome 2 often contained additional chromosomal gains (especially of 15), chromosomal losses (especially of 11 or X/Y), or were tetraploid (P

DFF45/ICAD can be directly processed by granzyme B during the induction of apoptosis.

Granzyme B (GzmB) is a component of cytotoxic lymphocyte granules that can rapidly initiate apoptosis in target cells. While several procaspases are cleaved and activated by GzmB, the absolute requirement of caspase activation for GzmB-induced apoptosis is controversial. In this report, we demonstrate that GzmB can initiate apoptosis in the absence of caspase-3 activity by directly cleaving DFF45/ICAD to liberate activated DFF40/CAD. DFF45/ICAD cleavage occurs less efficiently in cells that lack caspase-3 activity, suggesting that the caspases normally amplify the GzmB death signal. DFF45/ICAD-deficient mouse embryo fibroblasts are partially resistant to GzmB-induced death, demonstrating the biological importance of DFF45/ICAD for GzmB-mediated apoptosis.

Independent formation of DnaseI hypersensitive sites in the murine beta-globin locus control region.

Mammalian beta-globin loci are composed of multiple orthologous genes whose expression is erythroid specific and developmentally regulated. The expression of these genes both from the endogenous locus and from transgenes is strongly influenced by a linked 15-kilobase region of clustered DNaseI hypersensitive sites (HSs) known as the locus control region (LCR). The LCR encompasses 5 major HSs, each of which is highly homologous among humans, mice, and other mammals. To analyze the function of individual HSs in the endogenous murine beta-globin LCR, we have used homologous recombination in embryonic stem cells to produce 5 mouse lines, each of which is deficient for 1 of these major HSs. In this report, we demonstrate that deletion of the conserved region of 5'HS 1, 2, 3, 4, or 5/6 abolishes HS formation at the deletion site but has no influence on the formation of the remaining HSs in the LCR. Therefore, in the endogenous murine locus, there is no dominant or initiating site whose formation must precede the formation of the other HSs. This is consistent with the idea that HSs form autonomously. We discuss the implications of these findings for current models of beta-globin regulation.

Deletions within the mouse beta-globin locus control region preferentially reduce beta(min) globin gene expression.

The mouse beta-globin gene cluster is regulated, at least in part, by a locus control region (LCR) composed of several developmentally stable DNase I hypersensitive sites located upstream of the genes. In this report, we examine the level of expression of the beta(min) and beta(maj) genes in adult mice in which HS2, HS3, or HS5,6 has been either deleted or replaced by a selectable marker via homologous recombination in ES cells. Primer extension analysis of RNA extracted from circulating reticulocytes and HPLC analysis of globin chains from peripheral red blood cells revealed that all mutations that reduce the overall output of the locus preferentially decrease beta(min) expression over beta(maj). The implications of these findings for the mechanism by which the LCR controls expression of the beta(maj) and beta(min) promoters are discussed.