Erythroid Predominance in Bone Marrow Biopsies of Acute Myeloid Leukemia Patients after Decitabine Treatment Correlates with Mutation Profile and Complete Remission.

INTRODUCTION: Acute myeloid leukemia (AML) patients may receive hypomethylating agents such as decitabine (DAC) as part of their treatment. Not all patients respond to this therapy, and if they do, the clinical response may occur only after 3-6 courses of treatment. Hence, early biomarkers predicting response would be very useful. METHODS: We retrospectively analyzed a cohort of 22 AML patients who were treated with DAC. Histology of the bone marrow biopsy, pathogenic mutations, and methylation status were related to the treatment response. RESULTS: In 8/22 (36%) patients, an erythroid dominant response (EDR) pattern, defined as a ratio of myeloid cells/erythroid cells <1, was observed. In the remaining 14 cases, a myeloid predominance was preserved during treatment. No difference in the hypomethylating effect of DAC treatment was observed in patients with and without EDR, as global 5-methylcytosine levels dropped similarly in both groups. Mutational analysis by NGS using a panel of commonly mutated genes in AML showed that patients with an early EDR harbored on average less mutations, with U2AF1 mutations occurring more frequently, whereas RUNX1 mutations were underrepresented compared to non-EDR cases. Interestingly, the development of an EDR correlated with complete remission (7/8 cases with an EDR vs. only 2/14 cases without an EDR). CONCLUSION: We conclude that early histological bone marrow examination for the development of an EDR may be helpful to predict response in AML patients during treatment with DAC.

GFI1B and LSD1 repress myeloid traits during megakaryocyte differentiation.

The transcription factor Growth Factor Independence 1B (GFI1B) recruits Lysine Specific Demethylase 1 A (LSD1/KDM1A) to stimulate gene programs relevant for megakaryocyte and platelet biology. Inherited pathogenic GFI1B variants result in thrombocytopenia and bleeding propensities with varying intensity. Whether these affect similar gene programs is unknow. Here we studied transcriptomic effects of four patient-derived GFI1B variants (GFI1BT174N,H181Y,R184P,Q287*) in MEG01 megakaryoblasts. Compared to normal GFI1B, each variant affected different gene programs with GFI1BQ287* uniquely failing to repress myeloid traits. In line with this, single cell RNA-sequencing of induced pluripotent stem cell (iPSC)-derived megakaryocytes revealed a 4.5-fold decrease in the megakaryocyte/myeloid cell ratio in GFI1BQ287* versus normal conditions. Inhibiting the GFI1B-LSD1 interaction with small molecule GSK-LSD1 resulted in activation of myeloid genes in normal iPSC-derived megakaryocytes similar to what was observed for GFI1BQ287* iPSC-derived megakaryocytes. Thus, GFI1B and LSD1 facilitate gene programs relevant for megakaryopoiesis while simultaneously repressing programs that induce myeloid differentiation.

Susceptibility of pediatric acute lymphoblastic leukemia to STAT3 inhibition depends on p53 induction.

Advances in the clinical management of pediatric B-cell acute lymphoblastic leukemia (B-ALL) have dramatically improved outcomes for this disease. However, relapsed and high-risk disease still contribute to significant numbers of treatment failures. Development of new, broad range therapies is urgently needed for these cases. We previously reported the susceptibility of ETV6-RUNX1+ pediatric B-ALL to inhibition of signal transducer and activator of transcription 3 (STAT3) activity. In the present study, we demonstrate that pharmacological or genetic inhibition of STAT3 results in p53 induction and that CRISPR-mediated TP53 knockout substantially reverses susceptibility to STAT3 inhibition. Furthermore, we demonstrate that sensitivity to STAT3 inhibition in patient-derived xenograft (PDX) B-ALL samples is not restricted to any particular disease subtype, but rather depends on TP53 status, the only resistant samples being TP53 mutant. Induction of p53 following STAT3 inhibition is not directly dependent on MDM2 but correlates with degradation of MDM4. As such, STAT3 inhibition exhibits synergistic in vitro and in vivo anti-leukemia activity when combined with MDM2 inhibition. Taken together with the relatively low frequency of TP53 mutations in this disease, these data support the future development of combined STAT3/ MDM2 inhibition in the therapy of refractory and relapsed pediatric B-ALL.

