Multi-ancestry epigenome-wide analyses identify methylated sites associated with aortic augmentation index in TOPMed MESA

Despite the prognostic value of arterial stiffness (AS) and pulsatile hemodynamics (PH) for cardiovascular morbidity and mortality, epigenetic modifications that contribute to AS/PH remain unknown. To gain a better understanding of the link between epigenetics (DNA methylation) and AS/PH, we examined the relationship of eight measures of AS/PH with CpG sites and co-methylated regions using multi-ancestry participants from Trans-Omics for Precision Medicine (TOPMed) Multi-Ethnic Study of Atherosclerosis (MESA) with sample sizes ranging from 438 to 874. Epigenome-wide association analysis identified one genome-wide significant CpG (cg20711926-CYP1B1) associated with aortic augmentation index (AIx). Follow-up analyses, including gene set enrichment analysis, expression quantitative trait methylation analysis, and functional enrichment analysis on differentially methylated positions and regions, further prioritized three CpGs and their annotated genes (cg23800023-ETS1, cg08426368-TGFB3, and cg17350632-HLA-DPB1) for AIx. Among these, ETS1 and TGFB3 have been previously prioritized as candidate genes. Furthermore, both ETS1 and HLA-DPB1 have significant tissue correlations between Whole Blood and Aorta in GTEx, which suggests ETS1 and HLA-DPB1 could be potential biomarkers in understanding pathophysiology of AS/PH. Overall, our findings support the possible role of epigenetic regulation via DNA methylation of specific genes associated with AIx as well as identifying potential targets for regulation of AS/PH.

Depletion of tet2 results in age-dependent changes in DNA methylation and gene expression in a zebrafish model of myelodysplastic syndrome

Introduction

Myelodysplastic syndrome (MDS) is a heterogeneous group of clonal hematopoietic disorders characterized by ineffective hematopoiesis, cytopenias, and dysplasia. The gene encoding ten-eleven translocation 2 (tet2), a dioxygenase enzyme that catalyzes the conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine, is a recurrently mutated tumor suppressor gene in MDS and other myeloid malignancies. Previously, we reported a stable zebrafish line with a loss-of-function mutation in the tet2 gene. The tet2m/m-mutant zebrafish developed a pre-MDS state with kidney marrow dysplasia, but normal circulating blood counts by 11 months of age and accompanying anemia, signifying the onset of MDS, by 24 months of age.

Methods

In the current study, we collected progenitor cells from the kidney marrows of the adult tet2m/m and tet2wt/wt fish at 4 and 15 months of age and conducted enhanced reduced representation of bisulfite sequencing (ERRBS) and bulk RNA-seq to measure changes in DNA methylation and gene expression of hematopoietic stem and progenitor cells (HSPCs).

Results and discussion

A global increase in DNA methylation of gene promoter regions and CpG islands was observed in tet2m/m HSPCs at 4 months of age when compared with the wild type. Furthermore, hypermethylated genes were significantly enriched for targets of SUZ12 and the metal-response-element-binding transcription factor 2 (MTF2)-involved in the polycomb repressive complex 2 (PRC2). However, between 4 and 15 months of age, we observed a paradoxical global decrease in DNA methylation in tet2m/m HSPCs. Gene expression analyses identified upregulation of genes associated with mTORC1 signaling and interferon gamma and alpha responses in tet2m/m HSPCs at 4 months of age when compared with the wild type. Downregulated genes in HSPCs of tet2-mutant fish at 4 months of age were enriched for cell cycle regulation, heme metabolism, and interleukin 2 (IL2)/signal transducer and activator of transcription 5 (STAT5) signaling, possibly related to increased self-renewal and clonal advantage in HSPCs with tet2 loss of function. Finally, there was an overall inverse correlation between overall increased promoter methylation and gene expression.

Acute myeloid leukemia stratifies as two clinically relevant sphingolipidomic subtypes

Acute myeloid leukemia (AML) is an aggressive disease with complex and heterogeneous biology. Although several genomic classifications have been proposed, there is a growing interest in going beyond genomics to stratify AML. In this study, we profile the sphingolipid family of bioactive molecules in 213 primary AML samples and 30 common human AML cell lines. Using an integrative approach, we identify two distinct sphingolipid subtypes in AML characterized by a reciprocal abundance of hexosylceramide (Hex) and sphingomyelin (SM) species. The two Hex-SM clusters organize diverse samples more robustly than known AML driver mutations and are coupled to latent transcriptional states. Using transcriptomic data, we develop a machine-learning classifier to infer the Hex-SM status of AML cases in TCGA and BeatAML clinical repositories. The analyses show that the sphingolipid subtype with deficient Hex and abundant SM is enriched for leukemic stemness transcriptional programs and comprises an unappreciated high-risk subgroup with poor clinical outcomes. Our sphingolipid-focused examination of AML identifies patients least likely to benefit from standard of care and raises the possibility that sphingolipidomic interventions could switch the subtype of AML patients who otherwise lack targetable alternatives.

Abstract 3155: Gene expression profiles reveal distinct regulatory activities of transcription factors GATA1 and TAL1 upon AML relapse

