The cell type specific 5hmC landscape and dynamics of healthy human hematopoiesis and TET2-mutant pre-leukemia

Deep5hmC: predicting genome-wide 5-hydroxymethylcytosine landscape via a multimodal deep learning model

MOTIVATION

5-Hydroxymethylcytosine (5hmC), a crucial epigenetic mark with a significant role in regulating tissue-specific gene expression, is essential for understanding the dynamic functions of the human genome. Despite its importance, predicting 5hmC modification across the genome remains a challenging task, especially when considering the complex interplay between DNA sequences and various epigenetic factors such as histone modifications and chromatin accessibility.

RESULTS

Using tissue-specific 5hmC sequencing data, we introduce Deep5hmC, a multimodal deep learning framework that integrates both the DNA sequence and epigenetic features such as histone modification and chromatin accessibility to predict genome-wide 5hmC modification. The multimodal design of Deep5hmC demonstrates remarkable improvement in predicting both qualitative and quantitative 5hmC modification compared to unimodal versions of Deep5hmC and state-of-the-art machine learning methods. This improvement is demonstrated through benchmarking on a comprehensive set of 5hmC sequencing data collected at four developmental stages during forebrain organoid development and across 17 human tissues. Compared to DeepSEA and random forest, Deep5hmC achieves close to 4% and 17% improvement of Area Under the Receiver Operating Characteristic (AUROC) across four forebrain developmental stages, and 6% and 27% across 17 human tissues for predicting binary 5hmC modification sites; and 8% and 22% improvement of Spearman correlation coefficient across four forebrain developmental stages, and 17% and 30% across 17 human tissues for predicting continuous 5hmC modification. Notably, Deep5hmC showcases its practical utility by accurately predicting gene expression and identifying differentially hydroxymethylated regions (DhMRs) in a case-control study of Alzheimer's disease (AD). Deep5hmC significantly improves our understanding of tissue-specific gene regulation and facilitates the development of new biomarkers for complex diseases.

AVAILABILITY AND IMPLEMENTATION

Deep5hmC is available via https://github.com/lichen-lab/Deep5hmC.

Selective advantage of mutant stem cells in human clonal hematopoiesis is associated with attenuated response to inflammation and aging

Clonal hematopoiesis (CH) arises when hematopoietic stem cells (HSCs) acquire mutations, most frequently in the DNMT3A and TET2 genes, conferring a competitive advantage through mechanisms that remain unclear. To gain insight into how CH mutations enable gradual clonal expansion, we used single-cell multi-omics with high-fidelity genotyping on human CH bone marrow (BM) samples. Most of the selective advantage of mutant cells occurs within HSCs. DNMT3A- and TET2-mutant clones expand further in early progenitors, while TET2 mutations accelerate myeloid maturation in a dose-dependent manner. Unexpectedly, both mutant and non-mutant HSCs from CH samples are enriched for inflammatory and aging transcriptomic signatures, compared with HSCs from non-CH samples, revealing a non-cell-autonomous effect. However, DNMT3A- and TET2-mutant HSCs have an attenuated inflammatory response relative to wild-type HSCs within the same sample. Our data support a model whereby CH clones are gradually selected because they are resistant to the deleterious impact of inflammation and aging.

Ascorbate depletion increases quiescence and self-renewal potential in hematopoietic stem cells and multipotent progenitors

Ascorbate (vitamin C) limits hematopoietic stem cell (HSC) function and suppresses leukemia development by promoting the function of the Tet2 tumor suppressor. In humans, ascorbate is obtained from the diet while in mice it is synthesized in the liver. In this study, we show that deletion of the Slc23a2 ascorbate transporter severely depleted ascorbate from hematopoietic cells. Slc23a2 deficiency increased HSC reconstituting potential and self-renewal potential upon transplantation into irradiated mice. Slc23a2 deficiency also increased the reconstituting and self-renewal potential of multipotent hematopoietic progenitors (MPPs), conferring the ability to long-term reconstitute irradiated mice. Slc23a2-deficient HSCs and MPPs divided much less frequently than control HSCs and MPPs. Increased self-renewal and reconstituting potential were observed particularly in quiescent Slc23a2-deficient HSCs and MPPs. The effect of Slc23a2 deficiency on MPP self-renewal was not mediated by reduced Tet2 function. Ascorbate thus regulates quiescence and restricts self-renewal potential in HSCs and MPPs such that ascorbate depletion confers MPPs with long-term self-renewal potential.

