Open Chromatin Research Articles

Overcoming the Silencing of Doxycycline-Inducible Promoters in hiPSC-derived Cardiomyocytes

Background Human induced pluripotent stem cells (hiPSCs) are pivotal for studying human development, modeling diseases, and advancing regenerative medicine. Effective control of transgene expression is crucial to achieve temporal and quantitative precision in all of these contexts. The doxycycline (dox)-inducible OPTi-OX system, which integrates the Tet-On 3G transactivator and dox-responsive transgene at the hROSA26 and AAVS1 genomic safe harbors (GSHs), respectively, offers a promising solution. Yet, transgene silencing, particularly in hiPSC-derived cardiomyocytes (hiPSC-CMs), limits its utility. Methods To address this, we evaluated strategies to enhance dox-inducible transgene expression. We compared two promoters, TRE3VG and T11, for activity and stability, and investigated the addition of a Ubiquitous Chromatin Opening Element (UCOE) to reduce silencing. We also tested relocating the transgene cassette to the CLYBL GSH, and employed sodium butyrate (SB), a histone deacetylase inhibitor, to restore promoter activity. Transgene expression was assessed via flow cytometry and real-time quantitative PCR. Results TRE3VG exhibited higher activity than T11, but both were prone to silencing. UCOE did not enhance promoter activity in hiPSCs, but modestly reduced silencing in hiPSC-CMs. Targeting the CLYBL locus improved promoter activity compared to AAVS1 in both hiPSCs and hiPSC-CMs. SB restored activity in silenced inducible promoters within hiPSC-CMs, but compromised hiPSC viability. Unexpectedly, Tet-On 3G was silenced in some clones and could not be reactivated by SB. Conclusions These findings underscore the need for integrating multiple strategies, including careful GSH selection, improved cassette design, epigenetic modulation, and clone screening, to develop robust dox-inducible systems that retain functionality during hiPSC differentiation.

Linking Iris Cis-Regulatory Variants to Primary Angle-Closure Glaucoma Via Clinical Imaging and Multiomics.

To elucidate the genetic basis of primary angle-closure glaucoma (PACG) by identifying pathogenic tissue and critical tissue-specific variants. The correlations among PACG susceptibility, axial length (AL), and anterior chamber depth (ACD) were evaluated using meta-analyses. Propensity score matching was utilized on 2161 participants from the Handan Eye Study to determine the risk factors independent of ACD and AL for PACG. Subsequently, we employed the assay for transposase-accessible chromatin with sequencing (ATAC-seq) and allele-specific self-transcribing active regulatory region sequencing (STARR-seq) to screen 202 PACG genome-wide association study (GWAS) variants for chromatin accessibility and functional roles. The meta-analysis found that PACG susceptibility loci are not associated with ACD or AL. However, abnormal iris phenotypes emerged as significant independent risk factors for primary angle-closure disease (PACD), unrelated to ACD and AL. Substantial enrichment of PACG heritability was observed in the open chromatin regions of the human iris. Within the iris-relevant cellular context, 22 out of the 202 PACG GWAS variants could influence enhancer activity. Two variants in the iris open chromatin regions were implicated in the modulation of PLEKHA7 and C10orf53 expression. The downregulation of these two genes affects cytoskeletal organization. Our findings underscore the importance of the iris in the pathogenesis of PACG and identified iris-specific, enhancer-modulating variants that may influence disease risk. Our approach also provides a generalizable framework for studying ocular diseases from the perspective of enhancer-modulating variants.

Silent silencers: Leave no trace

The regulatory mechanisms of silencers have remained poorly understood. In this issue, Hofbauer etal.1 conduct a genome-wide screen in Drosophila melanogaster and reveal three silencer types that appear to work alone-without the need for combinatorial action, traditional chromatin marks, or open chromatin regions.

OCRClassifier: integrating statistical control chart into machine learning framework for better detecting open chromatin regions.