The purpose of this study is to identify key regulatory pathways that potentially drive abnormal gene expression program in relapsed Acute Myeloid Leukemia (AML) patients, by integrative computational analyses on multi-omics molecular profiles. Relapsed AML remains a clinical challenge. Epigenetic heterogeneity may contribute to transcriptional dysregulation and disease progression in AML. However, what specific transcriptional programs and potential regulatory mechanisms contribute to disease relapse are not yet well understood. To characterize the global transcriptional landscapes in relapsed AML, we integrated genomics data from two cohorts of matched diagnosis and relapse patient specimens. We identified 5,416 differentially expressed genes (DEGs) between diagnosis and relapse in Cohort I. Unsupervised clustering yielded three distinct DEG groups: group A, B and C genes that were predominantly (88%) down-regulated, divergently regulated, or predominantly (65%) up-regulated, respectively, upon relapse. The expression pattern of all DEGs separated the patients into two clusters, most robustly by Group B genes. Interestingly, the majority of DEGs did not associate with changes in gene promoter methylation. Similar patterns were observed in Cohort II. We used Binding Analysis for Regulation of Transcription (BART) to identify transcriptional regulators (TRs) that potentially regulated the DEGs not associated with DNA methylation changes, and assessed the differential expression of identified TRs during disease progression. PU.1 was identified as a potential TR for Group A genes and was down-regulated upon relapse. GATA1 and TAL1 were identified as regulating Group B genes and were up-regulated in patient cluster1 and down-regulated in cluster2, consistent with the expression pattern of Group B genes. RBBP5 was a top predicted TR for Group C genes and was up-regulated upon relapse. We next validated the potential functionality of those predicted factors. In NSG mice transplanted with a human AML specimen, TAL1 and GATA1 were downregulated in AML cells collected four weeks after chemotherapy treatment, and were inferred as TRs for the down-regulated genes, similar to the patient data. PU.1 was inferred as regulating the up-regulated genes. Furthermore, we found that the level of differential expression of TAL1, GATA1, and PU.1 in each patient specimen associated with the correlation of DEG profiles between the patient specimen and TR perturbation in human-derived hematopoietic cell lines. Our results support the possibility that in some AML patients, TRs with roles in hematopoiesis and leukemia might contribute to disease relapse. Further mechanistic studies deciphering the molecular and phenotypic events facilitated by these TRs will yield significant insight into disease biology and possible therapeutic targeting approaches in relapsed AML.

