Lee, Jeong‐Heon

The intestinal epithelium is exposed to diverse combinations of microbiota-derived compounds; however, the mechanisms by which the host integrates these signals remain poorly defined. Studying two highly abundant microbial metabolites, we identified the purine metabolite hypoxanthine as an effector metabolite that directly drives signaling and the short-chain fatty acid butyrate as a regulatory metabolite that conditions host responsiveness. Specifically, hypoxanthine activates the adenosine A1 receptor–TRPC4 axis in enterochromaffin (EC) cells, triggering calcium influx and serotonin release resulting in accelerated gastrointestinal transit locally and increased platelet activation systemically. In contrast, butyrate epigenetically upregulates specific G protein–coupled receptors and ion channels to enhance response to hypoxanthine and the neurotransmitters norepinephrine and dopamine. These findings define a cooperative signaling framework and highlight the role of EC cells as a distinct epithelial signaling hub that senses and integrates microbial metabolite signals to drive physiological responses. Our findings provide a mechanistic foundation for therapeutic strategies that leverage combinatorial microbial signaling.

Abstract Repetitive injury is hypothesized to lead to progressive tissue fibrosis and end-stage organ failure. Whether tissue-resident mesenchymal cell populations retain epigenetic memory of prior injuries that contribute to this pathological process is unknown. Here we used a genetic lineage labeling approach to mark the lung mesenchyme prior to injury, then performed multimodal analyses on isolated lung mesenchyme during the initiation, progression, and resolution of the fibrotic response. Our results demonstrate the remarkable epigenetic and transcriptional plasticity of the lung mesenchyme during fibrotic activation and de-activation. Despite this plasticity, we also find that the lung mesenchyme exhibits an enhanced fibrotic program upon reinjury. We identify RUNX1 as a critical driver of both fibrotic activation and fibrotic memory. Comparison of fresh isolated and cultured lung mesenchyme demonstrates that RUNX1 is spontaneously activated in standard culture conditions, previously masking these roles of RUNX1. Targeted knockdown of RUNX1 dampens fibrotic mesenchymal cell activation immediately after cell isolation, but with reduced efficacy after only days of culture, confirming its functional importance to both early activation and long-term memory. Collectively, our findings implicate RUNX1 in the initiation and memory of fibrotic mesenchymal cell activation that together prime enhanced mesenchymal cell responses upon repeated injury.

Obesity is an epidemic with myriad health effects, but little is understood regarding individual obese phenotypes and how they may respond to therapy. Epigenetic changes associated with obesity have been detected in blood, liver, pancreas, and adipose tissues. Previous work using human organoids found that dietary glucose hyperabsorption is a steadfast trait in cultures derived from some obese subjects, but detailed transcriptional or epigenomic features of the intestinal epithelia associated with this persistent phenotype are unknown. This study evaluated differentially expressed genes and relative chromatin accessibility in intestinal organoids established from donors classified as non-obese, obese, or obese hyperabsorptive by body mass index and glucose transport assays. Transcriptomic analysis indicated that obese hyperabsorptive subject organoids have significantly upregulated dietary nutrient absorption transcripts and downregulated type I interferon targets. Chromatin accessibility and transcription factor footprinting predicted that enhanced HNF4G binding may promote the obese hyperabsorption phenotype. Quantitative RT-PCR assessment in organoids representing a larger subject cohort suggested that intestinal epithelial expression of CUBN, GIP, SLC5A11, and SLC2A5 were highly correlated with hyperabsorption. Thus, the obese hyperabsorption phenotype was characterized by transcriptional changes that support increased nutrient uptake by intestinal epithelia, potentially driven by differentially accessible chromatin. Recognizing unique intestinal phenotypes in obesity provides a new perspective in considering therapeutic targets and options with which to manage the disease.