scANANSE gene regulatory network and motif analysis of single-cell clusters.

The recent development of single-cell techniques is essential to unravel complex biological systems. By measuring the transcriptome and the accessible genome on a single-cell level, cellular heterogeneity in a biological environment can be deciphered. Transcription factors act as key regulators activating and repressing downstream target genes, and together they constitute gene regulatory networks that govern cell morphology and identity. Dissecting these gene regulatory networks is crucial for understanding molecular mechanisms and disease, especially within highly complex biological systems. The gene regulatory network analysis software ANANSE and the motif enrichment software GimmeMotifs were both developed to analyse bulk datasets. We developed scANANSE, a software pipeline for gene regulatory network analysis and motif enrichment using single-cell RNA and ATAC datasets. The scANANSE pipeline can be run from either R or Python. First, it exports data from standard single-cell objects. Next, it automatically runs multiple comparisons of cell cluster data. Finally, it imports the results back to the single-cell object, where the result can be further visualised, integrated, and interpreted. Here, we demonstrate our scANANSE pipeline on a publicly available PBMC multi-omics dataset. It identifies well-known cell type-specific hematopoietic factors. Importantly, we also demonstrated that scANANSE combined with GimmeMotifs is able to predict transcription factors with both activating and repressing roles in gene regulation.

Acid ceramidase regulates innate immune memory.

Innate immune memory, also called "trained immunity," is a functional state of myeloid cells enabling enhanced immune responses. This phenomenon is important for host defense, but also plays a role in various immune-mediated conditions. We show that exogenously administered sphingolipids and inhibition of sphingolipid metabolizing enzymes modulate trained immunity. In particular, we reveal that acid ceramidase, an enzyme that converts ceramide to sphingosine, is a potent regulator of trained immunity. We show that acid ceramidase regulates the transcription of histone-modifying enzymes, resulting in profound changes in histone 3 lysine 27 acetylation and histone 3 lysine 4 trimethylation. We confirm our findings by identifying single-nucleotide polymorphisms in the region of ASAH1, the gene encoding acid ceramidase, that are associated with the trained immunity cytokine response. Our findings reveal an immunomodulatory effect of sphingolipids and identify acid ceramidase as a relevant therapeutic target to modulate trained immunity responses in innate immune-driven disorders.

Understanding blood development and leukemia using sequencing-based technologies and human cell systems.

Our current understanding of human hematopoiesis has undergone significant transformation throughout the years, challenging conventional views. The evolution of high-throughput technologies has enabled the accumulation of diverse data types, offering new avenues for investigating key regulatory processes in blood cell production and disease. In this review, we will explore the opportunities presented by these advancements for unraveling the molecular mechanisms underlying normal and abnormal hematopoiesis. Specifically, we will focus on the importance of enhancer-associated regulatory networks and highlight the crucial role of enhancer-derived transcription regulation. Additionally, we will discuss the unprecedented power of single-cell methods and the progression in using in vitro human blood differentiation system, in particular induced pluripotent stem cell models, in dissecting hematopoietic processes. Furthermore, we will explore the potential of ever more nuanced patient profiling to allow precision medicine approaches. Ultimately, we advocate for a multiparameter, regulatory network-based approach for providing a more holistic understanding of normal hematopoiesis and blood disorders.

Saponin-based adjuvants enhance antigen cross-presentation in human CD11c+ CD1c+ CD5- CD163+ conventional type 2 dendritic cells.