Citation Format: Zhenjia Wang, Yaseswini Neelamraju, Cem Meydan, Nicholas Dunham, Jorge Gandara, Tak Lee, Subhash Prajapati, Franck Rapaport, Caroline Sheridan, Paul Zumbo, Michael Becker, Lars Bullinger, Martin Carroll, Richard D’Andrea, Richard Dillon, Ross Levine, Christopher E. Mason, Ari Melnick, Donna Neuberg, Stefan Bekiranov, Chongzhi Zang, Francine E. Garrett-Bakelman. Gene expression profiles reveal distinct regulatory activities of transcription factors GATA1 and TAL1 upon AML relapse [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3155.

Harnessing the power of sphingolipids: Prospects for acute myeloid leukemia

Acute myeloid leukemia (AML) is an aggressive, heterogenous malignancy characterized by clonal expansion of bone marrow-derived myeloid progenitor cells. While our current understanding of the molecular and genomic landscape of AML has evolved dramatically and opened avenues for molecularly targeted therapeutics to improve upon standard intensive induction chemotherapy, curative treatments are elusive, particularly in older patients. Responses to current AML treatments are transient and incomplete, necessitating the development of novel treatment strategies to improve outcomes. To this end, harnessing the power of bioactive sphingolipids to treat cancer shows great promise. Sphingolipids are involved in many hallmarks of cancer of paramount importance in AML. Leukemic blast survival is influenced by cellular levels of ceramide, a bona fide pro-death molecule, and its conversion to signaling molecules such as sphingosine-1-phosphate and glycosphingolipids. Preclinical studies demonstrate the efficacy of therapeutics that target dysregulated sphingolipid metabolism as well as their combinatorial synergy with clinically-relevant therapeutics. Thus, increased understanding of sphingolipid dysregulation may be exploited to improve AML patient care and outcomes. This review summarizes the current knowledge of dysregulated sphingolipid metabolism in AML, evaluates how pro-survival sphingolipids promote AML pathogenesis, and discusses the therapeutic potential of targeting these dysregulated sphingolipid pathways.

Blocking UBE2N abrogates oncogenic immune signaling in acute myeloid leukemia

Dysregulation of innate immune signaling pathways is implicated in various hematologic malignancies. However, these pathways have not been systematically examined in acute myeloid leukemia (AML). We report that AML hematopoietic stem and progenitor cells (HSPCs) exhibit a high frequency of dysregulated innate immune-related and inflammatory pathways, referred to as oncogenic immune signaling states. Through gene expression analyses and functional studies in human AML cell lines and patient-derived samples, we found that the ubiquitin-conjugating enzyme UBE2N is required for leukemic cell function in vitro and in vivo by maintaining oncogenic immune signaling states. It is known that the enzyme function of UBE2N can be inhibited by interfering with thioester formation between ubiquitin and the active site. We performed in silico structure-based and cellular-based screens and identified two related small-molecule inhibitors UC-764864/65 that targeted UBE2N at its active site. Using these small-molecule inhibitors as chemical probes, we further revealed the therapeutic efficacy of interfering with UBE2N function. This resulted in the blocking of ubiquitination of innate immune- and inflammatory-related substrates in human AML cell lines. Inhibition of UBE2N function disrupted oncogenic immune signaling by promoting cell death of leukemic HSPCs while sparing normal HSPCs in vitro. Moreover, baseline oncogenic immune signaling states in leukemic cells derived from discrete subsets of patients with AML exhibited a selective dependency on UBE2N function in vitro and in vivo. Our study reveals that interfering with UBE2N abrogates leukemic HSPC function and underscores the dependency of AML cells on UBE2N-dependent oncogenic immune signaling states.

Validation of CRISPR targeting for proliferation and cytarabine resistance control genes in the acute myeloid leukemia cell line MOLM-13

Acute myeloid leukemia patients with FMS-like tyrosine kinase 3–internal tandem duplications and mixed lineage leukemia–protein AF9 fusion proteins suffer from poor clinical outcomes. The MOLM-13 acute myeloid leukemia cell line harbors both of these abnormalities and is used in CRISPR experiments to identify disease drivers. However, experimental observations may be biased or inconclusive in the absence of experimentally validated positive control genes. We validated sgRNAs for knockdown of TP53 for cell proliferation and for DCK knockdown and CDA upregulation for cytarabine resistance control genes in MOLM-13 cells. We have provided a detailed CRISPR protocol applicable to both gene knockdown or activation experiments and downstream leukemic phenotype analyses. Inclusion of these controls in CRISPR experiments will enhance the capacity to identify novel myeloid leukemia drivers in MOLM-13 cells.

CRISPR-Cas9 knockdown and upregulation approaches were used for changing the expression of control genes TP53, DCK, and CDA in the acute myeloid leukemia cell line MOLM-13. Lentivirus transduced cells were selected (sorted for green fluorescent protein positive cells or selected by antibiotic treatment), and knockdown or upregulation of genes in the stable cells was confirmed by real-time quantitative PCR and western blot analyses. EdU incorporation assay was used to measure cell proliferation. Colorimetric proliferation/survival assay and flow cytometric cell counting were used to measure cell growth and survival. Results from the experimental and control groups were compared using Student's t-test.

DNA methylation landscapes of 1538 breast cancers reveal a replication-linked clock, epigenomic instability and cis-regulation

DNA methylation is aberrant in cancer, but the dynamics, regulatory role and clinical implications of such epigenetic changes are still poorly understood. Here, reduced representation bisulfite sequencing (RRBS) profiles of 1538 breast tumors and 244 normal breast tissues from the METABRIC cohort are reported, facilitating detailed analysis of DNA methylation within a rich context of genomic, transcriptional, and clinical data. Tumor methylation from immune and stromal signatures are deconvoluted leading to the discovery of a tumor replication-linked clock with genome-wide methylation loss in non-CpG island sites. Unexpectedly, methylation in most tumor CpG islands follows two replication-independent processes of gain (MG) or loss (ML) that we term epigenomic instability. Epigenomic instability is correlated with tumor grade and stage, TP53 mutations and poorer prognosis. After controlling for these global trans-acting trends, as well as for X-linked dosage compensation effects, cis-specific methylation and expression correlations are uncovered at hundreds of promoters and over a thousand distal elements. Some of these targeted known tumor suppressors and oncogenes. In conclusion, this study demonstrates that global epigenetic instability can erode cancer methylomes and expose them to localized methylation aberrations in-cis resulting in transcriptional changes seen in tumors.

Frequent somatic TET2 mutations in chronic NK-LGL leukemia with distinct patterns of cytopenias

Chronic natural killer large granular lymphocyte (NK-LGL) leukemia, also referred to as chronic lymphoproliferative disorder of NK cells, is a rare disorder defined by prolonged expansion of clonal NK cells. Similar prevalence of STAT3 mutations in chronic T-LGL and NK-LGL leukemia is suggestive of common pathogenesis. We undertook whole-genome sequencing to identify mutations unique to NK-LGL leukemia. The results were analyzed to develop a resequencing panel that was applied to 58 patients. Phosphatidylinositol 3-kinase pathway gene mutations (PIK3CD/PIK3AP1) and TNFAIP3 mutations were seen in 5% and 10% of patients, respectively. TET2 was exceptional in that mutations were present in 16 (28%) of 58 patient samples, with evidence that TET2 mutations can be dominant and exclusive to the NK compartment. Reduced-representation bisulfite sequencing revealed that methylation patterns were significantly altered in TET2 mutant samples. The promoter of TET2 and that of PTPRD, a negative regulator of STAT3, were found to be methylated in additional cohort samples, largely confined to the TET2 mutant group. Mutations in STAT3 were observed in 19 (33%) of 58 patient samples, 7 of which had concurrent TET2 mutations. Thrombocytopenia and resistance to immunosuppressive agents were uniquely observed in those patients with only TET2 mutation (Games-Howell post hoc test, P = .0074; Fisher's exact test, P = .00466). Patients with STAT3 mutation, inclusive of those with TET2 comutation, had lower hematocrit, hemoglobin, and absolute neutrophil count compared with STAT3 wild-type patients (Welch's t test, P ≤ .015). We present the discovery of TET2 mutations in chronic NK-LGL leukemia and evidence that it identifies a unique molecular subtype.

An Esrrb and Nanog Cell Fate Regulatory Module Controlled by Feed Forward Loop Interactions

Cell fate decisions during development are governed by multi-factorial regulatory mechanisms including chromatin remodeling, DNA methylation, binding of transcription factors to specific loci, RNA transcription and protein synthesis. However, the mechanisms by which such regulatory “dimensions” coordinate cell fate decisions are currently poorly understood. Here we quantified the multi-dimensional molecular changes that occur in mouse embryonic stem cells (mESCs) upon depletion of Estrogen related receptor beta (Esrrb), a key pluripotency regulator. Comparative analyses of expression changes subsequent to depletion of Esrrb or Nanog, indicated that a system of interlocked feed-forward loops involving both factors, plays a central part in regulating the timing of mESC fate decisions. Taken together, our meta-analyses support a hierarchical model in which pluripotency is maintained by an Oct4-Sox2 regulatory module, while the timing of differentiation is regulated by a Nanog-Esrrb module.

Chemotherapy Induces Senescence-Like Resilient Cells Capable of Initiating AML Recurrence

Patients with acute myeloid leukemia (AML) frequently relapse after chemotherapy, yet the mechanism by which AML reemerges is not fully understood. Herein, we show that primary AML cells enter a senescence-like phenotype following chemotherapy in vitro and in vivo. This is accompanied by induction of senescence/inflammatory and embryonic diapause transcriptional programs, with downregulation of MYC and leukemia stem cell genes. Single-cell RNA sequencing suggested depletion of leukemia stem cells in vitro and in vivo, and enrichment for subpopulations with distinct senescence-like cells. This senescence effect was transient and conferred superior colony-forming and engraftment potential. Entry into this senescence-like phenotype was dependent on ATR, and persistence of AML cells was severely impaired by ATR inhibitors. Altogether, we propose that AML relapse is facilitated by a senescence-like resilience phenotype that occurs regardless of their stem cell status. Upon recovery, these post-senescence AML cells give rise to relapsed AMLs with increased stem cell potential. SIGNIFICANCE: Despite entering complete remission after chemotherapy, relapse occurs in many patients with AML. Thus, there is an urgent need to understand the relapse mechanism in AML and the development of targeted treatments to improve outcome. Here, we identified a senescence-like resilience phenotype through which AML cells can survive and repopulate leukemia.This article is highlighted in the In This Issue feature, p. 1307.

Cell-free DNA (cfDNA) and Exosome Profiling from a Year-Long Human Spaceflight Reveals Circulating Biomarkers

Liquid biopsies based on cell-free DNA (cfDNA) or exosomes provide a noninvasive approach to monitor human health and disease but have not been utilized for astronauts. Here, we profile cfDNA characteristics, including fragment size, cellular deconvolution, and nucleosome positioning, in an astronaut during a year-long mission on the International Space Station, compared to his identical twin on Earth and healthy donors. We observed a significant increase in the proportion of cell-free mitochondrial DNA (cf-mtDNA) inflight, and analysis of post-flight exosomes in plasma revealed a 30-fold increase in circulating exosomes and patient-specific protein cargo (including brain-derived peptides) after the year-long mission. This longitudinal analysis of astronaut cfDNA during spaceflight and the exosome profiles highlights their utility for astronaut health monitoring, as well as cf-mtDNA levels as a potential biomarker for physiological stress or immune system responses related to microgravity, radiation exposure, and the other unique environmental conditions of spaceflight.

Temporal Telomere and DNA Damage Responses in the Space Radiation Environment

Telomeres, repetitive terminal features of chromosomes essential for maintaining genome integrity, shorten with cell division, lifestyle factors and stresses, and environmental exposures, and so they provide a robust biomarker of health, aging, and age-related diseases. We assessed telomere length dynamics (changes over time) in three unrelated astronauts before, during, and after 1-year or 6-month missions aboard the International Space Station (ISS). Similar to our results for National Aeronautics and Space Administration's (NASA's) One-Year Mission twin astronaut (Garrett-Bakelman et al., 2019), significantly longer telomeres were observed during spaceflight for two 6-month mission astronauts. Furthermore, telomere length shortened rapidly after return to Earth for all three crewmembers and, overall, telomere length tended to be shorter after spaceflight than before spaceflight. Consistent with chronic exposure to the space radiation environment, signatures of persistent DNA damage responses were also detected, including mitochondrial and oxidative stress, inflammation, and telomeric and chromosomal aberrations, which together provide potential mechanistic insight into spaceflight-specific telomere elongation.