Abstract Background: DNMT3A and TET2 are epigenetic regulator genes commonly mutated in age-related clonal hematopoiesis (CH). Despite having opposed epigenetic functions, these mutations are associated with increased all-cause mortality and a low risk for hematological neoplasms. In this study, we make use of a natural inflammatory response occurring during COVID-19 to define the impact of these mutations on inflammation and outcomes. Methods: A cohort of 243 community-based patients (Olmsted County, MN) with COVID-19 were identified. CH status was determined on PBMC using error corrected sequencing (lower detection limit 0.5% VAF). In a subset of patients, we use multi-omics interrogation of the genome, transcriptome, and epigenome at single cell (sc) resolution along with DNA methylation, in PBMC, from patients with COVID-19 and CH. Results: We detected 97 CH mutations in 72 (29.6%) of 243 patients with COVID-19. The most frequent being DNMT3A (n=30, 30%) and TET2 (n=26, 28%). The presence of CH negatively impacted overall survival when restricted to TET2 and DNMT3A mutations. However, only DNMT3Amt CH retained age and comorbidity independent prognostication (HR 2.84, p = 0.022). We then assessed for changes in DNA methylation using Illumina MethylationEPIC array on patients with COVID-19 and DNMT3A (n=3) or TET2mt (n=4) CH. Site specific DNA methylation changes revealed significant loss of methylation in DNMT3Amt CH in comparison to TET2mt patients, with 10,944 hypomethylated and 1,160 hypermethylated sites. We found that actively transcribed states (Tx, TxWk) were more commonly hypomethylated in DNMT3Amt CH while enhancers (Enh) and promoters (TssA, TssAFlnk) were more commonly hypomethylated in TET2mt CH. We then conducted scDNA-seq in DNMT3Amt (n=3) and TET2mt (n=1) CH in COVID-19 and observed that while in TET2mt CH, mutations were largely restricted to classical and intermediate monocytes, in DNMT3Amt CH, mutations were seen in myeloid and lymphoid lineage cells. We performed scRNA-seq on patients with COVID-19 and DNMT3A (n=3) or TET2mt (n=6) CH. Differential gene expression analysis identified 1,569 upregulated genes in DNMT3Amt patients and 205 downregulated. Pathway analysis demonstrated an upregulation of genes involved in T cell survival, function, and inflammation in DNMT3Amt patients compared to TET2mt patients with a significant overexpression of GIMAP1, GIMAP4, IRF2, IL32 in multiple cell types. We correlated transcriptomic changes with chromatin accessibility using the 10X Genomics Multiome platform to profile both gene expression and open chromatin from the same samples. Analysis of global distribution of cut sites and differentially accessible peaks showed increased chromatin accessibility in DNMT3Amt CH, especially CD4+ T lymphocytes and NK cells. Conclusion: We demonstrate an age, sex, and comorbidity-independent adverse impact on mortality associated with DNMT3Amt CH in the context of COVID-19 and highlight transcriptomic and epigenetic differences between DNMT3A and TET2mt CH, with single cell resolution. Citation Format: Jenna A Fernandez, Wazim Mohammed Ismail, Moritz Binder, Terra Lasho, Susan M. Geyer, Amelia Mazzone, Christy M. Finke, Abhishek Mangaonkar, Jeong-Heon Lee, Liguo Wang, Vernadette A. Simon, Fariborz Rakhshan Rohakthar, Amik Munankarmy, Jonathan J. Harrington, Melissa R. Snyder, Keith D. Robertson, Alexandre Gaspar-Maia, Mrinal Patnaik. Divergent impacts of DNMT3A and TET2 mutant clonal hematopoiesis on COVID-19 related outcomes [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: DNA Methylation, Clonal Hematopoiesis, and Cancer; 2025 Feb 1-4; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2025;85(3 Suppl):Abstract nr A012.

Introduction: ASXL1-mutant Chronic Myelomonocytic Leukemia (CMML) is a high-risk myeloid neoplasm characterized by epigenetic remodeling and overexpression of leukemogenic driver genes. DLK1 (delta-like homolog 1) is a non-canonical Notch ligand, crucial for the proliferation of early hematopoietic progenitor cells and its expression is controlled by maternal imprinting. DLK1 expression is restricted to select progenitor populations, but was found to be overexpressed in several malignancies including ASXL1-mutant CMML. Here we investigate the molecular mechanisms that drive the aberrant re-expression of DLK1 and postulate potential targeted therapeutic interventions. Methods: We performed bulk RNA- and ChIP-seq (H3K4me1/3, H3K27ac/me3, H2AK119ub, H2BK20ac) on bone marrow mononuclear cells from 20 patients with CMML (11 ASXL1-mutant, 9 -wildtype). All libraries were sequenced on an Illumina NovaSeq instrument. Differential gene expression and histone modification occupancy between ASXL1-mutant and -wildtype patients were assessed. Single-cell RNA-seq was performed for a subset of patients to investigate the cell populations of interest for DLK1 expression. To assess DLK1 allelic expression, we identified a single nucleotide polymorphism (SNP) in the coding region of DLK1 (rs1802710) and used cases which were confirmed heterozygous by Sanger sequencing. DLK1 transcript levels were quantified by RT-PCR, and the relative abundance of alleles was measured using digital droplet PCR. We designed primers to amplify the region around our SNP of interest and utilized high mismatch discrimination locked nucleic acid (LNA) probes labeled with FAM (reference SNP allele) and HEX (alternate SNP allele) to quantify specific alleles. The number of RNA molecules originating from each allele was quantified by counting the number of droplets that fluoresced with the corresponding probe colors. Illumina 850K EPIC microarrays were used to evaluate DNA methylation. Drug studies were carried out in progenitor colony forming assays derived from ASXL1-mutant patients (bone marrow mononuclear cells) with the dual BRD4 and p300 / CBP inhibitor EP31670 (62 nM, 24-hour exposure), quantifying DLK1 expression by real-time RT-PCR. Results: This study included 20 patients with CMML (median age at diagnosis 74 years). All ASXL1 mutations were typical hotspot mutations (most commonly Gly646Trpfs*12) and the variant allele frequencies were compatible with heterozygosity (median 0.41). The spectrum of co-mutations was typical and included TET2 (23%), SRSF2 (16%), CBL (7%), among others. DLK1 was strongly overexpressed in ASXL1-mutant patients (log2 fold-change +4.0, FDR < 0.001). Allele-specific gene expression studies (n = 3) revealed the re-expression of the imprinted maternal DLK1 allele in ASXL1-mutant patients. Single-cell gene expression studies (n = 5) demonstrated that the DLK1 overexpression is restricted to the hematopoietic progenitor cell population. Differential methylation analysis between patients with ASXL1-mutant clonal hematopoiesis (n = 4) and matched controls (n = 8) identified two hypermethylated CpG islands in the DLK1 promoter region (cg18122767, Δß = +8.6%, p = 0.021 and cg15831296, Δß = +5.8%, p = 0.008). Histone modification analysis identified DLK1 as a p300 target gene by mapping H2BK20ac onto established ASXL1-mutant-specific enhancer regions (H3K4me1 + H3K27ac co-occupancy). The distal enhancer EH38D2760182 was predicted to bind p300 and to govern the expression of DLK1. ChIP-seq experiment in K-562 and Kasumi-1 cells showed the binding of BRD4 and p300 to EH38D2760182, respectively. Targeting BRD4 and p300 using EP31670 led to the downregulation of DLK1 expression (mean fold-change = 0.569, p < 0.001, n = 3) in progenitor colony forming assays. Conclusions: We observed the aberrant re-expression of the maternally imprinted DLK1 allele in patients with ASXL1-mutant CMML. Both DLK1 promoter hypermethylation and the genotype-specific activation of enhancers by p300 were identified as putative regulatory mechanisms. The distal enhancer EH38D2760182 was found to bind BRD4 and p300 and the pharmacologic disruption of BRD4 / p300 activity in model systems led to the downregulation of DLK1. The clinical-grade dual inhibitor of BRD and CBP / p300 EP31670 is currently being investigated in patients with ASXL1-mutant chronic myeloid neoplasms (NCT05488548).

Abstract DNMT3A and TET2 are epigenetic regulator genes commonly mutated in age-related clonal hematopoiesis (CH). Despite having opposed epigenetic functions, these mutations are associated with increased all-cause mortality and a low risk for progression to hematologic neoplasms. Although individual impacts on the epigenome have been described using different model systems, the phenotypic complexity in humans remains to be elucidated. Here, we make use of a natural inflammatory response occurring during coronavirus disease 2019 (COVID-19), to understand the association of these mutations with inflammatory morbidity (acute respiratory distress syndrome [ARDS]) and mortality. We demonstrate the age-independent, negative impact of DNMT3A mutant (DNMT3Amt) CH on COVID-19–related ARDS and mortality. Using single-cell proteogenomics we show that DNMT3A mutations involve myeloid and lymphoid lineage cells. Using single-cell multiomics sequencing, we identify cell-specific gene expression changes associated with DNMT3A mutations, along with significant epigenomic deregulation affecting enhancer accessibility, resulting in overexpression of interleukin-32 (IL-32), a proinflammatory cytokine that can result in inflammasome activation in monocytes and macrophages. Finally, we show with single-cell resolution that the loss of function of DNMT3A is directly associated with increased chromatin accessibility in mutant cells. Hence, we demonstrate the negative prognostic impact of DNMT3Amt CH on COVID-19–related ARDS and mortality. DNMT3Amt CH in the context of COVID-19, was associated with inflammatory transcriptional priming, resulting in overexpression of IL32. This overexpression was secondary to increased chromatic accessibility, specific to DNMT3Amt CH cells. DNMT3Amt CH can thus serve as a potential biomarker for adverse outcomes in COVID-19.

Abstract Background Malignant peripheral nerve sheath tumors (MPNSTs) can arise from atypical neurofibromas (ANF). Loss of the polycomb repressor complex 2 (PRC2) is a common event. Previous studies on PRC2-regulated genes in MPNST used genetic add-back experiments in highly aneuploid MPNST cell lines which may miss PRC2-regulated genes in NF1-mutant ANF-like precursor cells. A set of PRC2-regulated genes in human Schwann cells (SCs) has not been defined. We hypothesized that PRC2 loss has direct and indirect effects on gene expression resulting in MPNST, so we sought to identify PRC2-regulated genes in immortalized human Schwann cells (iHSCs). Methods We engineered NF1-deficient iHSCs with loss of function SUZ12 or EED mutations. RNA sequencing revealed 1327 differentially expressed genes to define PRC2-regulated genes. To investigate MPNST pathogenesis, we compared genes in iHSCs to consistent gene expression differences between ANF and MPNSTs. Chromatin immunoprecipitation sequencing was used to further define targets. Methylome and proteomic analyses were performed to further identify enriched pathways. Results We identified potential PRC2-regulated drivers of MPNST progression. Pathway analysis indicates many upregulated cancer-related pathways. We found transcriptional evidence for activated Notch and Sonic Hedgehog (SHH) signaling in PRC2-deficient iHSCs. Functional studies confirm that Notch signaling is active in MPNST cell lines, patient-derived xenografts, and transient cell models of PRC2 deficiency. A combination of MEK and γ-secretase inhibition shows synergy in MPNST cell lines. Conclusions We identified PRC2-regulated genes and potential drivers of MPNSTs. Our findings support the Notch pathway as a druggable target in MPNSTs. Our identification of PRC2-regulated genes and pathways could result in more novel therapeutic approaches.

Abstract Background DNMT3A and TET2 are epigenetic regulator genes commonly mutated in age-related clonal hematopoiesis (CH). Despite having opposing epigenetic functions, these mutations are associated with increased all-cause mortality and a low risk for progression to hematological neoplasms. While individual impacts on the epigenome have been described using different model systems, the phenotypic complexity in humans remains to be elucidated. Results Here we make use of a natural inflammatory response occurring during coronavirus disease 2019 (COVID-19), to understand the association of these mutations with inflammatory morbidity and mortality. We demonstrate the age-independent, negative impact of DNMT3A mutant CH on COVID-19-related cytokine release severity and mortality. Using single cell proteogenomics we show that DNMT3A mutations involve myeloid and lymphoid cells. Using single cell multiomics sequencing, we identify cell-specific gene expression changes associated with DNMT3A mutations, along with significant epigenomic deregulation affecting enhancer accessibility, resulting in overexpression of IL32, a proinflammatory cytokine that can result in inflammasome activation in monocytes and macrophages. Finally, we show with single cell resolution that the loss of function of DNMT3A is directly associated with increased chromatin accessibility in mutant cells. Conclusions We demonstrate the negative prognostic impact of DNMT3Amt CH on COVID-19 related inflammatory morbidity and mortality. DNMT3Amt CH involves myeloid and lymphoid cells and in the context of COVID-19, was associated with inflammatory transcriptional priming, resulting in overexpression of IL32. This overexpression was secondary to increased chromatic accessibility, specific to DNMT3Amt CH cells. DNMT3Amt CH can serve as a potential biomarker for adverse inflammatory outcomes.

Abstract Neurofibromatosis type 1 (NF1) syndrome is an autosomal dominant cancer predisposition syndrome defined by germline deletion of one NF1 allele. Somatic loss of heterozygosity in the remaining NF1 allele in Schwann lineage cells gives rise to benign tumors known as plexiform neurofibromas. These tumors have a high risk of transforming into malignant peripheral nerve sheath tumors (MPNST) through the sequential deletion of CDKN2A/B followed by recurrent mutations in SUZ12 or EED, constituting the inactivation of polycomb repressor complex 2 (PRC2). PRC2 is a transcriptional repressor complex involved in the silencing of several development and differentiation genes. Loss of this complex leads to a depletion in trimethylation marks on Histone H3 Lysine 27 (H3K27me3) and the subsequent gain of acetylation marks on the same lysine residue (H3K27Ac). These molecular hallmarks are shared with diffuse intrinsic pontine gliomas (DIPG) harboring the H3K27M mutation which prevents effective deposition of H3K27me3, thus having similar molecular phenotypes to MPNST. However, the full range of epigenetic consequences of these events are currently not well characterized. Additionally, given that MPNSTs originate from cells of neural crest origin and DIPGs from cells of neural origin, it is not know to what degree epigenetic changes stemming from H3K27me3-loss are general across tumors or are context-dependent based on cell-of-origin. In order to contribute to this effort, the Largaespada lab has engineered Schwann cell lines to harbor deletions in NF1 and either SUZ12 or EED and assayed the H3K27ac landscape of these cells in addition to patient-derived MPNST cell lines. Utilizing the ranked ordering of super enhancers (ROSE) algorithm, we have identified common super enhancers within each MPNST cell line. These super enhancer elements are located nearby a variety of putative oncogenes responsible for the maintenance of cell identity. These elements were then compared to publicly available data from DIPG primary and patient-derived samples to understand common regulatory events driven by loss of H3K27me3 as well as molecular signatures present only in MPNST or in DIPG that may be a result of PRC2 or H3K27M alone. We further characterized these super enhancers by assessing the expression levels of their associated genes to understand their role in transcriptional activation. These efforts provide a roadmap of the super enhancer landscape of both MPNSTs and DIPGs to investigate novel mechanisms of H3K27me3-dependent tumorigenesis while also contributing to knowledge of potential effects of H3K27me3-loss that dependent on the cell-of-origin for each respective tumor. Citation Format: Christopher M. Stehn, Minu Bhunia, Kyle Williams, Alex T. Larsson, Ariel Caride, Jeong-Heon Lee, Tamas Ordog, David Largaespada. Characterizing the super enhancer landscape of H3K27me3 loss in MPNST and DIPG [abstract]. In: Proceedings of the AACR Special Conference on Brain Cancer; 2023 Oct 19-22; Minneapolis, Minnesota. Philadelphia (PA): AACR; Cancer Res 2024;84(5 Suppl_1):Abstract nr B009.

Introduction: Clonal dendritic cell (DC) aggregates in chronic myelomonocytic leukemia (CMML) are seen in approximately 20-30% patients at diagnosis. Previously, these were thought to be comprised of CD123+ plasmacytoid DCs but our data has suggested that these are heterogenous and comprise of myeloid DCs, monocytes, and myeloid-derived suppressor cells. Furthermore, these aggregates uniformly express immune checkpoints such as indoleamine 2,3-dioxygenase-1 (IDO), associate with an IDO-specific systemic impact (increased trytophan catabolism), and expanded regulatory T cell population in CMML, suggesting a role in inducing immune tolerance. DC aggregates are also associated with increased risk of acute myeloid leukemia (AML) transformation and resistance to hypomethylating agents, highlighting that they are important for disease progression in CMML (Mangaonkar A et al. Blood Adv 2020). We performed this study to decipher the architecture of a CMML DC aggregate to understand the complex inter-cellular interactions, and assess their impact on the adaptive immune system. Methods: In this study, we used archival fresh-frozen paraffin-embedded (FFPE) bone marrow biopsy specimens obtained from CMML patients, and isolated monocytes (CD14+ magnetic bead sorting) and lymphocytes (CD14-ve fraction) from untreated CMML and age-matched normal control (NC) peripheral blood mononuclear cells (PBMCs). DCs were cultured from monocytes using a commercial protocol. Spatial architecture of the CMML BM microenvironment was interrogated through multiplex imaging of FFPE samples by co-detection by indexing (CODEX) using 30 markers (CCR6, CD107a, CD11c, CD20, CD21, CD31, CD3e, CD4, CD44, CD45, CD45RO CD56, CD68, CD8, E-Cadherin, FoxP3, HLA-DR, CD14, CD16, CD11b, CD33, CD15, CD66b, CD123, CD303, CD34, CD163 and IDO). Single-cell transcriptomic signatures were analyzed using the 10X single cell Fixed RNA profiling platform, and validated using bulk RNA seq (Illumina stranded mRNA prep) from isolated and cultured cells. Impact of DCs on lymphocytes was assessed via autologous mixed DC-T culture experiments. Results: BM FFPE samples from 5 CMML cases (3 with DC aggregates) were included for multiplex imaging via CODEX. Results confirm the significant interactions of CD8+ and CD4+ T cells with CD123+/HLA-DR+ DCs within a CMML DC aggregate (Fig. 1A). Upon quantitative assessment of cell-cell interactions, IDO+/CD8+T ( P=0.006), IDO+/MDSC ( P=0.005), and IDO+/macrophage ( P=0.001) interactions were significantly increased in CMML samples with versus without DC aggregates . Specific findings were validated using single-stain immunohistochemistry and multiplex immunoflourescence staining. Next, in order to assess transcriptomic signatures within BM at the single-cell level, we chose two cases; CMML (BM aspirate blast% 1, WBC 49.6 x 10 9/L, ASXL1/GATA2/SETBP1/SRSF2 mutations) and AML transformed from CMML (BM aspirate blast% 81, WBC 152 x 10 9/L, NPM1/DNMT3A/TET2 mutations). Approximately, 10,000 cells were obtained and analyzed with cell type annotation through human PBMC reference (Azimuth). Results show an expanded regulatory T cell population in AML compared to CMML (UMAP plot, Fig. 1B), suggestive of immune tolerance. In order to validate findings and perform functional experiments, we used sorted monocytes and cutured DCs from CMML and NC. DC phenotype was confirmed both through visual inspection and flow cytometry, and clonality was confirmed through Sanger sequencing (for previously identified mutations). Differential gene expression analysis through bulk RNAseq showed significant differences between CMML and NC lymphocytes with immunoinhibitory receptor, LILRB2 (log2fold change=7.05, adjusted P=0.0002) being one of the top five significantly upregulated genes, further confirming immune tolerance. Finally, mixed autologous DC-lymphocyte reaction (1:4 ratio) experiments demonstrated that DCs induce T cell tolerance in CMML (expansion of exhausted PD1/TIM3/LAG3+ CD8 T cells) when compared to NC. Conclusion: Clonal dendritic cell aggregates within the CMML bone marrow microenviroment interact with T cells and induce immune tolerance. This represents a novel mechanism of disease progression in CMML with potential therapeutic targets (IDO, LILRB2, and/or CD123).