Machine learning identifies cell-free DNA 5-hydroxymethylation biomarkers that detect occult colorectal cancer in PLCO Screening Trial subjects

Background

Colorectal cancer (CRC) is a leading cause of cancer-related mortality, and CRC detection through screening improves survival rates. A promising avenue to improve patient screening compliance is the development of minimally-invasive liquid biopsy assays that target CRC biomarkers on circulating cell-free DNA (cfDNA) in peripheral plasma. In this report, we identify cfDNA biomarker candidate genes bearing the epigenetic mark 5-hydroxymethylcytosine (5hmC) that diagnose occult CRC up to 36 months prior to clinical diagnosis using the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial samples.

Methods

Archived PLCO Trial plasma samples containing cfDNA were obtained from the National Cancer Institute (NCI) biorepositories. Study subjects included those who were diagnosed with CRC within 36 months of blood collection (i.e., case, n = 201) and those who were not diagnosed with any cancer during an average of 16.3 years of follow-up (i.e., controls, n = 402). Following the extraction of 3 - 8 ng cfDNA from less than 300 microliters plasma, we employed the sensitive 5hmC-Seal chemical labeling approach, followed by next-generation sequencing (NGS). We then conducted association studies and machine-learning modeling to analyze the genome-wide 5hmC profiles within training and validation groups that were randomly selected at a 2:1 ratio.

Results

Despite the technical challenges associated with the PLCO samples (e.g., limited plasma volumes, low cfDNA amounts, and long archival times), robust genome-wide 5hmC profiles were successfully obtained from these samples. Association analyses using the Cox proportional hazards models suggested several epigenetic pathways relevant to CRC development distinguishing cases from controls. A weighted Cox model, comprised of 32-associated gene bodies, showed predictive detection value for CRC as early as 24-36 months prior to overt tumor presentation, and a trend for increased predictive power was observed for blood samples collected closer to CRC diagnosis. Notably, the 5hmC-based predictive model showed comparable performance regardless of sex and self-reported race/ethnicity, and significantly outperformed risk factors such as age and obesity according to BMI (body mass index). Additionally, further improvement of predictive performance was achieved by combining the 5hmC-based model and risk factors for CRC.

Conclusions

An assay of 5hmC epigenetic signals on cfDNA revealed candidate biomarkers with the potential to predict CRC occurrence despite the absence of clinical symptoms or the availability of effective predictors. Developing a minimally-invasive clinical assay that detects 5hmC-modified biomarkers holds promise for improving early CRC detection and ultimately patient survival through higher compliance screening and earlier intervention. Future investigation to expand this strategy to prospectively collected samples is warranted.

Transient inhibition of 53BP1 increases the frequency of targeted integration in human hematopoietic stem and progenitor cells

Genome editing by homology directed repair (HDR) is leveraged to precisely modify the genome of therapeutically relevant hematopoietic stem and progenitor cells (HSPCs). Here, we present a new approach to increasing the frequency of HDR in human HSPCs by the delivery of an inhibitor of 53BP1 (named "i53") as a recombinant peptide. We show that the use of i53 peptide effectively increases the frequency of HDR-mediated genome editing at a variety of therapeutically relevant loci in HSPCs as well as other primary human cell types. We show that incorporating the use of i53 recombinant protein allows high frequencies of HDR while lowering the amounts of AAV6 needed by 8-fold. HDR edited HSPCs were capable of long-term and bi-lineage hematopoietic reconstitution in NSG mice, suggesting that i53 recombinant protein might be safely integrated into the standard CRISPR/AAV6-mediated genome editing protocol to gain greater numbers of edited cells for transplantation of clinically meaningful cell populations.

IDH1-mutant preleukemic hematopoietic stem cells can be eliminated by inhibition of oxidative phosphorylation

Rare preleukemic hematopoietic stem cells (pHSCs) harboring only the initiating mutations can be detected at the time of AML diagnosis. pHSCs are the origin of leukemia and a potential reservoir for relapse. Using primary human samples and gene-editing to model isocitrate dehydrogenase 1 (IDH1) mutant pHSCs, we show epigenetic, transcriptional, and metabolic differences between pHSCs and healthy hematopoietic stem cells (HSCs). We confirm that IDH1 driven clonal hematopoiesis is associated with cytopenia, suggesting an inherent defect to fully reconstitute hematopoiesis. Despite giving rise to multilineage engraftment, IDH1-mutant pHSCs exhibited reduced proliferation, blocked differentiation, downregulation of MHC Class II genes, and reprogramming of oxidative phosphorylation metabolism. Critically, inhibition of oxidative phosphorylation resulted in complete eradication of IDH1-mutant pHSCs but not IDH2-mutant pHSCs or wildtype HSCs. Our results indicate that IDH1-mutant preleukemic clones can be targeted with complex I inhibitors, offering a potential strategy to prevent development and relapse of leukemia.