Targeted detection and quantitation of histone modifications from 1,000 cells

Histone post-translational modifications (PTMs) create a powerful regulatory mechanism for maintaining chromosomal integrity in cells. Histone acetylation and methylation, the most widely studied histone PTMs, act in concert with chromatin-associated proteins to control access to genetic information during transcription. Alterations in cellular histone PTMs have been linked to disease states and have crucial biomarker and therapeutic potential. Traditional bottom-up mass spectrometry of histones requires large numbers of cells, typically one million or more. However, for some cell subtype-specific studies, it is difficult or impossible to obtain such large numbers of cells and quantification of rare histone PTMs is often unachievable. An established targeted LC-MS/MS method was used to quantify the abundance of histone PTMs from cell lines and primary human specimens. Sample preparation was modified by omitting nuclear isolation and reducing the rounds of histone derivatization to improve detection of histone peptides down to 1,000 cells. In the current study, we developed and validated a quantitative LC-MS/MS approach tailored for a targeted histone assay of 75 histone peptides with as few as 10,000 cells. Furthermore, we were able to detect and quantify 61 histone peptides from just 1,000 primary human stem cells. Detection of 37 histone peptides was possible from 1,000 acute myeloid leukemia patient cells. We anticipate that this revised method can be used in many applications where achieving large cell numbers is challenging, including rare human cell populations.

Abstract 651: Chemotherapy induces a transient senescence-like state capable of initiating AML recurrence

In acute myeloid leukemia (AML), chemotherapy is frequently followed by disease relapse, yet the mechanism by which AML reemerges is not fully understood. We hypothesized that chemotherapy-induced senescence (CIS) facilitates survival to genotoxic exposure, allowing AML cells to endure treatment in a transiently dormant state while retaining potential for leukemic repopulation.