BACKGROUND: Adjuvants are key for effective vaccination against cancer and chronic infectious diseases. Saponin-based adjuvants (SBAs) are unique among adjuvants in their ability to induce robust cell-mediated immune responses in addition to antibody responses. Recent preclinical studies revealed that SBAs induced cross-presentation and lipid bodies in otherwise poorly cross-presenting CD11b+ murine dendritic cells (DCs). METHOD: Here, we investigated the response of human DC subsets to SBAs with RNA sequencing and pathway analyses, lipid body induction visualized by laser scanning microscopy, antigen translocation to the cytosol, and antigen cross-presentation to CD8+ T cells. RESULTS: RNA sequencing of SBA-treated conventional type 1 DC (cDC1) and type 2 DC (cDC2) subsets uncovered that SBAs upregulated lipid-related pathways in CD11c+ CD1c+ cDC2s, especially in the CD5- CD163+ CD14+ cDC2 subset. Moreover, SBAs induced lipid bodies and enhanced endosomal antigen translocation into the cytosol in this particular cDC2 subset. Finally, SBAs enhanced cross-presentation only in cDC2s, which requires the CD163+ CD14+ cDC2 subset. CONCLUSIONS: These data thus identify the CD163+ CD14+ cDC2 subset as the main SBA-responsive DC subset in humans and imply new strategies to optimize the application of saponin-based adjuvants in a potent cancer vaccine.

Inducible MLL-AF9 Expression Drives an AML Program during Human Pluripotent Stem Cell-Derived Hematopoietic Differentiation.

A t(9;11)(p22;q23) translocation produces the MLL-AF9 fusion protein, which is found in up to 25% of de novo AML cases in children. Despite major advances, obtaining a comprehensive understanding of context-dependent MLL-AF9-mediated gene programs during early hematopoiesis is challenging. Here, we generated a human inducible pluripotent stem cell (hiPSC) model with a doxycycline dose-dependent MLL-AF9 expression. We exploited MLL-AF9 expression as an oncogenic hit to uncover epigenetic and transcriptomic effects on iPSC-derived hematopoietic development and the transformation into (pre-)leukemic states. In doing so, we observed a disruption in early myelomonocytic development. Accordingly, we identified gene profiles that were consistent with primary MLL-AF9 AML and uncovered high-confidence MLL-AF9-associated core genes that are faithfully represented in primary MLL-AF9 AML, including known and presently unknown factors. Using single-cell RNA-sequencing, we identified an increase of CD34 expressing early hematopoietic progenitor-like cell states as well as granulocyte-monocyte progenitor-like cells upon MLL-AF9 activation. Our system allows for careful chemically controlled and stepwise in vitro hiPSC-derived differentiation under serum-free and feeder-free conditions. For a disease that currently lacks effective precision medicine, our system provides a novel entry-point into exploring potential novel targets for personalized therapeutic strategies.

Differences in Immune Responses in Individuals of Indian and European Origin: Relevance for the COVID-19 Pandemic.

During the coronavirus disease 2019 (COVID-19) pandemic, large differences in susceptibility and mortality due to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection have been reported between populations in Europe and South Asia. While both host and environmental factors (including Mycobacterium bovis BCG vaccination) have been proposed to explain this, the potential biological substrate of these differences is unknown. We purified peripheral blood mononuclear cells from individuals living in India and the Netherlands at baseline and 10 to 12 weeks after BCG vaccination. We compared chromatin accessibility between the two populations at baseline, as well as gene transcription profiles and cytokine production capacities upon stimulation. The chromatin accessibility of genes important for adaptive immunity was higher in the Indians than in the Europeans, while the latter had more accessible chromatin regions in genes of the innate immune system. At the transcriptional level, we observed that the Indian volunteers displayed a more tolerant immune response to stimulation, in contrast to a more exaggerated response in the Europeans. BCG vaccination strengthened the tolerance program in the Indians but not in the Europeans. These differences may partly explain the different impact of COVID-19 on the two populations. IMPORTANCE In this study, we assessed the differences in immune responses in individuals from India and Europe. This aspect is of great relevance, because of the described differences in morbidity and mortality between India and Europe during the pandemic. We found a significant difference in chromatin accessibility in immune cells from the two populations, followed by a more balanced and effective response in individuals from India. These exciting findings represent a very important piece of the puzzle for understanding the COVID-19 pandemic at a global level.

