Characterization of the bone marrow niche in patients with chronic myeloid leukemia identifies CXCL14 as a new therapeutic option.

Although tyrosine kinase inhibitors (TKIs) are effective in treating chronic myeloid leukemia (CML), they often fail to eradicate the leukemia-initiating stem cells (LSCs), causing disease persistence and relapse. Evidence indicates that LSC persistence may be because of bone marrow (BM) niche protection; however, little is known about the underlying mechanisms. Herein, we molecularly and functionally characterize BM niches in patients with CML at diagnosis and reveal the altered niche composition and function in these patients. Long-term culture initiating cell assay showed that the mesenchymal stem cells from patients with CML displayed an enhanced supporting capacity for normal and CML BM CD34+CD38- cells. Molecularly, RNA sequencing detected dysregulated cytokine and growth factor expression in the BM cellular niches of patients with CML. Among them, CXCL14 was lost in the BM cellular niches in contrast to its expression in healthy BM. Restoring CXCL14 significantly inhibited CML LSC maintenance and enhanced their response to imatinib in vitro, and CML engraftment in vivo in NSG-SGM3 mice. Importantly, CXCL14 treatment dramatically inhibited CML engraftment in patient-derived xenografted NSG-SGM3 mice, even to a greater degree than imatinib, and this inhibition persisted in patients with suboptimal TKI response. Mechanistically, CXCL14 upregulated inflammatory cytokine signaling but downregulated mTOR signaling and oxidative phosphorylation in CML LSCs. Together, we have discovered a suppressive role of CXCL14 in CML LSC growth. CXCL14 might offer a treatment option targeting CML LSCs.

Characterization of therapy-related acute myeloid leukemia: increasing incidence and prognostic implications.

Studies of therapy-related AML (t-AML) are usually performed in selected cohorts and reliable incidence rates are lacking. In this study, we characterized, defined the incidence over time and studied prognostic implications in all t-AML patients diagnosed in Sweden between 1997 and 2015. Data were retrieved from nationwide population-based registries. In total, 6,779 AML patients were included in the study, of whom 686 (10%) had t-AML. The median age for t-AML was 71 years and 392 (57%) patients were females. During the study period, the incidence of t-AML almost doubled with a yearly increase in t-AML of 4.5% (95% confidence interval: 2.8%-6.2%), which contributed significantly to the general increase in AML incidence over the study period. t-AML solidly constituted over 10% of all AML cases during the later period of the study. Primary diagnoses with the largest increase in incidence and decrease in mortality rate during the study period (i.e., breast and prostate cancer) contributed significantly to the increased incidence of t-AML. In multivariable analysis, t-AML was associated with poorer outcome in cytogenetically intermediate- and adverse-risk cases but t-AML had no significant impact on outcome in favorable-risk AML, including core binding leukemias, acute promyelocytic leukemia and AML with mutated NPM1 without FLT3-ITD. We conclude that there is a strong increase in incidence in t-AML over time and that t-AML constitutes a successively larger proportion of the AML cases. Furthermore, we conclude that t-AML confers a poor prognosis in cytogenetically intermediate- and adverse-risk, but not in favorable-risk AML.

Targeting SAMHD1 with hydroxyurea in first-line cytarabine-based therapy of newly diagnosed acute myeloid leukaemia: Results from the HEAT-AML trial.