Here we show that, following treatment with cytarabine (AraC), primary AML cells exhibit hallmark senescence features, including morphological changes, senescence-associated transcriptomic alterations, and senescence-associated-beta-galactosidase (SA-β-gal) activity. To quantify SA-β-gal activity, we used a fluorogenic substrate that enabled sensitive cytometric quantification of the senescence marker in viable cells. Mechanistically, we demonstrate that genotoxic stress following AraC treatment activated the ATR kinase, a key regulator of DNA damage response. Inhibition of ATR abrogated CIS in leukemia cells and induced cell death, indicating that CIS mediates stress survival.

We next analyzed the gene expression profiles of the minor residual fraction of chemotherapy-surviving cells ex vivo and in patients after chemotherapy. We found that leukemia stem cell (LSC) gene signatures were largely depleted in chemotherapy-surviving cells. However, we did find that the CIS signature was positively enriched. To recapitulate relapse of disease, we established an AML relapse model using primary human specimens engrafted in NSG mice. Treatment of AML-engrafted NSG mice with a physiologically relevant regimen of AraC reduced leukemic burden to a small residual leukemia fraction 8 days after therapeutic initiation. This reduction in leukemic burden was reversed to full blown leukemia present at day 29. Gene expression analysis of purified human AML cells at nadir (day 8) showed significant enrichment for senescence signatures. This enrichment was reduced by day 29 post-treatment, as expression was partially reversed to untreated levels. While we did not observe an enrichment of LSC signatures at nadir, we found an enrichment of LSC gene signatures in 60% of xenografted AML cases after overt recurrence of disease, suggesting partial reversal of the senescence phenotype and recovery of the LSC phenotype. We observed a similar pattern, when we performed RNA-seq analysis in matched AML specimens from patients at diagnosis vs relapse. To examine if LSCs might be present in the CIS cell population, we performed scRNA-seq in primary AML before and after chemotherapy. Although chemotherapy killed a majority of AML cells, we did not find an enrichment of LSCs on single cell level compared to untreated control cells in most AML specimens. Instead, we found increased expression of the CIS gene signature. Next, to test if CIS cells possess self-renewal capability, we sorted for SA-β-gal+ cells after AraC exposure and seeded them in CFC assays. We found colony-forming potential of CIS cells, suggesting a capacity to repopulate leukemia. To confirm that leukemia can be initiated from CIS AML in vivo, we transplanted sorted CIS AML cells into NSG mice and found that these cells successfully repopulated leukemia.

Altogether, we propose a mechanism of AML relapse whereby AML cells survive chemotherapy via acquisition of a transient senescent-like state.

