Dnmt3a-null hematopoietic stem cells (HSCs) cannot sustain long-term hematopoiesis.
Cooperating c-Kit mutations drive leukemic transformation of Dnmt3a-null HSCs.
Genome sequencing studies of patient samples have implicated the involvement of various components of the epigenetic machinery in myeloid diseases, including the de novo DNA methyltransferase DNMT3A. We have recently shown that Dnmt3a is essential for hematopoietic stem cell differentiation. Here, we investigated the effect of loss of Dnmt3a on hematopoietic transformation by forcing the normally quiescent hematopoietic stem cells to divide in vivo. Mice transplanted with Dnmt3a-null bone marrow in the absence of wild-type support cells succumbed to bone marrow failure (median survival, 328 days) characteristic of myelodysplastic syndromes with symptoms including anemia, neutropenia, bone marrow hypercellularity, and splenomegaly with myeloid infiltration. Two out of 25 mice developed myeloid leukemia with >20% blasts in the blood and bone marrow. Four out of 25 primary mice succumbed to myeloproliferative disorders, some of which progressed to secondary leukemia after long latency. Exome sequencing identified cooperating c-Kit mutations found only in the leukemic samples. Ectopic introduction of c-Kit variants into a Dnmt3a-deficient background produced acute leukemia with a short latency (median survival, 67 days). Our data highlight crucial roles of Dnmt3a in normal and malignant hematopoiesis and suggest that a major role for this enzyme is to facilitate developmental progression of progenitor cells at multiple decision checkpoints.
Hematopoietic stem cell (HSC) fate decisions are controlled by signaling pathways, cues from the niche, and the actions of cell-autonomous regulators such as transcription factors, but they are now also recognized to be influenced by a significant epigenetic component. DNA methylation is one of the major epigenetic modifications in the vertebrate genome and is important for development, stem cell differentiation, and oncogenesis.1-3 DNA methylation is catalyzed by the DNA methyltransferase enzymes Dnmt1, Dnmt3a, and Dnmt3b.4-6 Genome sequencing studies of myeloid malignancies have identified recurrent DNMT3A somatic mutations in approximately 22%,7,8 10%,9,10 and 8%11,12 of patients with acute myeloid leukemia (AML), myelodysplastic syndromes (MDS), and myeloproliferative neoplasms (MPN), respectively, and are associated with poor prognosis.13 The most prevalent DNMT3A mutation is an R882H variant that produces a protein that acts as a dominant negative.14,15 DNMT3A nonsense and frameshift mutations are all predicted to result in truncated proteins that eliminate or shorten the methyltransferase domain or are associated with nonsense-mediated decay, suggesting loss of function.8
We had previously studied the role of Dnmt3a in HSC function.16,17 Inducible conditional knockout mice were generated by crossing Dnmt3afl/fl mice18 with the Mx1-CRE driver, with deletion of floxed Dnmt3a alleles induced by sequential injections of polyinosinic-polycytidylic acid (pIpC; henceforth referred to as Dnmt3a-KO). Serial transfer of Dnmt3a-KO HSCs with fresh wild-type (WT) whole bone marrow (WBM) competitor revealed a decline in the differentiated cell output from Dnmt3a-KO over multiple rounds of transplantation on a per-HSC basis.16 Although the striking accumulation of phenotypically defined HSCs in the bone marrow of Dnmt3a-KO HSC recipient mice was reminiscent of a myeloid abnormality, these mice did not develop overt disease even when aged to over 14 months posttransplant.16
In these competitive transplants, the presence of WT WBM may have suppressed malignant transformation of the mutant HSCs. Dnmt3a-null HSCs were less proliferative than counterpart control HSCs,16 suggesting that the cellular turnover threshold necessary to generate additional genetic and/or epigenetic lesions required for leukemogenesis was not achieved. To further understand the contribution of Dnmt3a loss of function in hematopoiesis, we performed noncompetitive transplantation of Dnmt3a-KO bone marrow. This forces the mutant HSCs to divide in vivo to regenerate the hematopoietic system following lethal irradiation and should uncover any predispositions to transformation.
Mice and transplantation
Animal procedures were approved by the Institutional Animal Care and Use Committee and conducted in accordance with Washington University institutional guidelines. All mice were C57Bl/6 background, distinguished by CD45.1 or CD45.2 alleles. Dnmt3afl/fl mice were originally obtained from the Beaudet laboratory at Baylor College of Medicine (with the consent of En Li) and crossed to Mx1-CRE mice. Deletion of floxed alleles was mediated by 6 intraperitoneal injections (300 μg per mouse) of pIpC (Sigma) in phosphate-buffered saline every other day. For primary noncompetitive transplantation, WBM was harvested from donor mice 4 weeks after the last pIpC injection. Recipient mice were transplanted with 1 × 106 unfractionated WBM cells by retro-orbital injection following a split dose of 10.5 Gy irradiation. For secondary transplantation, 1 × 106 WBM or spleen cells from primary diseased mice were transplanted into sublethally irradiated (6.0 Gy) mice. Peripheral blood counts were performed with a Hemavet 950 (Drew Scientific). Peripheral blood smears and bone marrow and spleen cytospins were stained with the Hema 3 stat pack (Fisher Scientific) and images captured with a Nikon Eclipse E200 microscope equipped with an Infinity 2 color camera (Lumenera) controlled by Infinity Capture software (Lumenera).
Cell purification and flow cytometry
Antibody staining was performed as previously described.19 The following gating strategies were used: HSCs (CD150+ CD48− Lineage− Sca-1+ c-Kit+/Flk2− CD34− Lineage− Sca-1+ c-Kit+), common myeloid progenitors (CMPs) (Lineage− Il7rα− Sca-1low c-Kit+ CD34+ FCγr−), granulocyte-monocyte progenitors (GMPs) (Lineage− Il7rα− Sca-1low c-Kit+ CD34+ FCγr+), megakaryocyte-erythroid progenitors (MEPs) (Lineage− Il7rα− Sca-1low c-Kit+ CD34− FCγr−). Lineage marker cocktail consisted of Gr-1, Mac-1, B220, Ter119, CD4, and CD8. The following antibody (clones) were used (eBioscience or BioLegend): Gr-1 (RB6-8C5), Mac-1 (M1/70), B220 (RA3-6B2), Ter119 (TER119), CD4 (GK1.5), CD8 (53-6.7), Sca-1 (D7), c-Kit (2B8), CD34 (RAM34), Flk2 (A2F10.1), CD150 (TC15-12F12.2), CD48 (HM48-1), CD45.1 (A20), CD45.2 (104), CD71 (R17217), and FcεR1 (MAR-1). Proliferation analysis was performed with the FITC Mouse Anti-Human Ki-67 Set (BD Pharmingen). Apoptosis analysis was performed with the Annexin V Apoptosis Detection Kit APC (eBioscience). Cell sorting and analysis was performed at the Siteman Cancer Center flow cytometry core and the Department of Pathology and Immunology flow cytometry core.
Methocult serial replating
One hundred HSCs were sorted directly into each well of 6-well plates containing Methocult M3434 medium (Stem Cell Technologies) and cultured in vitro at 37°C. Colony-forming units (CFUs) were scored after 7 days, then cells were collected, pooled, and replated at a density of 5000 cells per well.
Plasmids and viral transduction
Mouse c-Kit and c-KitD814V complementary DNAs (cDNAs) were a kind gift of Dr Michael Tomasson (Washington University in St. Louis). The c-KitV750M variant was generated with the QuikChange II XL Site-Directed Mutagenesis Kit (Agilent). All c-Kit variants were subcloned into the HIV-MND-IRES-GFP lentiviral vector as previously described.20 For lentiviral production, 293T cells were cotransfected with the packaging plasmids pMD.G, pRSV-Rev, and PMDLg plus the respective HIV-MND plasmid. Viral supernatant concentrated by centrifugation at 76 000g for 1.5 hours at 4°C. For lentiviral transduction, hematopoietic progenitors were enriched using CD117 microbeads (Miltenyi Biotec). The positive cell fraction was adjusted to 5 × 105 cells per mL in Stempro34 medium (Gibco) supplemented with l-glutamine (2 mM), murine stem cell factor (100 ng/mL), murine thrombopoietin (100 ng/mL), murine Flt3L (50 ng/mL), and murine interleukin-3 (5 ng/mL), and polybrene (4 μg/mL; Sigma), spin-infected with lentivirus at 250g for 2 hours, and transplanted into lethally irradiated mice (100 000 cells per mouse).
DNMT3A cDNA (Open Biosystems) was subcloned into MSCV-IRES-GFP (MIG) vector using Gateway recombination. MIG empty vector was used as a control in retroviral transduction experiments. The DNMT3AR882H variant was generated by site-directed mutagenesis as above. Retroviruses were packaged by cotransfection with pCL-Eco into 293T cells. For retroviral transduction, donor mice were treated with 5-fluorouracil (5-FU; 150 mg/kg; American Pharmaceutical Partners, Schaumburg, IL) 6 days prior to harvest. c-Kit+ cells were selected, spin-infected, and transplanted as above.
32D cell proliferation and western blot
32D cells transduced with lentiviruses were washed with Dulbecco’s modified Eagle’s medium (DMEM) to remove interleukin-3 (IL-3). A total of 50 000 cells were then plated into 24-well plates with DMEM containing 10% fetal bovine serum and 1× penicillin/streptomycin in the presence of either IL-3 or IL-3 and stem cell factor (SCF). Viable cells were counted using Cellometer (Nexcelom Biosciences).
For western blot, transduced 32D cells that were grown until the midexponential phase were washed with DMEM to remove IL-3 and SCF. These cells were resuspended with DMEM containing 10% fetal bovine serum in the absence of IL-3 and SCF for 12 hours. Using 20 µg of whole-cell extracts, western blot was performed to detect phosphorylated or nonphosphorylated Jnk1 and Jnk2 proteins.21
Quantitative real-time PCR
RNA was isolated using the RNAqueous kit (Ambion) and reverse transcribed with the SuperScript VILO kit (Life Technologies). cDNA input was standardized and real-time polymerase chain reaction (PCR) was performed with TaqMan master Mix (Applied Biosystems), 18 s-rRNA probe (VIC-MGB; Applied Biosystems), and a gene-specific probe (FAM-MGB; Applied Biosystems) on a StepOnePlus Real-Time PCR System (Life Technologies). Samples were normalized to 18S and fold change determined by the ∆∆Ct method.
Library construction, capture, and exome sequencing
Genomic DNA was prepared using a PureLink Genomic DNA Mini Kit (Invitrogen). Library construction, capture, and exome sequencing was performed by The Genome Institute at Washington University in St. Louis. Illumina paired-end small-insert multiplexed libraries were constructed according to the manufacturer’s recommendations (Illumina) with the following modifications: (1) 500 ng of native genomic DNA was fragmented using a Covaris E220 DNA Sonicator (Covaris) to range in size between 100 and 400 bp, (2) Illumina adapter-ligated library fragments were amplified in four 50-μL PCR reactions for 18 cycles, (3) solid-phase reversible immobilization bead cleanup was used for enzymatic purification throughout the library process, as well as final library size selection targeting 300-bp to 500-bp fragments. For capture hybridization, 10 indexed libraries were pooled prior to hybridization. Hybridization was performed with the Sure Select Mouse All Exon kit (Agilent). The capture target is 51 Mb. This pool of samples was sequenced on 2 lanes of Illumina HiSequation 2000 (Illumina). Sequence data were aligned to mouse reference sequence build mm9 using bwa version 0.5.9.22 Single-nucleotide variants were detected using the union of 3 callers: samtools version r963,23 VarScan 2.2.6, and Strelka v0.4.6.2.24 Indels were detected using the union of 4 callers: gatk-somatic-indel v5336,25 pindel v0.5,26 VarScan 2.2.6, and Strelka v0.4.6.2.24 Exome sequences are available at the National Center for Biotechnology Information’s sequence read archive under study accession SRP047158.
Student t test and 1-way analysis of variance were used for statistical comparisons where appropriate. Significance is indicated on the figures using the following convention: *P < .05, **P < .01, and ***P < .001. Error bars on all graphs represent the standard error of the mean (SEM) unless indicated otherwise.
Loss of Dnmt3a in HSCs leads to enhanced serial replating capacity
To further examine the kinetics of HSC regulation by Dnmt3a, we performed serial replating CFU assays. Four weeks after pIpC treatment, 100 control (Mx1-CRE:Dnmt3a+/+) or Dnmt3a-KO HSCs (Lineage− c-Kit+ Sca-1+ CD48− CD150+) were plated and scored for CFU after 7 days. A total of 5000 cells were then replated every 7 days. Dnmt3a-KO HSCs generated significantly more colonies after each successive round (Figure 1A), similar to their self-renewal advantage in vivo, which became more dramatic after serial transplantation. Fluorescence-activated cell sorter (FACS) analysis revealed that Dnmt3a-KO colonies possessed a more immature phenotype (Mac-1− Gr-1− c-Kit+ CD34−) after the third plating (Figure 1B). Third-plate Dnmt3a-KO colonies showed increased expression of the self-renewal regulators Meis1 and Evi1 and reduced expression of myeloid-specific factors such as CEBPα and Mpo (Figure 1C). These data suggested that cell turnover is required to enforce the molecular changes resulting from loss of Dnmt3a. In support of this, HSCs from Dnmt3a-KO mice given 2 prior injections of 5-FU 1 month apart (to force HSCs to divide in vivo) showed increased CFU ability compared with control Dnmt3a-KO mice (Figure 1D).
Enforced differentiation of Dnmt3a-KO HSCs leads to bone marrow failure resembling MDS
As it appeared cell turnover was necessary to manifest the Dnmt3a-KO HSC phenotype, we hypothesized that transplantation of Dnmt3a-KO WBM in the absence of WT support WBM would elucidate any predisposition to transformation. We performed noncompetitive transplantation of Dnmt3a-KO WBM 4 weeks following pIpC treatment. Two transplant cohorts were established independently using different donor mice. Control mice in the first cohort were Mx1-CRE-:Dnmt3afl/fl and in the second were Mx1-CRE+:Dnmt3a+/+ to control for interferon-mediated effects, nonspecific deletion of floxed Dnmt3a alleles, and long-term CRE expression in HSCs. As there were no significant differences between the cohorts (supplemental Figure 1, available on the Blood Web site), the final data are compiled from both transplants.
Mice transplanted with Dnmt3a-KO WBM succumbed to bone marrow failure (100% penetrance) ∼1 year posttransplant (median survival, 328 days). No hematopoietic-related mortality was observed from mice transplanted with control or Dnmt3a-heterozygous (Mx1-CRE+:Dnmt3afl/+ = Dnmt3a-HET) WBM when monitored out to 500 days posttransplant (Figure 2A). The majority of the moribund mice (19/25) demonstrated myeloid and erythroid pathologies most consistent with MDS27 (Figure 2B), including peripheral neutropenia and anemia (Figure 2C) and bone marrow hypercellularity (Figure 2D). These mice showed an accumulation of phenotypically defined HSCs in the marrow (Figure 2E-F), although not to the magnitude witnessed in competitive HSC transplants.16 These mice also developed splenomegaly (Figure 2G-H) with myeloid infiltration (Figure 2I-J). Two of the 25 mice developed AML with leukocytosis and >20% myeloid blasts in the blood and bone marrow, whereas 4 out of 25 mice developed diseases with intermediate peripheral anemia and mild leukocytosis with <5% blasts in the bone marrow, most consistent with a mixed MDS/MPN (supplemental Table 1).
We performed a detailed analysis of myeloid cell development in the bone marrow of MDS mice. Although the frequency of CMPs was not different, Dnmt3a-KO mice had increased GMPs and reduced MEPs (Figure 3A). To account for these differences, analysis of proliferation and apoptosis was performed (supplemental Figure 2A). Dnmt3a-KO mice showed increased Ki67+ myeloid progenitors (MPs), but not HSCs, KSL (c-Kit+ Sca-1+ Lineage−) progenitors, or mature (Gr-1+ Mac-1+) myeloid cells (Figure 3B). The frequency of cells undergoing apoptosis within each myeloid cell fraction was determined by Annexin V staining (supplemental Figure 2B). No differences were observed for the proportion of Annexin V+ cells within Dnmt3a-KO KSL or MP populations, but mature myeloid cells from Dnmt3a-KO bone marrow showed significantly higher rates of apoptosis compared with control cells (Figure 3C). Peripheral cytopenias with bone marrow hypercellularity, increased apoptosis, and increased cell proliferation are consistent with an MDS.28
To explore the mechanism leading to anemia, a comparison of erythroid progenitor development between control and Dnmt3a-KO MDS mice was performed. The maturation of erythroid progenitors can be discriminated based on staining for CD71 and Ter119.29 These are from least to most differentiated: proerythroblasts (CD71high Ter119med, stage I), basophilic erythroblasts (CD71high Ter119high, stage II), polychromatophilic erythroblasts (CD71med Ter119high, stage III), and orthochromatophilic erythroblasts (CD71low Ter119high, stage IV). Analysis revealed an accumulation of stage I and stage II progenitors in the spleens of Dnmt3a-KO MDS mice and a significant decrease in the numbers of stage IV progenitors (Figure 3E). Stage I progenitors were also increased in the bone marrow of Dnmt3a-KO MDS mice (Figure 3F). Annexin V staining did not reveal any significant differences in the proportion of apoptotic cells within each cell fraction between control and Dnmt3a-KO mice (Figure 3G-H). This suggests that the accumulation of early erythroid progenitors in the spleens of Dnmt3a-KO MDS mice represents a developmental arrest and not an increase to compensate for abortive loss of later progenitors.
The DNMT3AR882H variant drives myeloproliferation
In myeloid neoplasms, patients with DNMT3A mutations exist in the heterozygous state. The most common mutation occurs at amino acid 882, with over 60% of DNMT3A mutations in AML being the DNMT3AR882H variant.8 The resultant mutant protein acts as a dominant negative by inhibiting oligomerization of higher-order complexes and thereby blocking the ability of DNMT3A to form active tetramers.15 To explore the biological functions of the DNMT3AR882H variant in a genetically relevant system, we transduced Dnmt3a-HET bone marrow with retroviruses carrying green fluorescent protein (GFP)-only control (MIG), WT DNMT3A (MIG-DNMT3A), and DNMT3AR882H (MIG-R882H). To exacerbate the phenotype, we transduced posttransplant Dnmt3a-HET bone marrow.
Ectopic expression of WT DNMT3A had a negative effect on blood differentiation. Expression of DNMT3AR882H produced levels of GFP+ blood chimerism similar to those of control MIG cells (Figure 4A). However, DNMT3AR882H did have a significant effect on the lineage distribution of engrafted cells (Figure 4B), causing a shift toward myeloid differentiation (Figure 4C). Bone marrow analysis 1 year posttransplant showed an expansion of phenotypically defined HSCs in these mice (Figure 4D) and increased numbers of Gr-1+ Mac-1+ cells (Figure 4E), but no overt leukemia. Gene expression analysis confirmed upregulation of Meis1 and HoxA9 induced by DNMT3AR882H30 (Figure 4F). Together, these data suggest DNMT3AR882H drives myeloproliferation and may act as a dominant negative in vivo due to the phenocopy of Dnmt3a loss-of-function alleles.
The dominant clone establishes disease upon secondary transfer
Secondary transplantation of diseased Dnmt3a-KO bone marrow and/or spleen into sublethally irradiated mice was performed to determine if a dominant clone would prevail upon disease reconstitution (and to confirm the phenotypes were cell intrinsic). Secondary transplantation of primary AML quickly recapitulated AML (Figure 5A; median survival, 84 days). Transplantation of primary MDS/MPN revealed that in most secondary mice, AML evolved, albeit with long latency (Figure 5B; median survival 162 days). For Dnmt3a-KO MDS, secondary transfer largely recapitulated the primary disease (median survival, 142 days). Peripheral blood counts of the secondary recipients of individual tumors were strikingly similar to the original tumor they were transplanted with at the time of sacrifice (Figure 5C).
Exome sequencing identifies c-Kit mutations in leukemic transformation
The long latency suggested that acquisition of secondary genetic and/or epigenetic lesions subsequent to inactivation of Dnmt3a were required to produce disease. As we have previously performed extensive DNA methylation profiling,16 here we performed whole-exome sequencing to uncover potential cooperating genetic mutations (Figure 6A). We chose 5 representative primary MDS samples (bone marrow), 1 primary MDS/MPN that progressed to AML in secondary recipients (bone marrow and spleen), and 1 primary AML sample (bone marrow and spleen) for analysis. Diseased spleen cells in secondary mice were also sequenced for select cases to examine clonal evolution. As banked germline DNA from the original bone marrow donor mice was either not available or of insufficient quality/quantity for sequencing to definitively assign the somatic status of prospective mutations, the comparator DNAs for this analysis were derived from pooled bone marrow of recipient mice transplanted with control (Mx1-Cre-:Dnmt3afl/fl) bone marrow from littermate mice. Following capture hybridization, high-throughput sequencing of exome DNA was performed on an Illumina platform.
Initial analysis of primary disease samples identified 55 to 69 single-nucleotide variants and indels for 4 of the 5 MDS samples and 65 and 74 variants for the AML bone marrow and spleen samples, respectively, compared with the pooled reference bone marrow. As paired control DNA was not available to confirm the somatic basis of these mutations, we filtered out any known mouse SNPs identified in the Mouse Phenome Database (Jackson Laboratories) to exclude known strain polymorphisms. One of the 5 primary MDS bone marrow samples and the MDS/MPN spleen sample generated substantially more variant calls (1770 and 2405 respectively) than any other exome, suggesting some off-target capture. Any variants that were exclusively found only in these 2 samples were filtered out. These additional criteria resulted in identification of 16 to 31 high-quality coding variants per sample (supplemental Table 2), a number consistent with other mouse hematopoietic tumor sequencing studies.31 Comparison of these variants to The Cancer Genome Atlas AML database32 highlighted several genes that are also recurrently mutated in AML patients such as Fam5c and Smg1.
We were particularly interested in identifying mutations that were associated with leukemic transformation (present in Dnmt3a-KO AML, but not MDS). Five variants were exclusive to the AML sample, of which 4 were nonsynonymous: HgdP441L (amino acid change 441 P > L), Pet2V359A, IkbkbP551R, and c-KitV750M. Although the variant allele frequencies (VAFs) for these AML-specific variants (VAF range = 6.9% to 20.3%) suggested the pathogenic clone had not achieved complete clonal dominance in the primary sample, the VAFs approximated 50% in the secondary leukemic spleen cells (Figure 6B and supplemental Table 3). Exome sequencing results were confirmed for these variants by Sanger sequencing (Figure 6C). Of particular interest was the c-Kit variant given that AML and mastocytosis patients with DNMT3A mutations can also harbor KIT mutations,32,33 although we could not find evidence for c-Kit mutations in the other de novo AML sample or secondary AML derived from primary MDS/MPN.
c-Kit mutations functionally cooperate with loss of Dnmt3a in vivo
The functional significance of c-KitV750M (identified in our exome sequencing studies) and c-KitD814V (mouse homolog of KITD816V, most common variant found in DNMT3A-mutant AML and mastocytosis) were investigated using lentiviral expression models. Ectopic expression of c-Kit variants in the 32D mouse MP cell line lead to increased proliferation, both at baseline (Figure 7A) and in the presence of the c-Kit ligand SCF (Figure 7B). Analysis of these cells starved of IL-3 and SCF for 12 hours showed increased phosphorylation of Jnk1 and Jnk2 in 32D cells expressing c-KitV750M and c-KitD814V, indicative of ligand-independent activation of c-Kit signaling (Figure 7C). To determine the functional effects of c-Kit variants in primary bone marrow cells, hematopoietic progenitors were isolated and transduced with lentiviruses. Two days posttransduction, 100 GFP+ c-Kit+ Sca-1+ CD150+ cells were sorted and plated per well for CFU assay. Expression of c-KitD814V in a Dnmt3a-deficient background lead to explosive colony generation in the second plate (Figure 7D). c-KitD814V drove mast cell proliferation (Gr-1− Mac-1− c-Kit+ FcεR1+) in both genotypes, but the c-KitV750M variant seemed to have a specific effect in maintaining a more immature phenotype (Gr-1− Mac-1− c-Kit+ FcεR1− Sca-1+ CD150+) specifically in Dnmt3a-KO cells (Figure 7E).
To determine if these 2 pathways could cooperate in vivo, we transduced hematopoietic progenitor cells with c-Kit variants and transplanted into lethally irradiated mice. As previously reported, expression of c-KitD814V in WT cells lead to development of B-cell acute lymphoblastic leukemia (B-ALL).34 However, c-KitD814V in a Dnmt3a-KO background led to explosive leukemia with a much shorter latency (Figure 7F). Moreover, although many of the mice transplanted with Dnmt3a-KO c-KitD814V cells also succumbed to a B-ALL, 4 out of 13 (31%) developed mastocytosis with involvement of myeloid blasts, and 4 out of 13 (31%) mice developed a T-cell acute lymphoblastic leukemia (T-ALL; Figure 7G). Expression of the c-KitV750M mutant lead to a selective advantage seen only in a Dnmt3a-deficient background, with 2 out of 8 mice (25%) succumbing to B-ALL within 6 months posttransplant, while all mice transplanted with WT cells transduced with c-KitV750M (8/8) were alive and healthy at this time point. Of the surviving mice, analysis of the lineage distribution of transduced donor-derived (CD45.2+ GFP+) peripheral blood cells showed a relative increase in the proportion of B cells in Dnmt3a-KO c-KitV750M cells (Figure 7H). These data suggest that this particular mutation may only be competent for transformation in a Dnmt3a-deficient background.
Our data highlight crucial roles of Dnmt3a in normal and malignant hematopoiesis. Loss of Dnmt3a leads to developmental arrest of erythroid progenitors, enhanced proliferation of MPs, and increased apoptosis of mature myeloid cells. Although virtually all of the reported DNMT3A mutations in myeloid neoplasms are heterozygous,8,9,32 we did not observe a phenotype from Dnmt3a-heterozygous bone marrow. The most prevalent mutation in AML is the DNMT3AR882H variant, which acts as a dominant negative,14,15 making the patients effectively null for DNMT3A, akin to our genetic mouse model. The remainder of patients with heterozygous DNMT3A mutations typically present at an advanced age,8,32,35 suggesting the time to transformation may be dosage dependent in DNMT3A haploinsufficiency. The timeframe of our mouse experiments may not have been sufficient to observe this, and it is possible a phenotype for Dnmt3a-heterozygous bone marrow could be unveiled by extended transplantation or subjection to other stressors of HSC turnover such as serial 5-FU injection.
Activating mutations of KIT are associated with AML,36 germ cell tumors,37 and systemic mastocytosis.38 Exome sequencing identified acquisition of a c-Kit mutation (c-KitV750M) that was exclusively identified during leukemic transformation of Dnmt3a-KO cells. Characterization of this variant suggests it acts as a gain-of-function mutation, as it conferred a similar growth advantage and ligand-independent activation of downstream signaling pathways as the prototypical c-KitD814V variant. Functional studies showed synergism between Dnmt3a loss of function and c-Kit gain of function in generating an explosive leukemia with a much shorter latency than expression of c-KitD814V in WT bone marrow. Although DNMT3A and KIT have been shown to be concurrently mutated in mastocytosis33 and AML,32 Dnmt3a-KO c-KitV750M transplanted mice developed B-ALL. This is likely due to ectopic introduction of c-KitV750M into a lineage-biased HSC or progenitor cell. Similar to the BCR-ABL oncogene,39 the type of malignancy generated is likely dependent on the nature of the initiating cell that successfully integrates the virus, as evidenced by the spectrum of different diseases generated from c-KitD814V in Dnmt3a-KO transplanted mice (B-ALL, T-ALL, mastocytosis). To fully understand the relationship between Dnmt3a and activating c-Kit mutations in myeloid leukemogenesis, a genetic mouse model of regulatable stem cell–specific induction of c-Kit activation in Dnmt3a-KO HSCs will be required. However, here we show for the first time that these pathways can cooperate to accelerate transformation in vivo. This Dnmt3a/c-Kit disease model resembles the classical “two-hit” model of leukemogenesis40 in which one mutation inhibits differentiation (Dnmt3a loss-of-function), while another drives proliferation (c-Kit gain-of-function). Although we did not identify other mutations typically associated with DNMT3A-mutant AML patients (eg, FLT3-ITD, NPM1c), only approximately 10% of primary mice (2/25) developed de novo AML, and sequencing of many more samples would be required to fully ascribe to repertoire of cooperating mutations in this mouse model. As cooperating DNMT3A and KIT mutations are relatively rare in AML, the observation that this Dnmt3a mouse model selected for this combination in vivo likely represents stochastic selection and evolution at the level of the HSC.
We show that enforced differentiation of Dnmt3a-KO HSCs by noncompetitive transplantation leads to bone marrow failure resembling MDS. This is markedly distinct from our previous competitive transplantation outcomes and shows the critical requirement of cell division to acquire secondary genetic and epigenetic lesions necessary to promote disease. Such systems present a unique opportunity to study the sequence of early events leading to HSC transformation following Dnmt3a loss of function.
The online version of this article contains a data supplement.
There is an Inside Blood Commentary on this article in this issue.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
The authors thank Dr Michael Tomasson (Washington University School of Medicine [WUSM]) for c-Kit plasmids, Dr Timothy Ley and the Genome Institute (WUSM) for exome sequencing, Dr Eric Duncavage (WUSM) for pathology imaging, and the Alvin J. Siteman Cancer Center at WUSM for the use of the Siteman flow cytometry core, which provided cell sorting and analysis.
The Siteman Cancer Center is supported in part by the National Institutes of Health, National Cancer Institute Cancer Center Support Grant CA91842. This work was supported by the National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases (DK084259), the Children’s Discovery Institute, the Mallinckrodt Foundation, Alex's Lemonade Stand Foundation, and an American Society of Hematology Scholar award (all to G.A.C.).
Contribution: H.C., C.M., A.K., E.L.O., E.E., A.M., and G.A.C. designed and performed experiments; H.C., C.A.M., J.H., J.M.K., and G.A.C. analyzed data; and G.A.C. wrote and edited the paper.
Conflict-of-interest disclosure: The authors declare no competing financial interests.
The current affiliation for J.M.K. is St. Jude Children’s Research Hospital, Memphis, TN, 38105.
Correspondence: Grant A. Challen, Washington University School of Medicine, 660 Euclid Ave, St. Louis, MO 63110; e-mail: gchallen@DOM.wustl.edu.