Transcription factor HNF1ß controls a transcriptional network regulating kidney cell structure and tight junction integrity.

Mutations in the hepatocyte nuclear factor (HNF)1ß gene (HNF1B) cause autosomal dominant tubulointerstitial kidney disease, a rare and heterogeneous disease characterized by renal cysts and/or malformation, maturity-onset diabetes of the young, hypomagnesemia, and hypokalemia. The electrolyte disturbances may develop in the distal part of the nephron, which is important for fine-tuning of Mg2+ and Ca2+ reabsorption. Therefore, we aimed to study the transcriptional network directed by HNF1ß in the distal part of the nephron. We combined HNF1ß chromatin immunoprecipitation-sequencing and mRNA expression data to identify direct targets of HNF1ß in a renal distal convoluted tubule cell line (mpkDCT). Gene Ontology term pathway analysis demonstrated enrichment of cell polarity, cell-cell junction, and cytoskeleton pathways in the dataset. Genes directly and indirectly regulated by HNF1ß within these pathways included members of the apical and basolateral polarity complexes including Crumbs protein homolog 3 (Crb3), partitioning defective 6 homolog-ß (Pard6b), and LLGL Scribble cell polarity complex component 2 (Llgl2). In monolayers of mouse inner medullary collecting duct 3 cells expressing dominant negative Hnf1b, tight junction integrity was compromised, as observed by reduced transepithelial electrical resistance values and increased permeability for fluorescein (0.4 kDa) compared with wild-type cells. Expression of dominant negative Hnf1b also led to a decrease in height (30%) and an increase in surface (58.5%) of cells grown on membranes. Moreover, three-dimensional spheroids formed by cells expressing dominant negative Hnf1b were reduced in size compared with wild-type spheroids (30%). Together, these findings demonstrate that HNF1ß directs a transcriptional network regulating tight junction integrity and cell structure in the distal part of the nephron.NEW & NOTEWORTHY Genetic defects in transcription factor hepatocyte nuclear factor (HNF)1ß cause a heterogeneous disease characterized by electrolyte disturbances, kidney cysts, and diabetes. By combining RNA-sequencing and HNF1ß chromatin immunoprecipitation-sequencing data, we identified new HNF1ß targets that were enriched for cell polarity pathways. Newly discovered targets included members of polarity complexes Crb3, Pard6b, and Llgl2. Functional assays in kidney epithelial cells demonstrated decreased tight junction integrity and a loss of typical cuboidal morphology in mutant Hnf1b cells.

Regulome analysis in B-acute lymphoblastic leukemia exposes Core Binding Factor addiction as a therapeutic vulnerability.

The ETV6-RUNX1 onco-fusion arises in utero, initiating a clinically silent pre-leukemic state associated with the development of pediatric B-acute lymphoblastic leukemia (B-ALL). We characterize the ETV6-RUNX1 regulome by integrating chromatin immunoprecipitation- and RNA-sequencing and show that ETV6-RUNX1 functions primarily through competition for RUNX1 binding sites and transcriptional repression. In pre-leukemia, this results in ETV6-RUNX1 antagonization of cell cycle regulation by RUNX1 as evidenced by mass cytometry analysis of B-lineage cells derived from ETV6-RUNX1 knock-in human pluripotent stem cells. In frank leukemia, knockdown of RUNX1 or its co-factor CBFß results in cell death suggesting sustained requirement for RUNX1 activity which is recapitulated by chemical perturbation using an allosteric CBFß-inhibitor. Strikingly, we show that RUNX1 addiction extends to other genetic subtypes of pediatric B-ALL and also adult disease. Importantly, inhibition of RUNX1 activity spares normal hematopoiesis. Our results suggest that chemical intervention in the RUNX1 program may provide a therapeutic opportunity in ALL.