Citation Format: Meng Li, Matt Teater, Cem Meydan, Francine E. Garrett-Bakelman, Christopher Chin, Martin Carroll, Ari Melnick, Cihangir Duy. Chemotherapy induces a transient senescence-like state capable of initiating AML recurrence [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 651.

CD97 is a critical regulator of acute myeloid leukemia stem cell function

Despite significant efforts to improve therapies for acute myeloid leukemia (AML), clinical outcomes remain poor. Understanding the mechanisms that regulate the development and maintenance of leukemic stem cells (LSCs) is important to reveal new therapeutic opportunities. We have identified CD97, a member of the adhesion class of G protein-coupled receptors (GPCRs), as a frequently up-regulated antigen on AML blasts that is a critical regulator of blast function. High levels of CD97 correlate with poor prognosis, and silencing of CD97 reduces disease aggressiveness in vivo. These phenotypes are due to CD97's ability to promote proliferation, survival, and the maintenance of the undifferentiated state in leukemic blasts. Collectively, our data credential CD97 as a promising therapeutic target on LSCs in AML.

Rational Targeting of Cooperating Layers of the Epigenome Yields Enhanced Therapeutic Efficacy against AML

Disruption of epigenetic regulation is a hallmark of acute myeloid leukemia (AML), but epigenetic therapy is complicated by the complexity of the epigenome. Herein, we developed a long-term primary AML ex vivo platform to determine whether targeting different epigenetic layers with 5-azacytidine and LSD1 inhibitors would yield improved efficacy. This combination was most effective in TET2 ^mut AML, where it extinguished leukemia stem cells and particularly induced genes with both LSD1-bound enhancers and cytosine-methylated promoters. Functional studies indicated that derepression of genes such as GATA2 contributes to drug efficacy. Mechanistically, combination therapy increased enhancer-promoter looping and chromatin-activating marks at the GATA2 locus. CRISPRi of the LSD1-bound enhancer in patient-derived TET2 ^mut AML was associated with dampening of therapeutic GATA2 induction. TET2 knockdown in human hematopoietic stem/progenitor cells induced loss of enhancer 5-hydroxymethylation and facilitated LSD1-mediated enhancer inactivation. Our data provide a basis for rational targeting of cooperating aberrant promoter and enhancer epigenetic marks driven by mutant epigenetic modifiers. SIGNIFICANCE: Somatic mutations of genes encoding epigenetic modifiers are a hallmark of AML and potentially disrupt many components of the epigenome. Our study targets two different epigenetic layers at promoters and enhancers that cooperate to aberrant gene silencing, downstream of the actions of a mutant epigenetic regulator.This article is highlighted in the In This Issue feature, p. 813.

The NASA Twins Study: A multidimensional analysis of a year-long human spaceflight

To understand the health impact of long-duration spaceflight, one identical twin astronaut was monitored before, during, and after a 1-year mission onboard the International Space Station; his twin served as a genetically matched ground control. Longitudinal assessments identified spaceflight-specific changes, including decreased body mass, telomere elongation, genome instability, carotid artery distension and increased intima-media thickness, altered ocular structure, transcriptional and metabolic changes, DNA methylation changes in immune and oxidative stress-related pathways, gastrointestinal microbiota alterations, and some cognitive decline postflight. Although average telomere length, global gene expression, and microbiome changes returned to near preflight levels within 6 months after return to Earth, increased numbers of short telomeres were observed and expression of some genes was still disrupted. These multiomic, molecular, physiological, and behavioral datasets provide a valuable roadmap of the putative health risks for future human spaceflight.

CD99 is a therapeutic target on disease stem cells in myeloid malignancies

Acute myeloid leukemia (AML) and the myelodysplastic syndromes (MDS) are initiated and sustained by self-renewing malignant stem cells; thus, eradication of AML and MDS stem cells is required for cure. We identified CD99 as a cell surface protein frequently overexpressed on AML and MDS stem cells. Expression of CD99 allows for prospective separation of leukemic stem cells (LSCs) from functionally normal hematopoietic stem cells in AML, and high CD99 expression on AML blasts enriches for functional LSCs as demonstrated by limiting dilution xenotransplant studies. Monoclonal antibodies (mAbs) targeting CD99 induce the death of AML and MDS cells in a SARC family kinase-dependent manner in the absence of immune effector cells or complement, and anti-CD99 mAbs exhibit antileukemic activity in AML xenografts. These data establish CD99 as a marker of AML and MDS stem cells, as well as a promising therapeutic target in these disorders.

Cooperative Epigenetic Remodeling by TET2 Loss and NRAS Mutation Drives Myeloid Transformation and MEK Inhibitor Sensitivity

Mutations in epigenetic modifiers and signaling factors often co-occur in myeloid malignancies, including TET2 and NRAS mutations. Concurrent Tet2 loss and Nras^G12D expression in hematopoietic cells induced myeloid transformation, with a fully penetrant, lethal chronic myelomonocytic leukemia (CMML), which was serially transplantable. Tet2 loss and Nras mutation cooperatively led to decrease in negative regulators of mitogen-activated protein kinase (MAPK) activation, including Spry2, thereby causing synergistic activation of MAPK signaling by epigenetic silencing. Tet2/Nras double-mutant leukemia showed preferential sensitivity to MAPK kinase (MEK) inhibition in both mouse model and patient samples. These data provide insights into how epigenetic and signaling mutations cooperate in myeloid transformation and provide a rationale for mechanism-based therapy in CMML patients with these high-risk genetic lesions.

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 Idh2^R140Q or Dnmt3a^R882H single-mutant and Idh2^R140Q/Dnmt3a^R882H 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.

Combination Targeted Therapy to Disrupt Aberrant Oncogenic Signaling and Reverse Epigenetic Dysfunction in IDH2- and TET2-Mutant Acute Myeloid Leukemia

Genomic studies in acute myeloid leukemias (AML) have identified mutations that drive altered DNA methylation, including TET2 and IDH2 Here, we show that models of AML resulting from TET2 or IDH2 mutations combined with FLT3^ITD mutations are sensitive to 5-azacytidine or to the IDH2 inhibitor AG-221, respectively. 5-azacytidine and AG-221 treatment induced an attenuation of aberrant DNA methylation and transcriptional output and resulted in a reduction in leukemic blasts consistent with antileukemic activity. These therapeutic benefits were associated with restoration of leukemic cell differentiation, and the normalization of hematopoiesis was derived from mutant cells. By contrast, combining AG-221 or 5-azacytidine with FLT3 inhibition resulted in a reduction in mutant allele burden, progressive recovery of normal hematopoiesis from non-mutant stem-progenitor cells, and reversal of dysregulated DNA methylation and transcriptional output. Together, our studies suggest combined targeting of signaling and epigenetic pathways can increase therapeutic response in AML.Significance: AMLs with mutations in TET2 or IDH2 are sensitive to epigenetic therapy through inhibition of DNA methyltransferase activity by 5-azacytidine or inhibition of mutant IDH2 through AG-221. These inhibitors induce a differentiation response and can be used to inform mechanism-based combination therapy. Cancer Discov; 7(5); 494-505. ©2017 AACR.See related commentary by Thomas and Majeti, p. 459See related article by Yen et al., p. 478This article is highlighted in the In This Issue feature, p. 443.

Mutant IDH: a targetable driver of leukemic phenotypes linking metabolism, epigenetics and transcriptional regulation

Aberrant epigenomic programming is a hallmark of acute myeloid leukemia. This is partially due to somatic mutations that perturb cytosine methylation, histone post-translational modifications and transcription factors. Remarkably, mutations in the IDH1 and IDH2 genes perturb the epigenome through all three of these mechanisms. Mutant IDH enzymes produce high levels of the oncometabolite (R)-2-hydroxyglutarate that competitively inhibits dioxygenase enzymes that modify methylcytosine to hydroxymethylcytosine and histone tail methylation. The development of IDH mutant specific inhibitors may now enable the therapeutic reprogramming of both layers of the epigenome spontaneously to revert the malignant phenotype of these leukemias and improve clinical outcome for acute myeloid leukemia patients with IDH mutations.

Abstract LB-073: Epigenome evolution in relapsed acute myeloid leukemia

Acute myeloid leukemia (AML) is a predominantly fatal hematopoietic malignancy with high inter-patient and intra-patient genetic and epigenetic heterogeneity. The prognosis of relapsed AML remains dismal, yet the epigenetic basis of relapse is still unclear. Here we investigated whether and how the epigenome evolution impacts AML progression with biological and clinical relevance.

Methods: We obtained clinical annotation and AML specimens from 138 patients with paired diagnosis and relapsed samples. We used normal bone marrow (NBM) as epigenetic/transcriptomic controls and patients’ matched germline DNA as genetic controls. We then performed DNA methylation sequencing (ERRBS), RNA-seq, and Exome-seq. For one patient with 5 serial time points, we performed whole genome sequencing (WGS), ERRBS, and single cell RNA-seq. We measured the epigenetic allele burden using a compositional entropy-based approach (Methclone) and methylation heterogeneity using epipolymorphism.

Results: We found that diagnosis stage epigenetic allele burden (ΔS < -90) was linked to an inferior clinical outcome (p = 0.0064, log-rank test of relapse-free survival). The higher significance in promoter regions implies the functional impact of epigenetic dynamics. Promoter epiallele shift was associated with more differential expression events (p = 3.8 × 10−6, Wilcoxon signed-rank test) and promoter epiallele diversity is significantly associated with single cell resolution transcriptional heterogeneity (p < 2.2 × 10−16, ANOVA test). The global methylation heterogeneity is decreased from diagnosis to relapse, indicating a selective impact of chemotherapy on epigenetic variability (p = 0.0056, paired Wilcoxon test).

We investigated epigenetic allele burden progression from diagnosis to relapse by classifying patients into three clusters using K-means clustering: those with 1) decreased, 2) stable, or 3) increased abundance of epiallele burden. No association was seen between epigenetic clusters and patterns of genetic evolution, and the genetic abundance is higher in Cluster 3 than Cluster 1 (p = 0.048, Wilcoxon test), indicating divergent paths of genetic and epigenetic evolution. We next examined differential expression in the epigenetic cluster samples at diagnosis compared to NBM. Cluster 1 specific genes were enriched for cell cycle processes, while Cluster 3 genes were enriched for immune responses (p < 0.001, gene ontology hypergeometric tests). Integrating WGS and ERRBS data showed that epiallele burden is more dynamic than somatic mutations; a significant increase in epiallele burden preceded a major increase of somatic mutational abundance.