BACKGROUND: Treatment of newly diagnosed acute myeloid leukaemia (AML) is based on combination chemotherapy with cytarabine (ara-C) and anthracyclines. Five-year overall survival is below 30%, which has partly been attributed to cytarabine resistance. Preclinical data suggest that the addition of hydroxyurea potentiates cytarabine efficacy by increasing ara-C triphosphate (ara-CTP) levels through targeted inhibition of SAMHD1. OBJECTIVES: In this phase 1 trial, we evaluated the feasibility, safety and efficacy of the addition of hydroxyurea to standard chemotherapy with cytarabine/daunorubicin in newly diagnosed AML patients. METHODS: Nine patients were enrolled and received at least two courses of ara-C (1 g/m2 /2 h b.i.d. d1-5, i.e., a total of 10 g/m2 per course), hydroxyurea (1-2 g d1-5) and daunorubicin (60 mg/m2 d1-3). The primary endpoint was safety; secondary endpoints were complete remission rate and measurable residual disease (MRD). Additionally, pharmacokinetic studies of ara-CTP and ex vivo drug sensitivity assays were performed. RESULTS: The most common grade 3-4 toxicity was febrile neutropenia (100%). No unexpected toxicities were observed. Pharmacokinetic analyses showed a significant increase in median ara-CTP levels (1.5-fold; p = 0.04) in patients receiving doses of 1 g hydroxyurea. Ex vivo, diagnostic leukaemic bone marrow blasts from study patients were significantly sensitised to ara-C by a median factor of 2.1 (p = 0.0047). All nine patients (100%) achieved complete remission, and all eight (100%) with validated MRD measurements (flow cytometry or real-time quantitative polymerase chain reaction [RT-qPCR]) had an MRD level <0.1% after two cycles of chemotherapy. Treatment was well-tolerated, and median time to neutrophil recovery >1.0 × 109 /L and to platelet recovery >50 × 109 /L after the start of cycle 1 was 19 days and 22 days, respectively. Six of nine patients underwent allogeneic haematopoietic stem-cell transplantation (allo-HSCT). With a median follow-up of 18.0 (range 14.9-20.5) months, one patient with adverse risk not fit for HSCT experienced a relapse after 11.9 months but is now in second complete remission. CONCLUSION: Targeted inhibition of SAMHD1 by the addition of hydroxyurea to conventional AML therapy is safe and appears efficacious within the limitations of the small phase 1 patient cohort. These results need to be corroborated in a larger study.

AML displays increased CTCF occupancy associated with aberrant gene expression and transcription factor binding.

CCTC-binding factor (CTCF) is a key regulator of gene expression through organization of the chromatin structure. Still, it is unclear how CTCF binding is perturbed in leukemia or in cancer in general. We studied CTCF binding by chromatin immunoprecipitation sequencing in cells from patients with acute myeloid leukemia (AML) and in normal bone marrow (NBM) in the context of gene expression, DNA methylation, and azacitidine exposure. CTCF binding was increased in AML compared with NBM. Aberrant CTCF binding was enriched for motifs for key myeloid transcription factors such as CEBPA, PU.1, and RUNX1. AML with TET2 mutations was characterized by a particularly strong gain of CTCF binding, highly enriched for gain in promoter regions, while AML in general was enriched for changes at enhancers. There was a strong anticorrelation between CTCF binding and DNA methylation. Gain of CTCF occupancy was associated with increased gene expression; however, the genomic location (promoter vs distal regions) and enrichment of motifs (for repressing vs activating cofactors) were decisive for the gene expression pattern. Knockdown of CTCF in K562 cells caused loss of CTCF binding and transcriptional repression of genes with changed CTCF binding in AML, as well as loss of RUNX1 binding at RUNX1/CTCF-binding sites. In addition, CTCF knockdown caused increased differentiation. Azacitidine exposure caused major changes in CTCF occupancy in AML patient cells, partly by restoring a CTCF-binding pattern similar to NBM. We conclude that AML displays an aberrant increase in CTCF occupancy that targets key genes for AML development and impacts gene expression.

A novel CD34-specific T-cell engager efficiently depletes acute myeloid leukemia and leukemic stem cells in vitro and in vivo.

Less than a third of patients with acute myeloid leukemia (AML) are cured by chemotherapy and/or hematopoietic stem cell transplantation, highlighting the need to develop more efficient drugs. The low efficacy of standard treatments is associated with inadequate depletion of CD34+ blasts and leukemic stem cells, the latter a drug-resistant subpopulation of leukemia cells characterized by the CD34+CD38- phenotype. To target these drug-resistant primitive leukemic cells better, we have designed a CD34/CD3 bi-specific T-cell engager (BTE) and characterized its anti-leukemia potential in vitro, ex vivo and in vivo. Our results show that this CD34-specific BTE induces CD34-dependent T-cell activation and subsequent leukemia cell killing in a dose-dependent manner, further corroborated by enhanced T-cell-mediated killing at the singlecell level. Additionally, the BTE triggered efficient T-cell-mediated depletion of CD34+ hematopoietic stem cells from peripheral blood stem cell grafts and CD34+ blasts from AML patients. Using a humanized AML xenograft model, we confirmed that the CD34-specific BTE had in vivo efficacy by depleting CD34+ blasts and leukemic stem cells without side effects. Taken together, these data demonstrate that the CD34-specific BTE has robust antitumor effects, supporting development of a novel treatment modality with the aim of improving outcomes of patients with AML and myelodysplastic syndromes.