Preleukemic single-cell landscapes reveal mutation-specific mechanisms and gene programs predictive of AML patient outcomes

Acute myeloid leukemia (AML) and myeloid neoplasms develop through acquisition of somatic mutations that confer mutation-specific fitness advantages to hematopoietic stem and progenitor cells. However, our understanding of mutational effects remains limited to the resolution attainable within immunophenotypically and clinically accessible bulk cell populations. To decipher heterogeneous cellular fitness to preleukemic mutational perturbations, we performed single-cell RNA sequencing of eight different mouse models with driver mutations of myeloid malignancies, generating 269,048 single-cell profiles. Our analysis infers mutation-driven perturbations in cell abundance, cellular lineage fate, cellular metabolism, and gene expression at the continuous resolution, pinpointing cell populations with transcriptional alterations associated with differentiation bias. We further develop an 11-gene scoring system (Stem11) on the basis of preleukemic transcriptional signatures that predicts AML patient outcomes. Our results demonstrate that a single-cell-resolution deep characterization of preleukemic biology has the potential to enhance our understanding of AML heterogeneity and inform more effective risk stratification strategies.

TET2-mediated DNA hydroxymethylation of TGFB1 is related to selective intrauterine growth restriction in monochorionic twin pregnancies

INTRODUCTION

Selective intrauterine growth restriction (sIUGR), which specifically occurs in monochorionic (MC) twins, usually has a poor prognosis and the underlying mechanisms are not well understood. It is an ideal model for exploring epigenetic-modified mechanisms for fetal development in MCDA twins due to eliminating the interference of different heritable backgrounds and intrauterine environments among individuals.

METHODS

The levels of ten-eleven translocation 2 (TET2) and its upstream and downstream targets miR-29b-3p and transforming growth factor beta 1 (TGFB1) were determined using RT‒qPCR, western blotting, and immunohistochemistry. Using TET2 overexpression and knockdown methods, we investigated the role of TET2 in trophoblast functions. The regulatory relationships among TET2, miR-29b-3p, and TGFB1 were explored by cell migration assay, invasion assay, apoptotic ratio assays, Western blot, hMeDIP-qPCR and dual-luciferase assay.

RESULTS

A consistent upregulation of TET2 and TGFB1 was observed in the smaller placental shares compared to the larger placental shares in sIUGR. Gain-of-function studies of TET2 in trophoblasts showed decreased cell invasion and increased apoptosis, whereas loss-of-function studies of TET2 rescued this effect. Mechanistic studies revealed that miR-29b-3p and TGFB1 were the upstream factor and downstream target of TET2, respectively. Furthermore, miR-29b-3p/TET2/TGFB1-smad was identified as a unique axis that regulates trophoblast invasion, migration, and apoptosis in a DNA hydroxymethylation-dependent manner.

DISCUSSION

We elucidated the functional roles of TET2 and DNA hydroxymethylation in trophoblasts and identified a novel DNA regulatory mechanism, providing a basis for further exploration of DNA epigenetic regulatory patterns in sIUGR.

RUNX1 loss renders hematopoietic and leukemic cells dependent on IL-3 and sensitive to JAK inhibition

Disease-initiating mutations in the transcription factor RUNX1 occur as germline and somatic events that cause leukemias with particularly poor prognosis. However, the role of RUNX1 in leukemogenesis is not fully understood, and effective therapies for RUNX1-mutant leukemias remain elusive. Here, we used primary patient samples and a RUNX1-KO model in primary human hematopoietic cells to investigate how RUNX1 loss contributes to leukemic progression and to identify targetable vulnerabilities. Surprisingly, we found that RUNX1 loss decreased proliferative capacity and stem cell function. However, RUNX1-deficient cells selectively upregulated the IL-3 receptor. Exposure to IL-3, but not other JAK/STAT cytokines, rescued RUNX1-KO proliferative and competitive defects. Further, we demonstrated that RUNX1 loss repressed JAK/STAT signaling and rendered RUNX1-deficient cells sensitive to JAK inhibitors. Our study identifies a dependency of RUNX1-mutant leukemias on IL-3/JAK/STAT signaling, which may enable targeting of these aggressive blood cancers with existing agents.