Longitudinal single-cell transcriptomics reveals distinct patterns of recurrence in acute myeloid leukemia.

BACKGROUND: Acute myeloid leukemia (AML) is a heterogeneous and aggressive blood cancer that results from diverse genetic aberrations in the hematopoietic stem or progenitor cells (HSPCs) leading to the expansion of blasts in the hematopoietic system. The heterogeneity and evolution of cancer blasts can render therapeutic interventions ineffective in a yet poorly understood patient-specific manner. In this study, we investigated the clonal heterogeneity of diagnosis (Dx) and relapse (Re) pairs at genetic and transcriptional levels, and unveiled the underlying pathways and genes contributing to recurrence. METHODS: Whole-exome sequencing was used to detect somatic mutations and large copy number variations (CNVs). Single cell RNA-seq was performed to investigate the clonal heterogeneity between Dx-Re pairs and amongst patients. RESULTS: scRNA-seq analysis revealed extensive expression differences between patients and Dx-Re pairs, even for those with the same -presumed- initiating events. Transcriptional differences between and within patients are associated with clonal composition and evolution, with the most striking differences in patients that gained large-scale copy number variations at relapse. These differences appear to have significant molecular implications, exemplified by a DNMT3A/FLT3-ITD patient where the leukemia switched from an AP-1 regulated clone at Dx to a mTOR signaling driven clone at Re. The two distinct AML1-ETO pairs share genes related to hematopoietic stem cell maintenance and cell migration suggesting that the Re leukemic stem cell-like (LSC-like) cells evolved from the Dx cells. CONCLUSIONS: In summary, the single cell RNA data underpinned the tumor heterogeneity not only amongst patient blasts with similar initiating mutations but also between each Dx-Re pair. Our results suggest alternatively and currently unappreciated and unexplored mechanisms leading to therapeutic resistance and AML recurrence.

Dietary methionine starvation impairs acute myeloid leukemia progression.

Targeting altered tumor cell metabolism might provide an attractive opportunity for patients with acute myeloid leukemia (AML). An amino acid dropout screen on primary leukemic stem cells and progenitor populations revealed a number of amino acid dependencies, of which methionine was one of the strongest. By using various metabolite rescue experiments, nuclear magnetic resonance-based metabolite quantifications and 13C-tracing, polysomal profiling, and chromatin immunoprecipitation sequencing, we identified that methionine is used predominantly for protein translation and to provide methyl groups to histones via S-adenosylmethionine for epigenetic marking. H3K36me3 was consistently the most heavily impacted mark following loss of methionine. Methionine depletion also reduced total RNA levels, enhanced apoptosis, and induced a cell cycle block. Reactive oxygen species levels were not increased following methionine depletion, and replacement of methionine with glutathione or N-acetylcysteine could not rescue phenotypes, excluding a role for methionine in controlling redox balance control in AML. Although considered to be an essential amino acid, methionine can be recycled from homocysteine. We uncovered that this is primarily performed by the enzyme methionine synthase and only when methionine availability becomes limiting. In vivo, dietary methionine starvation was not only tolerated by mice, but also significantly delayed both cell line and patient-derived AML progression. Finally, we show that inhibition of the H3K36-specific methyltransferase SETD2 phenocopies much of the cytotoxic effects of methionine depletion, providing a more targeted therapeutic approach. In conclusion, we show that methionine depletion is a vulnerability in AML that can be exploited therapeutically, and we provide mechanistic insight into how cells metabolize and recycle methionine.

CBX2 shapes chromatin accessibility promoting AML via p38 MAPK signaling pathway.