A risk score based on real-world data to predict early death in acute promyelocytic leukemia.

With increasingly effective treatments, early death (ED) has become the predominant reason for therapeutic failure in patients with acute promyelocytic leukemia (APL). To better prevent ED, patients with high-risk of ED must be identified. Our aim was to develop a score that predicts the risk of ED in a real-life setting. We used APL patients in the populationbased Swedish AML Registry (n=301) and a Portuguese hospital-based registry (n=129) as training and validation cohorts, respectively. The cohorts were comparable with respect to age (median, 54 and 53 years) and ED rate (19.6% and 18.6%). The score was developed by logistic regression analyses, risk-per-quantile assessment and scoring based on ridge regression coefficients from multivariable penalized logistic regression analysis. White blood cell count, platelet count and age were selected by this approach as the most significant variables for predicting ED. The score identified low-, high- and very high-risk patients with ED risks of 4.8%, 20.2% and 50.9% respectively in the training cohort and with 6.7%, 25.0% and 36.0% as corresponding values for the validation cohort. The score identified an increased risk of ED already at sub-normal and normal white blood cell counts and, consequently, it was better at predicting ED risk than the Sanz score (AUROC 0.77 vs. 0.64). In summary, we here present an externally validated and population-based risk score to predict ED risk in a real-world setting, identifying patients with the most urgent need of aggressive ED prevention. The results also suggest that increased vigilance for ED is already necessary at sub-normal/normal white blood cell counts.

Clinical and genomic characterization of patients diagnosed with the provisional entity acute myeloid leukemia with BCR-ABL1, a Swedish population-based study.

Acute myeloid leukemia (AML) with t(9;22)(q34;q11), also known as AML with BCR-ABL1, is a rare, provisional entity in the WHO 2016 classification and is considered a high-risk disease according to the European LeukemiaNet 2017 risk stratification. We here present a retrospective, population-based study of this disease entity from the Swedish Acute Leukemia Registry. By strict clinical inclusion criteria we aimed to identify genetic markers further distinguishing AML with t(9;22) as a separate entity. Twenty-five patients were identified and next-generation sequencing using a 54-gene panel was performed in 21 cases. Interestingly, no mutations were found in NPM1, FLT3, or DNMT3A, three frequently mutated genes in AML. Instead, RUNX1 was the most commonly mutated gene, with aberrations present in 38% of the cases compared to around 10% in de novo AML. Additional mutations were identified in genes involved in RNA splicing (SRSF2, SF3B1) and chromatin regulation (ASXL1, STAG2, BCOR, BCORL1). Less frequently, mutations were found in IDH2, NRAS, TET2, and TP53. The mutational landscape exhibited a similar pattern as recently described in patients with chronic myeloid leukemia (CML) in myeloid blast crisis (BC). Despite the concomitant presence of BCR-ABL1 and RUNX1 mutations in our cohort, both features of high-risk AML, the RUNX1-mutated cases showed a superior overall survival compared to RUNX1 wildtype cases. Our results suggest that the molecular characteristics of AML with t(9;22)/BCR-ABL1 and CML in myeloid BC are similar and do not support a distinction of the two disease entities based on their underlying molecular alterations.

Prognostic scoring systems and comorbidities in chronic myelomonocytic leukaemia: a nationwide population-based study.

Outcomes in chronic myelomonocytic leukaemia (CMML) are highly variable and may be affected by comorbidity. Therefore, prognostic models and comorbidity indices are important tools to estimate survival and to guide clinicians in individualising treatment. In this nationwide population-based study, we assess comorbidities and for the first time validate comorbidity indices in CMML. We also compare the prognostic power of: the revised International Prognostic Scoring System (IPSS-R), CMML-specific prognostic scoring system (CPSS), MD Anderson Prognostic Scoring System (MDAPS) and Mayo score. In this cohort of 337 patients with CMML, diagnosed between 2009 and 2015, the median overall survival was 21·3 months. Autoimmune conditions were present in 25% of the patients, with polymyalgia rheumatica and Hashimoto's thyroiditis being most common. Of the tested comorbidity indices: the Charlson Comorbidity Index (CCI), Haematopoietic cell transplantation-specific Comorbidity Index (HCT-CI) and Myelodysplastic Syndrome-Specific Comorbidity Index (MDS-CI), CCI had the highest C-index (0·62) and was the only comorbidity index independently associated with survival in multivariable analyses. When comparing the prognostic power of the scoring systems, the CPSS had the highest C-index (0·69). In conclusion, using 'real-world' data we found that the CCI and CPSS have the best prognostic power and that autoimmune conditions are overrepresented in CMML.