Summary: Our results indicate that epigenetic dynamics may provide leukemia cells greater evolutionary fitness via transcriptional adaptation and is associated with clinical outcome. This provides an alternative mechanism of AML resilience during progression and a potential predictor of relapse.

Citation Format: Sheng Li, Francine E. Garrett-Bakelman, Stephen S. Chung, Todd Hricik, Franck Rapaport, Jay Patel, Richard Dillon, Priyanka Vijay, Anna L. Brown, Alexander E. Perl, Joy Connon, Mathijs A. Sanders, Peter J.M. Valk, Lars Bullinger, Selina Luger, Michael W. Becker, Ian D. Lewis, Luen Bik To, Richard J. D’Andrea, David Grimwade, Ruud Delwel, Bob Löwenberg, Hartmut Döhner, Konstanze Döhner, Monica L. Guzman, Duane C. Hassane, Gail J. Roboz, Martin Carroll, Christopher Y. Park, Donna S. Neuberg, Ross L. Levine, Ari M. Melnick, Christopher E. Mason. Epigenome evolution in relapsed acute myeloid leukemia. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-073.

Distinct evolution and dynamics of epigenetic and genetic heterogeneity in acute myeloid leukemia

Genetic heterogeneity contributes to clinical outcome and progression of most tumors. Yet, little is known regarding allelic diversity for epigenetic compartments and almost no data exists for acute myeloid leukemia (AML). Here we examined epigenetic heterogeneity as assessed by cytosine methylation within defined genomic loci with four CpGs (epigenetic alleles), somatic mutations and transcriptomes of AML patient samples at serial time points. We observe that epigenetic allele burden is linked to inferior outcome and varies considerably during disease progression. Epigenetic and genetic allelic burden and patterning follow different patterns and kinetics during disease progression. We observed a subset of AMLs with high epiallele and low somatic mutation burden at diagnosis, a subset with high somatic mutation and lower epiallele burdens at diagnosis, and a subset with a mixed profile, suggesting distinct modes of tumor heterogeneity. Genes linked to promoter-associated epiallele shifts during tumor progression display increased single-cell transcriptional variance and differential expression, suggesting functional impact on gene regulation. Thus, genetic and epigenetic heterogeneity can occur with distinct kinetics, each likely able to impact biological and clinical features of tumors.

Dnmt3a regulates myeloproliferation and liver-specific expansion of hematopoietic stem and progenitor cells

DNMT3A mutations are observed in myeloid malignancies, including myeloproliferative neoplasms (MPN), myelodysplastic syndromes (MDS), and acute myeloid leukemia (AML). Transplantation studies have elucidated an important role for Dnmt3a in stem cell self-renewal and in myeloid differentiation. Here we investigated the impact of conditional hematopoietic Dnmt3a loss on disease phenotype in primary mice. Mx1-Cre-mediated Dnmt3a ablation led to the development of a lethal, fully penetrant myeloproliferative neoplasm with myelodysplasia (MDS/MPN) characterized by peripheral cytopenias and by marked, progressive hepatomegaly. We detected expanded stem/progenitor populations in the liver of Dnmt3a-ablated mice. The MDS/MPN induced by Dnmt3a ablation was transplantable, including the marked hepatomegaly. Homing studies showed that Dnmt3a-deleted bone marrow cells preferentially migrated to the liver. Gene expression and DNA methylation analyses of progenitor cell populations identified differential regulation of hematopoietic regulatory pathways, including fetal liver hematopoiesis transcriptional programs. These data demonstrate that Dnmt3a ablation in the hematopoietic system leads to myeloid transformation in vivo, with cell autonomous aberrant tissue tropism and marked extramedullary hematopoiesis (EMH) with liver involvement. Hence, in addition to the established role of Dnmt3a in regulating self-renewal, Dnmt3a regulates tissue tropism and limits myeloid progenitor expansion in vivo.