Cellular Composition and 5hmC Signature Predict the Treatment Response of AML Patients to Azacitidine Combined with Chemotherapy

Azacitidine (AZA) is a DNA methyltransferase inhibitor and epigenetic modulator that can be an effective agent in combination with chemotherapy for patients with high-risk acute myeloid leukemia (AML). However, biological factors driving the therapeutic response of such hypomethylating agent (HMA)-based therapies remain unknown. Herein, the transcriptome and/or genome-wide 5-hydroxymethylcytosine (5hmC) is characterized for 41 patients with high-risk AML from a phase 1 clinical trial treated with AZA epigenetic priming followed by high-dose cytarabine and mitoxantrone (AZA-HiDAC-Mito). Digital cytometry reveals that responders have elevated Granulocyte-macrophage-progenitor-like (GMP-like) malignant cells displaying an active cell cycle program. Moreover, the enrichment of natural killer (NK) cells predicts a favorable outcome in patients receiving AZA-HiDAC-Mito therapy or other AZA-based therapies. Comparing 5hmC profiles before and after five-day treatment of AZA shows that AZA exposure induces dose-dependent 5hmC changes, in which the magnitude correlates with overall survival (p = 0.015). An extreme gradient boosting (XGBoost) machine learning model is developed to predict the treatment response based on 5hmC levels of 11 genes, achieving an area under the curve (AUC) of 0.860. These results suggest that cellular composition markedly impacts the treatment response, and showcase the prospect of 5hmC signatures in predicting the outcomes of HMA-based therapies in AML.

Loss of TET2 in human hematopoietic stem cells alters the development and function of neutrophils

Somatic mutations commonly occur in hematopoietic stem cells (HSCs). Some mutant clones outgrow through clonal hematopoiesis (CH) and produce mutated immune progenies shaping host immunity. Individuals with CH are asymptomatic but have an increased risk of developing leukemia, cardiovascular and pulmonary inflammatory diseases, and severe infections. Using genetic engineering of human HSCs (hHSCs) and transplantation in immunodeficient mice, we describe how a commonly mutated gene in CH, TET2, affects human neutrophil development and function. TET2 loss in hHSCs produce a distinct neutrophil heterogeneity in bone marrow and peripheral tissues by increasing the repopulating capacity of neutrophil progenitors and giving rise to low-granule neutrophils. Human neutrophils that inherited TET2 mutations mount exacerbated inflammatory responses and have more condensed chromatin, which correlates with compact neutrophil extracellular trap (NET) production. We expose here physiological abnormalities that may inform future strategies to detect TET2-CH and prevent NET-mediated pathologies associated with CH.

Hacking hematopoiesis - emerging tools for examining variant effects

Hematopoiesis is a continuous process of blood and immune cell production. It is orchestrated by thousands of gene products that respond to extracellular signals by guiding cell fate decisions to meet the needs of the organism. Although much of our knowledge of this process comes from work in model systems, we have learned a great deal from studies on human genetic variation. Considerable insight has emerged from studies on presumed monogenic blood disorders, which continue to provide key insights into the mechanisms critical for hematopoiesis. Furthermore, the emergence of large-scale biobanks and cohorts has uncovered thousands of genomic loci associated with blood cell traits and diseases. Some of these blood cell trait-associated loci act as modifiers of what were once thought to be monogenic blood diseases. However, most of these loci await functional validation. Here, we discuss the validation bottleneck and emerging methods to more effectively connect variant to function. In particular, we highlight recent innovations in genome editing, which have paved the path forward for high-throughput functional assessment of loci. Finally, we discuss existing barriers to progress, including challenges in manipulating the genomes of primary hematopoietic cells.

Myelodysplastic Syndrome associated TET2 mutations affect NK cell function and genome methylation