Combinatorial molecule screening identified a novel diterpene and the BET inhibitor CPI-203 as differentiation inducers of primary acute myeloid leukemia cells.

Combination treatment has proven effective for patients with acute promyelocytic leukemia, exemplifying the importance of therapy targeting multiple components of oncogenic regulation for a successful outcome. However, recent studies have shown that the mutational complexity of acute myeloid leukemia (AML) precludes the translation of molecular targeting into clinical success. Here, as a complement to genetic profiling, we used unbiased, combinatorial in vitro drug screening to identify pathways that drive AML and to develop personalized combinatorial treatments. First, we screened 513 natural compounds on primary AML cells and identified a novel diterpene (H4) that preferentially induced differentiation of FLT3 wild-type AML, while FLT3-ITD/mutations conferred resistance. The samples responding to H4, displayed increased expression of myeloid markers, a clear decrease in the nuclear-cytoplasmic ratio and the potential of re-activation of the monocytic transcriptional program reducing leukemia propagation in vivo. By combinatorial screening using H4 and molecules with defined targets, we demonstrated that H4 induces differentiation by the activation of the protein kinase C (PKC) signaling pathway, and in line with this, activates PKC phosphorylation and translocation of PKC to the cell membrane. Furthermore, the combinatorial screening identified a bromo- and extra-terminal domain (BET) inhibitor that could further improve H4-dependent leukemic differentiation in FLT3 wild-type monocytic AML. These findings illustrate the value of an unbiased, multiplex screening platform for developing combinatorial therapeutic approaches for AML.

The transcriptome-wide landscape of molecular subtype-specific mRNA expression profiles in acute myeloid leukemia.

Molecular classification of acute myeloid leukemia (AML) aids prognostic stratification and clinical management. Our aim in this study is to identify transcriptome-wide mRNAs that are specific to each of the molecular subtypes of AML. We analyzed RNA-sequencing data of 955 AML samples from three cohorts, including the BeatAML project, the Cancer Genome Atlas, and a cohort of Swedish patients to provide a comprehensive transcriptome-wide view of subtype-specific mRNA expression. We identified 729 subtype-specific mRNAs, discovered in the BeatAML project and validated in the other two cohorts. Using unique proteomics data, we also validated the presence of subtype-specific mRNAs at the protein level, yielding a rich collection of potential protein-based biomarkers for the AML community. To enable the exploration of subtype-specific mRNA expression by the broader scientific community, we provide an interactive resource to the public.

Secondary Acute Myeloid Leukemia and the Role of Allogeneic Stem Cell Transplantation in a Population-Based Setting.

Secondary AML (s-AML), including AML with an antecedent hematologic disorder (AHD-AML) and therapy-related AML (t-AML), constitutes a large proportion of patients with AML and is considered to confer a dismal prognosis. The role of allogeneic hematopoietic cell transplantation (HCT) in patients with s-AML and the extent to which HCT is performed in these patients has been little studied to date. We used the population-based Swedish AML Registry comprising 3337 intensively treated adult patients over a 17-year period to study the role of HCT within the group of patients with s-AML as well as compared with patients with de novo AML. HCT was performed in 576 patients (22%) with de novo AML, in 74 patients (17%) with AHD-AML, and in 57 patients (20%) with t-AML. At 5 years after diagnosis, there were no survivors among patients with previous myeloproliferative neoplasms who did not undergo HCT, and corresponding survival for patients with antecedent myelodysplastic syndromes and t-AML was and 2% and 4%, respectively. HCT was compared with chemotherapy consolidation in s-AML using 3 models: (1) a 200-day landmark analysis, in which HCT was favorable compared with conventional consolidation (P = .04, log-rank test); (2) a multivariable Cox regression with HCT as a time-dependent variable, in which the hazard ratio for mortality was 0.73 (95% confidence interval, 0.64 to 0.83) for HCT and favored HCT in all subgroups; and (3) a propensity score matching analysis, in which the 5-year overall survival (OS) and relapse-free survival in patients with s-AML in first complete remission (CR1) was 48% and 43%, respectively, for patients undergoing HCT versus 20% and 21%, respectively, for those receiving chemotherapy consolidation (P = .01 and .02, respectively, log-rank test). Our observational data suggest that HCT improves survival and offers the only realistic curative treatment option in patients with s-AML.

