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.

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 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.

Histone‐lysine N‐methyltransferase SETD2 (SETD2), the sole histone methyltransferase that catalyzes trimethylation of lysine 36 on histone H3 (H3K36me3), is often mutated in clear cell renal cell carcinoma (ccRCC). SETD2 mutation and/or loss of H3K36me3 is linked to metastasis and poor outcome in ccRCC patients. Epithelial‐to‐mesenchymal transition (EMT) is a major pathway that drives invasion and metastasis in various cancer types. Here, using novel kidney epithelial cell lines isogenic for SETD2, we discovered that SETD2 inactivation drives EMT and promotes migration, invasion, and stemness in a transforming growth factor‐beta‐independent manner. This newly identified EMT program is triggered in part through secreted factors, including cytokines and growth factors, and through transcriptional reprogramming. RNA‐seq and assay for transposase‐accessible chromatin sequencing uncovered key transcription factors upregulated upon SETD2 loss, including SOX2, POU2F2 (OCT2), and PRRX1, that could individually drive EMT and stemness phenotypes in SETD2 wild‐type (WT) cells. Public expression data from SETD2 WT/mutant ccRCC support the EMT transcriptional signatures derived from cell line models. In summary, our studies reveal that SETD2 is a key regulator of EMT phenotypes through cell‐intrinsic and cell‐extrinsic mechanisms that help explain the association between SETD2 loss and ccRCC metastasis.

The snATAC + snRNA platform allows epigenomic profiling of open chromatin and gene expression with single-cell resolution. The most critical assay step is to isolate high-quality nuclei to proceed with droplet-base single nuclei isolation and barcoding. With the increasing popularity of multiomic profiling in various fields, there is a need for optimized and reliable nuclei isolation methods, mainly for human tissue samples. Herein we compared different nuclei isolation methods for cell suspensions, such as peripheral blood mononuclear cells (PBMC, n = 18) and a solid tumor type, ovarian cancer (OC, n = 18), derived from debulking surgery. Nuclei morphology and sequencing output parameters were used to evaluate the quality of preparation. Our results show that NP-40 detergent-based nuclei isolation yields better sequencing results than collagenase tissue dissociation for OC, significantly impacting cell type identification and analysis. Given the utility of applying such techniques to frozen samples, we also tested frozen preparation and digestion (n = 6). A paired comparison between frozen and fresh samples validated the quality of both specimens. Finally, we demonstrate the reproducibility of scRNA and snATAC + snRNA platform, by comparing the gene expression profiling of PBMC. Our results highlight how the choice of nuclei isolation methods is critical for obtaining quality data in multiomic assays. It also shows that the measurement of expression between scRNA and snRNA is comparable and effective for cell type identification.

Abstract Understanding the function of non-coding genomic sequence variants represents a challenge for biomedicine. Many diseases are products of gene-by-environment interactions with complex mechanisms. This study addresses these themes by mechanistic characterization of non-coding variants that influence gene expression only after drug or hormone exposure. Using glucocorticoid signaling as a model system, we integrated genomic, transcriptomic, and epigenomic approaches to unravel mechanisms by which variant function could be revealed by hormones or drugs. Specifically, we identified cis-regulatory elements and 3D interactions underlying ligand-dependent associations between variants and gene expression. One-quarter of the glucocorticoid-modulated variants that we identified had already been associated with clinical phenotypes. However, their affected genes were ‘unmasked’ only after glucocorticoid exposure and often with function relevant to the disease phenotypes. These diseases involved glucocorticoids as risk factors or therapeutic agents and included autoimmunity, metabolic and mood disorders, osteoporosis and cancer. For example, we identified a novel breast cancer risk gene, MAST4, with expression that was repressed by glucocorticoids in cells carrying the risk genotype, repression that correlated with MAST4 expression in breast cancer and treatment outcomes. These observations provide a mechanistic framework for understanding non-coding genetic variant-chemical environment interactions and their role in disease risk and drug response.

AbstractAlternative polyadenylation (APA) is a common genetic regulatory mechanism that generates distinct 3′ ends for RNA transcripts. Changes in APA have been associated with multiple biological processes and disease phenotypes. However, the role of hormones and their drug analogs in APA remains largely unknown. In this study, we investigated transcriptome‐wide the impact of glucocorticoids on APA in 30 human B‐lymphoblastoid cell lines. We found that glucocorticoids could regulate APA for a subset of genes, possibly by changing the expression of 142 RNA‐binding proteins, some with known APA‐regulating properties. Interestingly, genes with glucocorticoid‐mediated APA were enriched in viral translation‐related pathways, while genes with glucocorticoid‐mediated expression were enriched in interferon and interleukin pathways, suggesting that glucocorticoid‐mediated APA might result in functional consequences distinct from gene expression. For example, glucocorticoids, a pharmacotherapy for severe COVID‐19, were found to change the APA but not the expression of LY6E, an important antiviral inhibitor in coronavirus diseases. Glucocorticoid‐mediated APA was also cell‐type‐specific, suggesting an action of glucocorticoids that may be unique to immune regulation. We also observed evidence for genotype‐dependent glucocorticoid‐mediated APA (referred to as pharmacogenomic‐alterative polyadenylation quantitative trait loci), providing potential functional mechanisms for a series of common genetic variants that had previously been associated with immune disorders, but without a clear mechanism. In summary, this study reports a series of observations regarding the impact of glucocorticoids on APA, raising the possibility that this mechanism might have implications for both disease pathophysiology and drug therapy.

Abstract The human genome encodes large numbers of non-coding RNAs, including divergent antisense transcripts at transcription start sites (TSSs). However, molecular mechanisms by which divergent antisense transcription is regulated have not been detailed. Here, we report a novel ZWC complex composed of ZC3H4, WDR82 and CK2 that suppresses divergent antisense transcription. The ZWC complex preferentially localizes at TSSs of active genes through direct interactions of ZC3H4 and WDR82 subunits with the S5p RNAPII C-terminal domain. ZC3H4 depletion leads to increased divergent antisense transcription, especially at genes that naturally produce divergent antisense transcripts. We further demonstrate that the ZWC complex phosphorylates the previously uncharacterized N-terminal acidic domain of SPT5, a subunit of the transcription-elongation factor DSIF, and that this phosphorylation is responsible for suppressing divergent antisense transcription. Our study provides evidence that the newly identified ZWC-DSIF axis regulates the direction of transcription during the transition from early to productive elongation.