Open chromatin regions (OCRs) play a crucial role in transcriptional regulation and gene expression. In recent years, there has been a growing interest in using plasma cell-free DNA (cfDNA) sequencing data to detect OCRs. By analyzing the characteristics of cfDNA fragments and their sequencing coverage, researchers can differentiate OCRs from non-OCRs. However, the presence of noise and variability in cfDNA-seq data poses challenges for the training data used in the noise-tolerance learning-based OCR estimation approach, as it contains numerous noisy labels that may impact the accuracy of the results. For current methods of detecting OCRs, they rely on statistical features derived from typical open and closed chromatin regions to determine whether a region is OCR or non-OCR. However, there are some atypical regions that exhibit statistical features that fall between the two categories, making it difficult to classify them definitively as either open or closed chromatin regions (CCRs). These regions should be considered as partially open chromatin regions (pOCRs). In this paper, we present OCRClassifier, a novel framework that combines control charts and machine learning to address the impact of high-proportion noisy labels in the training set and classify the chromatin open states into three classes accurately. Our method comprises two control charts. We first design a robust Hotelling T2 control chart and create new run rules to accurately identify reliable OCRs and CCRs within the initial training set. Then, we exclusively utilize the pure training set consisting of OCRs and CCRs to create and train a sensitized T2 control chart. This sensitized T2 control chart is specifically designed to accurately differentiate between the three categories of chromatin states: open, partially open, and closed. Experimental results demonstrate that under this framework, the model exhibits not only excellent performance in terms of three-class classification, but also higher accuracy and sensitivity in binary classification compared to the state-of-the-art models currently available.

Transposable Elements Contribute to the Regulation of Long Noncoding RNAs in Drosophila melanogaster

Background: Transposable elements (TEs) and noncoding sequences are major components of the genome, yet their functional contributions to long noncoding RNAs (lncRNAs) are not well understood. Although many lncRNAs originating from TEs (TE-lncRNAs) have been identified across various organisms, their characteristics and regulatory roles, particularly in insects, remain largely unexplored. This study integrated multi-omics data to investigate TE-lncRNAs in D. melanogaster, focusing on the influence of transposons across different omics levels. Results: We identified 16,118 transposons overlapping with lncRNA sequences that constitute 2119 TE-lncRNAs (40.4% of all lncRNAs) using 256 public RNA-seq samples and 15 lncRNA-seq samples of Drosophila S2 cells treated with heavy metals. Of these, 67.2% of TE-lncRNAs contain more than one TE. The LTR/Gypsy family was the most common transposon insertion. Transposons preferred to insert into promoters, transcription starting sites, and intronic regions, especially in chromosome ends. Compared with lncRNAs, TE-lncRNAs showed longer lengths, a lower conservation, and lower levels but a higher specificity of expression. Multi-omics data analysis revealed positive correlations between transposon insertions and chromatin openness at the pre-transcriptional level. Notably, a total of 516 TE-lncRNAs provided transcriptional factor binding sites through transposon insertions. The regulatory network of a key transcription factor was rewired by transposons, potentially recruiting other transcription factors to exert regulatory functions under heavy metal stress. Additionally, 99 TE-lncRNAs were associated with m6A methylation modification sites, and 115 TE-lncRNAs potentially provided candidate small open reading frames through transposon insertions. Conclusions: Our data analysis demonstrated that TEs contribute to the regulation of lncRNAs. TEs not only promote the transcriptional regulation of lncRNAs, but also facilitate their post-transcriptional and epigenetic regulation.

A multi-ancestry genome-wide association study identifies novel candidate loci in the RARB gene associated with hypertensive disorders of pregnancy