Enhanced reduced representation bisulfite sequencing for assessment of DNA methylation at base pair resolution

DNA methylation pattern mapping is heavily studied in normal and diseased tissues. A variety of methods have been established to interrogate the cytosine methylation patterns in cells. Reduced representation of whole genome bisulfite sequencing was developed to detect quantitative base pair resolution cytosine methylation patterns at GC-rich genomic loci. This is accomplished by combining the use of a restriction enzyme followed by bisulfite conversion. Enhanced Reduced Representation Bisulfite Sequencing (ERRBS) increases the biologically relevant genomic loci covered and has been used to profile cytosine methylation in DNA from human, mouse and other organisms. ERRBS initiates with restriction enzyme digestion of DNA to generate low molecular weight fragments for use in library preparation. These fragments are subjected to standard library construction for next generation sequencing. Bisulfite conversion of unmethylated cytosines prior to the final amplification step allows for quantitative base resolution of cytosine methylation levels in covered genomic loci. The protocol can be completed within four days. Despite low complexity in the first three bases sequenced, ERRBS libraries yield high quality data when using a designated sequencing control lane. Mapping and bioinformatics analysis is then performed and yields data that can be easily integrated with a variety of genome-wide platforms. ERRBS can utilize small input material quantities making it feasible to process human clinical samples and applicable in a range of research applications. The video produced demonstrates critical steps of the ERRBS protocol.

An optimized algorithm for detecting and annotating regional differential methylation

BACKGROUND

DNA methylation profiling reveals important differentially methylated regions (DMRs) of the genome that are altered during development or that are perturbed by disease. To date, few programs exist for regional analysis of enriched or whole-genome bisulfate conversion sequencing data, even though such data are increasingly common. Here, we describe an open-source, optimized method for determining empirically based DMRs (eDMR) from high-throughput sequence data that is applicable to enriched whole-genome methylation profiling datasets, as well as other globally enriched epigenetic modification data.

RESULTS

Here we show that our bimodal distribution model and weighted cost function for optimized regional methylation analysis provides accurate boundaries of regions harboring significant epigenetic modifications. Our algorithm takes the spatial distribution of CpGs into account for the enrichment assay, allowing for optimization of the definition of empirical regions for differential methylation. Combined with the dependent adjustment for regional p-value combination and DMR annotation, we provide a method that may be applied to a variety of datasets for rapid DMR analysis. Our method classifies both the directionality of DMRs and their genome-wide distribution, and we have observed that shows clinical relevance through correct stratification of two Acute Myeloid Leukemia (AML) tumor sub-types.

CONCLUSIONS

Our weighted optimization algorithm eDMR for calling DMRs extends an established DMR R pipeline (methylKit) and provides a needed resource in epigenomics. Our method enables an accurate and scalable way of finding DMRs in high-throughput methylation sequencing experiments. eDMR is available for download at http://code.google.com/p/edmr/.

methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles

DNA methylation is a chemical modification of cytosine bases that is pivotal for gene regulation, cellular specification and cancer development. Here, we describe an R package, methylKit, that rapidly analyzes genome-wide cytosine epigenetic profiles from high-throughput methylation and hydroxymethylation sequencing experiments. methylKit includes functions for clustering, sample quality visualization, differential methylation analysis and annotation features, thus automating and simplifying many of the steps for discerning statistically significant bases or regions of DNA methylation. Finally, we demonstrate methylKit on breast cancer data, in which we find statistically significant regions of differential methylation and stratify tumor subtypes. methylKit is available at http://code.google.com/p/methylkit.

Base-pair resolution DNA methylation sequencing reveals profoundly divergent epigenetic landscapes in acute myeloid leukemia

We have developed an enhanced form of reduced representation bisulfite sequencing with extended genomic coverage, which resulted in greater capture of DNA methylation information of regions lying outside of traditional CpG islands. Applying this method to primary human bone marrow specimens from patients with Acute Myelogeneous Leukemia (AML), we demonstrated that genetically distinct AML subtypes display diametrically opposed DNA methylation patterns. As compared to normal controls, we observed widespread hypermethylation in IDH mutant AMLs, preferentially targeting promoter regions and CpG islands neighboring the transcription start sites of genes. In contrast, AMLs harboring translocations affecting the MLL gene displayed extensive loss of methylation of an almost mutually exclusive set of CpGs, which instead affected introns and distal intergenic CpG islands and shores. When analyzed in conjunction with gene expression profiles, it became apparent that these specific patterns of DNA methylation result in differing roles in gene expression regulation. However, despite this subtype-specific DNA methylation patterning, a much smaller set of CpG sites are consistently affected in both AML subtypes. Most CpG sites in this common core of aberrantly methylated CpGs were hypermethylated in both AML subtypes. Therefore, aberrant DNA methylation patterns in AML do not occur in a stereotypical manner but rather are highly specific and associated with specific driving genetic lesions.

Acute myeloid leukemias (AML) are a group of malignancies that originate in the bone marrow. While many different genetic lesions have been linked to the different forms of this disease, it is also clear that these genetic lesions are not always sufficient to cause AML. DNA methylation plays a role in gene expression regulation, and abnormal distribution of DNA methylation has been observed in many cancers, including AML. Here we demonstrate that changes in DNA methylation in AML are not uniform across all AML subtypes, but rather they display unique patterns, which are closely linked to the underlying genetic lesions of each of the different forms of AML. Furthermore, these unique patterns of DNA methylation have different impacts on gene expression regulation in each AML subtype.