Myelodysplastic syndromes (MDS) are clonal hematopoietic disorders, representing high risk of progression to acute myeloid leukaemia, and frequently associated to somatic mutations, notably in the epigenetic regulator TET2. Natural Killer (NK) cells play a role in the anti-leukemic immune response via their cytolytic activity. Here we show that patients with MDS clones harbouring mutations in the TET2 gene are characterised by phenotypic defects in their circulating NK cells. Remarkably, NK cells and MDS clones from the same patient share the TET2 genotype, and the NK cells are characterised by increased methylation of genomic DNA and reduced expression of Killer Immunoglobulin-like receptors (KIR), perforin, and TNF-α. In vitro inhibition of TET2 in NK cells of healthy donors reduces their cytotoxicity, supporting its critical role in NK cell function. Conversely, NK cells from patients treated with azacytidine (#NCT02985190; https://clinicaltrials.gov/ ) show increased KIR and cytolytic protein expression, and IFN-γ production. Altogether, our findings show that, in addition to their oncogenic consequences in the myeloid cell subsets, TET2 mutations contribute to repressing NK-cell function in MDS patients.

Clinical and molecular characteristics associated with Vitamin C deficiency in myeloid malignancies; real world data from a prospective cohort

Vitamin C is an essential vitamin that acts as a co-factor for many enzymes involved in epigenetic regulation in humans. Low vitamin C levels in hematopoietic stem cells (HSC) promote self-renewal and vitamin C supplementation retards leukaemogenesis in vitamin C-deficient mouse models. Studies on vitamin C levels in patients with myeloid malignancies are limited. We thus conducted a retrospective analysis on a prospective cohort of patients with myeloid malignancies on whom plasma vitamin C levels were measured serially at diagnosis and during treatment. Baseline characteristics including hematological indices, cytogenetics, and molecular mutations are described in this cohort. Among 64 patients included in our study, 11 patients (17%) had low vitamin C levels. We noted a younger age at diagnosis for patients with myeloid malignancies who had low plasma vitamin C levels. Patients with low plasma vitamin C levels were more likely to have acute myeloid leukemia compared to other myeloid malignancies. Low vitamin C levels were associated with ASXL1 mutations. Our study calls for further multi-institutional studies to understand the relevance of low plasma vitamin C level in myeloid neoplasms, the role of vitamin C deficiency in leukemogenesis, and the potential benefit of vitamin C supplementation.

Viewing AML through a New Lens: Technological Advances in the Study of Epigenetic Regulation

Epigenetic modifications, such as histone modifications and DNA methylation, are essential for ensuring the dynamic control of gene regulation in every cell type. These modifications are associated with gene activation or repression, depending on the genomic context and specific type of modification. In both cases, they are deposited and removed by epigenetic modifier proteins. In acute myeloid leukemia (AML), the function of these proteins is perturbed through genetic mutations (i.e., in the DNA methylation machinery) or translocations (i.e., MLL-rearrangements) arising during leukemogenesis. This can lead to an imbalance in the epigenomic landscape, which drives aberrant gene expression patterns. New technological advances, such as CRISPR editing, are now being used to precisely model genetic mutations and chromosomal translocations. In addition, high-precision epigenomic editing using dCas9 or CRISPR base editing are being used to investigate the function of epigenetic mechanisms in gene regulation. To interrogate these mechanisms at higher resolution, advances in single-cell techniques have begun to highlight the heterogeneity of epigenomic landscapes and how these impact on gene expression within different AML populations in individual cells. Combined, these technologies provide a new lens through which to study the role of epigenetic modifications in normal hematopoiesis and how the underlying mechanisms can be hijacked in the context of malignancies such as AML.

Reduced PU.1 Expression Collaborates with Tet2 Loss to Trigger Myeloid Leukemogenesis

The leukemic transformation of hematopoietic stem and progenitor cells in the setting of Tet2 deficiency is driven by PU.1 gene network loss through complementary reduction in PU.1 expression and hypermethylation of ETS loci at the enhancers of PU.1 target genes. See related article by Aivalioti et al., p. 444 (6).

Clonal Hematopoiesis: Role in Hematologic and Non-Hematologic Malignancies

Hematopoietic stem cells (HSCs) ensure the coordinated and balanced production of all hematopoietic cell types throughout life. Aging is associated with a gradual decline of the self-renewal and regenerative potential of HSCs and with the development of clonal hematopoiesis. Clonal hematopoiesis of indeterminate potential (CHIP) defines the clonal expansion of genetically variant hematopoietic cells bearing one or more gene mutations and/or structural variants (such as copy number alterations). CHIP increases exponentially with age and is associated with cancers, including hematologic neoplasia, cardiovascular and other diseases. The presence of CHIP consistently increases the risk of hematologic malignancy, particularly in individuals who have CHIP in association with peripheral blood cytopenia.