Genetic factors related to pregnancy-related traits are understudied, especially in ancestrally diverse cohorts. To assess maternal contributions to hypertensive disorders of pregnancy (HDP, we performed a multi-ancestry genome-wide association study (GWAS) of HDP in data from the North Carolina-based Personalized Environment and Genes Study (PEGS) cohort with validation in the UK Biobank (UKBB). The GWAS revealed two maternal loci associated with HDP at the genome-wide significance level. The lead independent variants were rs114954125 on chromosome 2 (near LRP1B; OR (95% CI): 2.96 (2.02,4.34); P=2.82 x 10-8) and rs61176331 on chromosome 3 (on RARB; OR (95% CI): 3.08 (2.12,4.48); P=3.52 x 10-9). We validated the associations near RARB with a meta-analysis of PEGS and the UK Biobank. We also identified cis-eQTLs in the candidate region associated with decreased RARB expression in macrophage cells exposed to Salmonella. Chromatin mapping in FUMA identified a significant interaction within chromosome 3's enhancer and open chromatin regions, with strong effects observed for RARB and H3P10 gene regulation in mesendoderm cells, mesenchymal stem cells, and trophoblast-like stem cells. We applied existing polygenic scores (PGS) for preeclampsia and gestational hypertension and found the scores were significantly associated with HDP in PEGS. The findings demonstrate the power of multi-ancestry studies for genetic discovery and highlight the relationship between immune response, regulation, and HDP and the utility of PGS for risk prediction.

Phospho-regulation of ASCL1-mediated chromatin opening during cellular reprogramming.

The proneural transcription factor ASCL1 regulates neurogenesis and drives somatic cell reprogramming into neurons. However, not all cell types can be reprogrammed by ASCL1, raising the questions of what provides competence and how we can overcome barriers to enable directed differentiation. Here, we investigate how levels of ASCL1 and its phosphorylation modulate its activity over progressive lineage restriction of embryonic stem cells. We find that inhibition of ASCL1 phosphorylation enhances reprogramming of both mesodermal and neuroectodermal cells, while pluripotent cells remain refractory to ASCL1-directed neuronal differentiation. By performing RNA-seq and ATAC-seq in neuroectoderm, we find that un(der)phosphorylated ASCL1 causes increased chromatin accessibility at sites proximal to neuronal genes, accompanied by their increased expression. Combined analysis of protein stability and proneural function of phosphomutant and phosphomimetic ASCL1 reveals that protein stability plays only a marginal role in regulating activity, while changes in amino acid charge cannot fully explain enhanced activity of the serine-proline mutant variants of ASCL1. Our work provides new insights into proneural factor activity and regulation and suggests ways to optimize reprogramming protocols in cancer and regenerative medicine.

The co-receptor Neuropilin-1 enhances proliferation in inv(16) acute myeloid leukemia via VEGF signaling.

Oncogenic programs regulate the proliferation and maintenance of cancer stem cells, and can define pharmacologic dependencies. In acute myeloid leukemia (AML) with the chromosome inversion 16 (inv(16)), the fusion oncoproteinCBFβ::MYH11 regulates pathways associated with leukemia stem cell activity. Here we demonstrate that expression of Neuropilin-1 (NRP1) is regulated by the fusion oncoprotein, and promotes AML expansion. Mechanistically, we show that the NRP1 locus has open chromatin in inv(16) AML, and that CBFβ::MYH11 modulates the local function of the transcription factors ERG, GATA2 and RUNX1 to sustain NRP1 levels. We found that ERG activates NRP1 expression, and that CBFβ::MYH11 knockdown represses ERG expression, thereby allowing the repressive activity of GATA2/RUNX1 at three NRP1 enhancers. Functionally, we demonstrate that NRP1 enhances the expansion of leukemic cells in vitro and in mice, and that this activity is dependent on its VEGFR-associated FV/FVIII domain. Finally, we show that treatment with VEGF inhibitor axitinib reduces AML cell growth and delays median leukemia latency in vivo. Our findings reveal that the NRP1/VEGF axis mediates proliferation in inv(16) AML blasts, and suggest that targeting NRP1 function could be promising in combination AML therapy.

Aberrant expression of histone H2B variants reshape chromatin and alter oncogenic gene expression programs.

Chromatin architecture governs DNA accessibility and gene expression. Thus, any perturbations to chromatin can significantly alter gene expression programs and promote disease. Prior studies demonstrate that every amino acid in a histone is functionally significant, and that even a single amino acid substitution can drive specific cancers. We previously observed that naturally occurring H2B variants are dysregulated during the epithelial to mesenchymal transition (EMT) in bronchial epithelial cells. Naturally occurring H2B variants differ from canonical H2B by only a few amino acids, yet single amino acid changes in other histone variants (e.g., H3.3) can drive cancer. We therefore hypothesized that H2B variants might function like oncohistones, and investigated how they modify chromatin architecture, dynamics, and function. We find that H2B variants are frequently dysregulated in many cancers, and correlate with patient prognosis. Despite high sequence similarity, mutations in each H2B variant tend to occur at specific "hotspots" in cancer. Some H2B variants cause tighter DNA wrapping around nucleosomes, leading to more compact chromatin structures and reduced transcription factor accessibility to nucleosomal DNA. They also altered genome-wide accessibility to oncogenic regulatory elements and genes, with concomitant changes in oncogenic gene expression programs. Although we did not observe changes in cell proliferation or migration in vitro , our Gene Ontology (GO) analyses of ATAC-seq peaks and RNA-seq data indicated significant changes in oncogenic pathways. These findings suggest that H2B variants may influence early-stage, cancer-associated regulatory mechanisms, potentially setting the stage for oncogenesis later on. Thus, H2B variant expression could serve as an early cancer biomarker, and H2B variants might be novel therapeutic targets.

Elevated EGR1 Binding at Enhancers in Excitatory Neurons Correlates with Neuronal Subtype-Specific Epigenetic Regulation.

Brain development and neuronal cell specification are accompanied with epigenetic changes to achieve diverse gene expression regulation. Interacting with cell-type specific epigenetic marks, transcription factors bind to different sets of cis-regulatory elements in different types of cells. Currently, it remains largely unclear how cell-type specific gene regulation is achieved for neurons. In this study, we generated epigenetic maps to perform comparative histone modification analysis between excitatory and inhibitory neurons. We found that neuronal cell-type specific histone modifications are enriched in super enhancer regions containing abundant EGR1 motifs. Further CUT&RUN data validated that more EGR1 binding sites can be detected in excitatory neurons and primarily located in enhancers. Integrative analysis revealed that EGR1 binding is strongly correlated with various epigenetic markers for open chromatin regions and associated with distinct gene pathways with neuronal subtype-specific functions. In inhibitory neurons, the majority of genomic regions hosting EGR1 binding sites become accessible at early embryonic stages. In contrast, the super enhancers in excitatory neurons hosting EGR1 binding sites gained their accessibility during postnatal stages. This study highlights the significance of transcription factor binding to enhancer regions, which may play a crucial role in establishing cell-type specific gene regulation in neurons.

Rewired chromatin structure and epigenetic gene dysregulation during HTLV-1 infection to leukemogenesis.

Human T-cell leukemia virus type 1 (HTLV-1) broadly impacts host genes, affecting the infected cell population and inducing the development of a disease with a poor prognosis, adult T-cell leukemia-lymphoma (ATL). This study aimed to provide a comprehensive epigenomic characterization of the infected cell population and evaluated the transcriptome and chromatin structures of peripheral blood cells in HTLV-1-infected individuals using RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin sequencing (ATAC-seq). The infected cells showed significant changes in gene expression patterns from the polyclonal stage and before ATL onset while demonstrating similarities to tumor-forming ATL cells. These similarities were a result of large-scale open chromatin changes, supporting the independent early formation of epigenomic aberrations as an underlying mechanism for later clonal propagation. This study also demonstrated that HTLV-1 Tax directly affects the host chromatin structure, thereby developing fundamental epigenomic characteristics. Several Tax target genes, including the RASGRP3-ERK pathway, were recognized, indicating an impact on signaling pathways. This genome-wide variability in chromatin structural property is a novel feature of HTLV-1 infection and may contribute to pathogenic mechanisms. In addition, it has crucial implications for better understanding the impact of HTLV-1 on the host genome and identifying novel therapeutic targets.

ZBTB48 is a priming factor regulating B-cell-specific CIITA expression.

The class-II transactivator (CIITA) is the master regulator of MHC class-II gene expression and hence the adaptive immune response. Three cell type-specific promoters (pI, pIII, and pIV) are involved in the regulation of CIITA expression, which can be induced by IFN-γ in non-immune cells. While key regulatory elements have been identified within these promoters, our understanding of the transcription factors regulating CIITA expression is incomplete. Here, we demonstrate that the telomere-binding protein and transcriptional activator ZBTB48 directly binds to both critical activating elements within the B-cell-specific promoter CIITA pIII. ZBTB48 knockout impedes the CIITA/MHC-II expression program induced in non-APC cells by IFN-γ, and loss of ZBTB48 in mice silences MHC-II expression in pro-B and immature B cells. Transcriptional regulation of CIITA by ZBTB48 is enabled by ZBTB48-dependent chromatin opening at CIITA pIII upstream of activating H3K4me3 marks. We conclude that ZBTB48 primes CIITA pIII by acting as a molecular on-off-switch for B-cell-specific CIITA expression.

Integrated single-cell analysis reveals distinct epigenetic-regulated cancer cell states and a heterogeneity-guided core signature in tamoxifen-resistant breast cancer.

Inter- and intra-tumor heterogeneity is considered a significant factor contributing to the development of endocrine resistance in breast cancer. Recent advances in single-cell RNA sequencing (scRNA-seq) and single-cell ATAC sequencing (scATAC-seq) allow us to explore inter- and intra-tumor heterogeneity at single-cell resolution. However, such integrated single-cell analysis has not yet been demonstrated to characterize the transcriptome and chromatin accessibility in breast cancer endocrine resistance. In this study, we conducted an integrated analysis combining scRNA-seq and scATAC-seq on more than 80,000 breast tissue cells from two normal tissues (NTs), three primary tumors (PTs), and three tamoxifen-treated recurrent tumors (RTs). A variety of cell types among breast tumor tissues were identified, PT- and RT-specific cancer cell states (CSs) were defined, and a heterogeneity-guided core signature (HCS) was derived through such integrated analysis. Functional experiments were performed to validate the oncogenic role of BMP7, a key gene within the core signature. We observed a striking level of cell-to-cell heterogeneity among six tumor tissues and delineated the primary to recurrent tumor progression, underscoring the significance of these single-cell level tumor cell clusters classified from scRNA-seq data. We defined nine CSs, including five PT-specific, three RT-specific, and one PT-RT-shared CSs, and identified distinct open chromatin regions of CSs, as well as a HCS of 137 genes. In addition, we predicted specific transcription factors (TFs) associated with the core signature and novel biological/metabolism pathways that mediate the communications between CSs and the tumor microenvironment (TME). We finally demonstrated that BMP7 plays an oncogenic role in tamoxifen-resistant breast cancer cells through modulating MAPK signaling pathways. Our integrated single-cell analysis provides a comprehensive understanding of the tumor heterogeneity in tamoxifen resistance. We envision this integrated single-cell epigenomic and transcriptomic measure will become a powerful approach to unravel how epigenetic factors and the tumor microenvironment govern the development of tumor heterogeneity and to uncover potential therapeutic targets that circumvent heterogeneity-related failures.

NSC95397 Is a Novel HIV-1 Latency-Reversing Agent.

The latent viral reservoir represents one of the major barriers to curing HIV-1. Focus on the "kick and kill" (also called "shock and kill") approach, in which virus expression is reactivated, and then cells producing virus are selectively depleted, has led to the discovery of many latency-reversing agents (LRAs) that have furthered our understanding of the mechanisms driving HIV-1 latency and latency reversal. Thus far, individual compounds have yet to be robust enough to work as a therapy, highlighting the importance of identifying new compounds that target novel pathways and synergize with known LRAs. In this study, we identified a promising LRA, NSC95397, from a screen of ~4250 compounds. We validated that NSC95397 reactivates latent viral transcription and protein expression from cells with unique integration events and across different latency models. Co-treating cells with NSC95397 and known LRAs demonstrated that NSC95397 synergizes with different drugs under both standard normoxic and physiological hypoxic conditions. NSC95397 does not globally increase open chromatin, and bulk RNA sequencing revealed that NSC95397 does not greatly increase cellular transcription. Instead, NSC95397 downregulates pathways key to metabolism, cell growth, and DNA repair-highlighting the potential of these pathways in regulating HIV-1 latency. Overall, we identified NSC95397 as a novel LRA that does not largely alter global transcription, shows potential for synergy with known LRAs, and may act through novel pathways not previously recognized for their ability to modulate HIV-1 latency.

Mechanisms driving epigenetic and transcriptional responses of microglia in a neurodegenerative lysosomal storage disorder model.

Lysosomal dysfunction is causally linked to neurodegeneration in many lysosomal storage disorders (LSDs) and is associated with various age-related neurodegenerative diseases 1,2 , but there is limited understanding of the mechanisms by which altered lysosomal function leads to changes in gene expression that drive pathogenic cellular phenotypes. To investigate this question, we performed systematic imaging, transcriptomic, and epigenetic studies of major brain cell types in Sgsh null (KO) mice, a preclinical mouse model for Sanfilippo syndrome (Mucopolysaccharidosis Type IIIA, MPS-IIIA) 3,4 . MPS-IIIA is a neurodegenerative LSD caused by homozygous loss-of-function (LoF) mutations in SGSH which results in severe early-onset developmental, behavioral, and neurocognitive impairment 5-15 . Electron microscopy, immunohistochemistry, and single-nucleus RNA-sequencing analysis revealed microglia as the cell type exhibiting the most dramatic phenotypic alterations in Sgsh KO mice. Further temporal analysis of microglia gene expression showed dysregulation of genes associated with lysosomal function and immune signaling pathways beginning early in the course of the disease. Sgsh deficiency similarly resulted in increases in open chromatin and histone acetylation at thousands of putative microglia-specific enhancers associated with upregulated genes but had much less impact on the epigenetic landscapes of neurons or oligodendrocytes. We provide evidence for dominant and context-dependent roles of members of the MITF/TFE family as major drivers of microglia-specific epigenetic and transcriptional changes resulting from lysosomal stress that are dependent on collaborative interactions with PU.1/ETS and C/EBP transcription factors. Lastly, we show that features of the transcriptomic and epigenetic alterations observed in murine Sgsh deficiency are also observed in microglia derived from mouse models of age-related neurodegeneration and in human Alzheimer's disease patients, revealing common and disease-specific transcriptional mechanisms associated with disease-associated microglia phenotypes.

Deciphering the tissue-specific functional effect of Alzheimer risk SNPs with deep genome annotation

Alzheimer’s disease (AD) is a highly heritable brain dementia, along with substantial failure of cognitive function. Large-scale genome-wide association studies (GWASs) have led to a set of SNPs significantly associated with AD and related traits. GWAS hits usually emerge as clusters where a lead SNP with the highest significance is surrounded by other less significant neighboring SNPs. Although functionality is not guaranteed even with the strongest associations in GWASs, lead SNPs have historically been the focus of the field, with the remaining associations inferred to be redundant. Recent deep genome annotation tools enable the prediction of function from a segment of a DNA sequence with significantly improved precision, which allows in-silico mutagenesis to interrogate the functional effect of SNP alleles. In this project, we explored the impact of top AD GWAS hits around APOE region on chromatin functions and whether it will be altered by the genetic context (i.e., alleles of neighboring SNPs). Our results showed that highly correlated SNPs in the same LD block could have distinct impacts on downstream functions. Although some GWAS lead SNPs showed dominant functional effects regardless of the neighborhood SNP alleles, several other SNPs did exhibit enhanced loss or gain of function under certain genetic contexts, suggesting potential additional information hidden in the LD blocks.

Evolution of 3D Chromatin Folding.

Studies examining the evolution of genomes have focused mainly on sequence conservation. However, the inner working of a cell implies tightly regulated crosstalk between complex gene networks controlled by small dispersed regulatory elements of physically contacting DNA regions. How these different levels of chromatin organization crosstalk in different species underpins the potential for genome evolutionary plasticity. We review the evolution of chromatin organization across the Animal Tree of Life. We introduce general aspects of the mode and tempo of genome evolution to later explore the multiple layers of genome organization. We argue that both genome and chromosome size modulate patterns of chromatin folding and that chromatin interactions facilitate the formation of lineage-specific chromosomal reorganizations, especially in germ cells. Overall, analyzing the mechanistic forces involved in the maintenance of chromatin structure and function of the germ line is critical for understanding genome evolution, maintenance, and inheritance.

KDM4A facilitates neuropathic pain and microglial M1 polarization by regulating BDNF in a rat model of brachial plexus avulsion

BackgroundMany patients with brachial plexus avulsion (BPA) suffer from neuropathic pain, but the mechanism remains elusive. Modifications of histones, the proteins responsible for organizing DNA, may play an important role...

EPCO-26. LONGITUDINAL SINGLE-CELL ANALYSES IDENTIFY DRIVERS OF GENETIC, EPIGENOMIC, AND CELLULAR EVOLUTION IN IDH-MUTANT GLIOMA

Abstract Despite standard treatment, IDH-mutant gliomas inevitably recur. Previous studies suggest that therapeutic resistance may result from a combination of intratumoral cellular heterogeneity, epigenetic evolution, and acquired genetic alterations. However, how these multilayered molecular features interact to influence the evolutionary paths of IDH-mutant gliomas remains incompletely understood. To chart this evolution, 75 glioma samples were longitudinally collected from 35 patients (n = 13 oligodendroglioma, n = 22 astrocytoma) and profiled using single nucleus RNA sequencing, ATAC sequencing, and bulk DNA sequencing. Analysis of 331,016 nuclei (snRNA) revealed eight tumor microenvironment cell types and five pan-IDH-mutant malignant cellular states. The malignant states were distributed along a cellular hierarchy of stem-like populations (neural progenitor cell-like, oligodendrocyte progenitor cell-like, and undifferentiated) and more differentiated populations (mesenchymal-like and astrocyte-like) with cycling cells enriched in the stem-like populations. Joint open chromatin accessibility data (snATAC) was available for a subset of malignant cells (71,088 nuclei) and state-specific differentially accessible peaks supported that these states are epigenetically encoded. Across both subtypes, higher tumor grade was associated with reduced malignant cell differentiation and increased cycling populations. Between the two time points, we identified that there was a longitudinal increase in the cycling, undifferentiated, and mesenchymal-like populations with a corresponding decrease in the astrocyte-like population. Longitudinal genetic analysis revealed that 19 of 35 tumors acquired at least one of the following genetic alterations: treatment-associated hypermutation, cell cycle alteration, or large changes in copy number alterations. Importantly, tumors that acquired these key genetic alterations demonstrated significant shifts towards reduced differentiation and increased cycling populations while those tumors that did not had more stable malignant profiles. Collectively, our results suggest a common cellular hierarchy across IDH-mutant gliomas with a shift towards reduced differentiation, increased cycling populations during disease progression that is driven by acquired genetic alterations.

Aryl-hydrocarbon receptor in smooth muscle cells protect against dioxin induced adverse remodeling of atherosclerosis.

Environmental exposure to dioxin has been linked to increased myocardial infarction. Smooth muscle cells (SMC) in the coronary vasculature play a critical role in atherosclerotic plaque remodeling due to their phenotypic plasticity, however, the detailed mechanism linking dioxin exposure to adverse SMC modulation is not well understood. Single-cell RNA and ATAC sequencing and histological analyses were performed on the aorta from mouse models of atherosclerosis exposed to 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) or control. Primary human coronary artery SMC (HCASMC) treated in culture with TCDD were used to perform RNA-Seq, ATAC-Seq, and functional phenotypic assays. ChIP-Seq was performed with antibodies against Aryl-hydrocarbon receptor (AHR) and TCF21, two of known SMC modulating transcription factors. Modulated SMC were the most transcriptionally responsive cell type to dioxin in the atherosclerotic aorta. Dioxin accelerated disease phenotype by promoting a modulated SMC phenotype early, resulting in increased lesion size, migration of SMC, and macrophage recruitment to the lesion. We found C3 expressing modulated SMCs to be likely contributing to the increased macrophage recruitment and inflammation. Analysis of the RNA-Seq data from HCASMC treated with TCDD showed differential enrichment of biological pathways related to cell migration, localization, and inflammation. Furthermore, ATAC-Seq data showed a significant activation for pathways regulating vascular development, cell migration, inflammation, and apoptosis. With TCDD treatment, there was also enrichment of AHR ChIP-Seq peaks, while the TCF21 enrichment decreased significantly. The SMC-specific Ahr knockout resulted in increased oxidative stress in SMC, increased lesion size and macrophage content, and loss of SMC lineage cells in the lesion cap when exposed to TCDD, consistent with a more vulnerable plaque phenotype. Dioxin adversely remodels atherosclerotic plaque by accelerating the SMC- phenotypic modulation, and increasing inflammation and oxidative stress resulting in increased macrophage recruitment and lesion size. Dioxin may adversely affect the SMC phenotype and disease state by affecting the TCF21 occupancy in the open chromatin regions. Furthermore, we observed that SMC-specific deletion of Ahr in mice resulted in worsening of dioxin mediated SMC modulation and atherosclerosis, suggesting that Ahr in SMC confers protection against dioxin by promoting a stable plaque phenotype and reducing dioxin induced oxidative stress. Exposure to dioxin, an environmental pollutant present in tobacco smoke and air pollution, accelerates smooth muscle cell modulation, and atherosclerosis.Dioxin exposure leads to inflammatory smooth muscle cell phenotype characterized by complement pathway activation and increased macrophage recruitment to plaqueAryl-hydrocarbon receptor in SMC protects against oxidative stress, and promotes a stable plaque phenotype.