BACKGROUND: The dynamic epigenome and proteins specialized in the interpretation of epigenetic marks critically contribute to leukemic pathogenesis but also offer alternative therapeutic avenues. Targeting newly discovered chromatin readers involved in leukemogenesis may thus provide new anticancer strategies. Accumulating evidence suggests that the PRC1 complex member CBX2 is overexpressed in solid tumors and promotes cancer cell survival. However, its role in leukemia is still unclear. METHODS: We exploited reverse genetic approaches to investigate the role of CBX2 in human leukemic cell lines and ex vivo samples. We also analyzed phenotypic effects following CBX2 silencing using cellular and molecular assays and related functional mechanisms by ATAC-seq and RNA-seq. We then performed bioinformatic analysis of ChIP-seq data to explore the influence of histone modifications in CBX2-mediated open chromatin sites. Lastly, we used molecular assays to determine the contribution of CBX2-regulated pathways to leukemic phenotype. RESULTS: We found CBX2 overexpressed in leukemia both in vitro and ex vivo samples compared to CD34+ cells. Decreased CBX2 RNA levels prompted a robust reduction in cell proliferation and induction of apoptosis. Similarly, sensitivity to CBX2 silencing was observed in primary acute myeloid leukemia samples. CBX2 suppression increased genome-wide chromatin accessibility followed by alteration of leukemic cell transcriptional programs, resulting in enrichment of cell death pathways and downregulation of survival genes. Intriguingly, CBX2 silencing induced epigenetic reprogramming at p38 MAPK-associated regulatory sites with consequent deregulation of gene expression. CONCLUSIONS: Our results identify CBX2 as a crucial player in leukemia progression and highlight a potential druggable CBX2-p38 MAPK network in AML.

Direct targeted therapy for MLL-fusion-driven high-risk acute leukaemias.

BACKGROUND: Improving the poor prognosis of infant leukaemias remains an unmet clinical need. This disease is a prototypical fusion oncoprotein-driven paediatric cancer, with MLL (KMT2A)-fusions present in most cases. Direct targeting of these driving oncoproteins represents a unique therapeutic opportunity. This rationale led us to initiate a drug screening with the aim of discovering drugs that can block MLL-fusion oncoproteins. METHODS: A screen for inhibition of MLL-fusion proteins was developed that overcomes the traditional limitations of targeting transcription factors. This luciferase reporter-based screen, together with a secondary western blot screen, was used to prioritize compounds. We characterized the lead compound, disulfiram (DSF), based on its efficient ablation of MLL-fusion proteins. The consequences of drug-induced MLL-fusion inhibition were confirmed by cell proliferation, colony formation, apoptosis assays, RT-qPCR, in vivo assays, RNA-seq and ChIP-qPCR and ChIP-seq analysis. All statistical tests were two-sided. RESULTS: Drug-induced inhibition of MLL-fusion proteins by DSF resulted in a specific block of colony formation in MLL-rearranged cells in vitro, induced differentiation and impeded leukaemia progression in vivo. Mechanistically, DSF abrogates MLL-fusion protein binding to DNA, resulting in epigenetic changes and down-regulation of leukaemic programmes setup by the MLL-fusion protein. CONCLUSION: DSF can directly inhibit MLL-fusion proteins and demonstrate antitumour activity both in vitro and in vivo, providing, to our knowledge, the first evidence for a therapy that directly targets the initiating oncogenic MLL-fusion protein.

Presence of mutant p53 increases stem cell frequency and is associated with reduced binding to classic TP53 binding sites in cell lines and primary AMLs.