Cancer-specific changes in DNA methylation reveal aberrant silencing and activation of enhancers in leukemia.

Acute myeloid leukemia (AML) is characterized by an impaired differentiation process leading to an accumulation of immature blasts in the blood. One feature of cytogenetically normal AML is alterations to the DNA methylome. We analyzed 57 AML patients with normal karyotype by using Illumina's 450k array and showed that aberrant DNA methylation is significantly altered at enhancer regions and that the methylation levels at specific enhancers predict overall survival of AML patients. The majority of sites that become differentially methylated in AML occur in regulatory elements of the human genome. Hypermethylation associates with enhancer silencing. In addition, chromatin immunoprecipitation sequencing analyses showed that a subset of hypomethylated sites correlate with enhancer activation, indicated by increased H3K27 acetylation. DNA hypomethylation is therefore not sufficient for enhancer activation. Some sites of hypomethylation occur at weak/poised enhancers marked with H3K4 monomethylation in hematopoietic progenitor cells. Other hypomethylated regions occur at sites inactive in progenitors and reflect the de novo acquisition of AML-specific enhancers. Altered enhancer dynamics are reflected in the gene expression of enhancer target genes, including genes involved in oncogenesis and blood cell development. This study demonstrates that histone variants and different histone modifications interact with aberrant DNA methylation and cause perturbed enhancer activity in cytogenetically normal AML that contributes to a leukemic transcriptome.

Prognostic scoring systems for myelodysplastic syndromes (MDS) in a population-based setting: a report from the Swedish MDS register.

The myelodysplastic syndromes (MDS) have highly variable outcomes and prognostic scoring systems are important tools for risk assessment and to guide therapeutic decisions. However, few population-based studies have compared the value of the different scoring systems. With data from the nationwide Swedish population-based MDS register we validated the International Prognostic Scoring System (IPSS), revised IPSS (IPSS-R) and the World Health Organization (WHO) Classification-based Prognostic Scoring System (WPSS). We also present population-based data on incidence, clinical characteristics including detailed cytogenetics and outcome from the register. The study encompassed 1329 patients reported to the register between 2009 and 2013, 14% of these had therapy-related MDS (t-MDS). Based on the MDS register, the yearly crude incidence of MDS in Sweden was 2·9 per 100 000 inhabitants. IPSS-R had a significantly better prognostic power than IPSS (P < 0·001). There was a trend for better prognostic power of IPSS-R compared to WPSS (P = 0·05) and for WPSS compared to IPSS (P = 0·07). IPSS-R was superior to both IPSS and WPSS for patients aged ≤70 years. Patients with t-MDS had a worse outcome compared to de novo MDS (d-MDS), however, the validity of the prognostic scoring systems was comparable for d-MDS and t-MDS. In conclusion, population-based studies are important to validate prognostic scores in a 'real-world' setting. In our nationwide cohort, the IPSS-R showed the best predictive power.

Clonal hematopoiesis and blood-cancer risk inferred from blood DNA sequence.

BACKGROUND: Cancers arise from multiple acquired mutations, which presumably occur over many years. Early stages in cancer development might be present years before cancers become clinically apparent. METHODS: We analyzed data from whole-exome sequencing of DNA in peripheral-blood cells from 12,380 persons, unselected for cancer or hematologic phenotypes. We identified somatic mutations on the basis of unusual allelic fractions. We used data from Swedish national patient registers to follow health outcomes for 2 to 7 years after DNA sampling. RESULTS: Clonal hematopoiesis with somatic mutations was observed in 10% of persons older than 65 years of age but in only 1% of those younger than 50 years of age. Detectable clonal expansions most frequently involved somatic mutations in three genes (DNMT3A, ASXL1, and TET2) that have previously been implicated in hematologic cancers. Clonal hematopoiesis was a strong risk factor for subsequent hematologic cancer (hazard ratio, 12.9; 95% confidence interval, 5.8 to 28.7). Approximately 42% of hematologic cancers in this cohort arose in persons who had clonality at the time of DNA sampling, more than 6 months before a first diagnosis of cancer. Analysis of bone marrow-biopsy specimens obtained from two patients at the time of diagnosis of acute myeloid leukemia revealed that their cancers arose from the earlier clones. CONCLUSIONS: Clonal hematopoiesis with somatic mutations is readily detected by means of DNA sequencing, is increasingly common as people age, and is associated with increased risks of hematologic cancer and death. A subset of the genes that are mutated in patients with myeloid cancers is frequently mutated in apparently healthy persons; these mutations may represent characteristic early events in the development of hematologic cancers. (Funded by the National Human Genome Research Institute and others.).