With an overall 5%-10% incidence rate in acute myeloid leukemia (AML), the occurrence of TP53 mutations is low compared with that in solid tumors. However, when focusing on high-risk groups including secondary AML (sAML) and therapy-related AMLs, the frequency of mutations reaches up to 35%. Mutations may include loss of heterozygosity (LOH) or deletion of the 17p allele, but are mostly missense substitutions that are located in the DNA-binding domain. Despite elaborate research on the effects of TP53 mutations in solid tumors, in hematological malignancies, the effects of TP53 mutations versus loss of TP53 remain unclear and under debate. Here, we compared the cellular effects of a TP53 mutant and loss of TP53 in human hematopoietic stem and progenitor cells (HSPCs). We found that when expressing TP53 mutant or loss of TP53 using siRNA, CD34+/CD38- cells have a significantly enhanced replating potential, which could not be demonstrated for the CD34+/CD38+ population. Using RNA-sequencing analysis, we found a loss of expression of p53 target genes in cells with TP53 knockdown. In contrast, an increased expression of a large number of genes was observed when expressing TP53 mutant, resulting in an increase in expression of genes involved in megakaryocytic differentiation, plasma membrane binding, and extracellular structure organization. When binding of p53 wild type and p53 mutant was compared in cell lines, we found that mutant p53 binds to a large number of binding sites genomewide, contrary to wild-type p53, for which binding is restricted to genes with a p53 binding motif. These findings were verified in primary AMLs with and without mutated TP53. In conclusion, in our models, we identified overlapping effects of TP53 mutant and loss of TP53 on in vitro stem cell properties but distinct effects on DNA binding and gene expression.

Dynamics of broad H3K4me3 domains uncover an epigenetic switch between cell identity and cancer-related genes.

Broad domains of H3K4 methylation have been associated with consistent expression of tissue-specific, cell identity, and tumor suppressor genes. Here, we identified broad domain-associated genes in healthy human thymic T cell populations and a collection of T cell acute lymphoblastic leukemia (T-ALL) primary samples and cell lines. We found that broad domains are highly dynamic throughout T cell differentiation, and their varying breadth allows the distinction between normal and neoplastic cells. Although broad domains preferentially associate with cell identity and tumor suppressor genes in normal thymocytes, they flag key oncogenes in T-ALL samples. Moreover, the expression of broad domain-associated genes, both coding and noncoding, is frequently deregulated in T-ALL. Using two distinct leukemic models, we showed that the ectopic expression of T-ALL oncogenic transcription factor preferentially impacts the expression of broad domain-associated genes in preleukemic cells. Finally, an H3K4me3 demethylase inhibitor differentially targets T-ALL cell lines depending on the extent and number of broad domains. Our results show that the regulation of broad H3K4me3 domains is associated with leukemogenesis, and suggest that the presence of these structures might be used for epigenetic prioritization of cancer-relevant genes, including long noncoding RNAs.

Characterization of a genomic region 8 kb downstream of GFI1B associated with myeloproliferative neoplasms.

A genomic locus 8 kb downstream of the transcription factor GFI1B (Growth Factor Independence 1B) predisposes to clonal hematopoiesis and myeloproliferative neoplasms. One of the most significantly associated polymorphisms in this region is rs621940-G. GFI1B auto-represses GFI1B, and altered GFI1B expression contributes to myeloid neoplasms. We studied whether rs621940-G affects GFI1B expression and growth of immature cells. GFI1B ChIP-seq showed clear binding to the rs621940 locus. Preferential binding of various hematopoietic transcription factors to either the rs621940-C or -G allele was observed, but GFI1B showed no preference. In gene reporter assays the rs621940 region inhibited GFI1B promoter activity with the G-allele having less suppressive effects compared to the C-allele. However, CRISPR-Cas9 mediated deletion of the locus in K562 cells did not alter GFI1B expression nor auto-repression. In healthy peripheral blood mononuclear cells GFI1B expression did not differ consistently between the rs621940 alleles. Long range and targeted deep sequencing did not detect consistent effects of rs621940-G on allelic GFI1B expression either. Finally, we observed that myeloid colony formation was not significantly affected by either rs621940 allele in 193 healthy donors. Together, these findings show no evidence that rs621940 or its locus affect GFI1B expression, auto-repression or growth of immature myeloid cells.