Incidence and prognostic significance of isolated trisomies in adult acute myeloid leukemia: A population-based study from the Swedish AML registry.

OBJECTIVES AND METHODS: To ascertain the incidence/clinical implications of isolated autosomal trisomies in adult acute myeloid leukemia (AML), all such cases were retrieved from the Swedish AML Registry. RESULTS: Of the 3179 cytogenetically informative AMLs diagnosed January 1997-May 2015, 246 (7.7%) had isolated trisomies. The frequency increased by age (2.4% at age 18-60 years vs. 23% at >60 years; P<.0001); the median age was 69 years. The five most common were +8 (4.0%), +13 (0.9%), +11 (0.8%), +21 (0.7%), and +4 (0.5%). Age and gender, types of AML and treatment, and complete remission and early death rates did not differ between the single trisomy and the intermediate risk (IR) groups or among cases with isolated gains of chromosomes 4, 8, 11, 13, or 21. The overall survival (OS) was similar in the single trisomy (median 1.6 years) and IR groups (1.7 years; P=.251). The OS differed among the most frequent isolated trisomies; the median OS was 2.5 years for +4, 1.9 years for +21, 1.5 years for +8, 1.1 years for +11, and 0.8 years for +13 (P=.013). CONCLUSION: AML with single trisomies, with the exception of +13, should be grouped as IR.

Anti-leukaemic effects induced by APR-246 are dependent on induction of oxidative stress and the NFE2L2/HMOX1 axis that can be targeted by PI3K and mTOR inhibitors in acute myeloid leukaemia cells.

The small molecule APR-246 (PRIMA-1(MET) ) is a novel drug that restores the activity of mutated and unfolded TP53 protein. However, the mechanisms of action and potential off-target effects are not fully understood. Gene expression profiling in TP53 mutant KMB3 acute myeloid leukaemia (AML) cells showed that genes which protected cells from oxidative stress to be the most up-regulated. APR-246 exposure also induced reactive oxygen species (ROS) formation and depleted glutathione in AML cells. The genes most up-regulated by APR-246, confirmed by quantitative real time polymerase chain reaction, were heme oxygenase-1 (HMOX1, also termed HO-1), SLC7A11 and RIT1. Up-regulation of HMOX1, a key regulator of cellular response to ROS, was independent of TP53 mutational status. NFE2L2 (also termed Nrf2), a master regulator of HMOX1 expression, showed transcriptional up-regulation and nuclear translocation by APR-246. Down-regulation of NFE2L2 by siRNA in AML cells significantly increased the antitumoural effects of APR-246. The PI3K inhibitor wortmannin and the mTOR inhibitor rapamycin inhibited APR-246-induced nuclear translocation of NFE2L2 and counteracted the protective cellular responses to APR-246, resulting in synergistic cell killing together with APR-246. In conclusion, ROS induction is important for antileukaemic activities of APR-246 and inhibiting the protective response of the Nrf-2/HMOX1 axis using PI3K inhibitors, enhances the antileukaemic effects.

Analysis of the DNA methylome and transcriptome in granulopoiesis reveals timed changes and dynamic enhancer methylation.

In development, epigenetic mechanisms such as DNA methylation have been suggested to provide a cellular memory to maintain multipotency but also stabilize cell fate decisions and direct lineage restriction. In this study, we set out to characterize changes in DNA methylation and gene expression during granulopoiesis using 4 distinct cell populations ranging from the oligopotent common myeloid progenitor stage to terminally differentiated neutrophils. We observed that differentially methylated sites (DMSs) generally show decreased methylation during granulopoiesis. Methylation appears to change at specific differentiation stages and overlap with changes in transcription and activity of key hematopoietic transcription factors. DMSs were preferentially located in areas distal to CpG islands and shores. Also, DMSs were overrepresented in enhancer elements and enriched in enhancers that become active during differentiation. Overall, this study depicts in detail the epigenetic and transcriptional changes that occur during granulopoiesis and supports the role of DNA methylation as a regulatory mechanism in blood cell differentiation.