Histone Marks H3K27ac Research Articles

A functional variant rs912304 for late-onset T1D risk contributes to islet dysfunction by regulating proinflammatory cytokine-responsive gene STXBP6 expression

BackgroundOur previous genome‑wide association studies (GWAS) have suggested rs912304 in 14q12 as a suggestive risk variant for type 1 diabetes (T1D). However, the association between this risk region and T1D subgroups and related clinical risk features, the underlying causal functional variant(s), putative candidate gene(s), and related mechanisms are yet to be elucidated.MethodsWe assessed the association between variant rs912304 and T1D, as well as islet autoimmunity and islet function, stratified by the diagnosed age of 12. We used epigenome bioinformatics analyses, dual luciferase reporter assays, and expression quantitative trait loci (eQTL) analyses to prioritize the most likely functional variant and potential causal gene. We also performed functional experiments to evaluate the role of the causal gene on islet function and its related mechanisms.ResultsWe identified rs912304 as a risk variant for T1D subgroups with diagnosed age ≥ 12 but not < 12. This variant is associated with residual islet function but not islet-specific autoantibody positivity in T1D individuals. Bioinformatics analysis indicated that rs912304 is a functional variant exhibiting spatial overlaps with enhancer active histone marks (H3K27ac and H3K4me1) and open chromatin status (ATAC-seq) in the human pancreas and islet tissues. Luciferase reporter gene assays and eQTL analyses demonstrated that the biallelic sites of rs912304 had differential allele-specific enhancer activity in beta cell lines and regulated STXBP6 expression, which was defined as the most putative causal gene based on Open Targets Genetics, GTEx v8 and Tiger database. Moreover, Stxbp6 was upregulated by T1D-related proinflammatory cytokines but not high glucose/fat. Notably, Stxbp6 over-expressed INS-1E cells exhibited decreasing insulin secretion and increasing cell apoptosis through Glut1 and Gadd45β, respectively.ConclusionsThis study expanded the genomic landscape regarding late-onset T1D risk and supported islet function mechanistically connected to T1D pathogenesis.

LAP2alpha facilitates myogenic gene expression by preventing nucleoplasmic lamin A/C from spreading to active chromatin regions.

A-type lamins form a filamentous meshwork beneath the nuclear membrane that anchors large heterochromatic genomic regions at the nuclear periphery. A-type lamins also exist as a dynamic, non-filamentous pool in the nuclear interior, where they interact with lamin-associated polypeptide 2 alpha (LAP2α). Both proteins associate with largely overlapping euchromatic genomic regions in the nucleoplasm, but the functional significance of this interaction is poorly understood. Here, we report that LAP2α relocates towards regions containing myogenic genes in the early stages of muscle differentiation, possibly facilitating efficient gene regulation, while lamins A and C mostly associate with genomic regions away from these genes. Strikingly, upon depletion of LAP2α, A-type lamins spread across active chromatin and accumulate at regions of active H3K27ac and H3K4me3 histone marks in the vicinity of myogenic genes whose expression is impaired in the absence of LAP2α. Reorganization of A-type lamins on chromatin is accompanied by depletion of the active chromatin mark H3K27ac and a significantly impaired myogenic differentiation. Thus, the interplay of LAP2α and A-type lamins is crucial for proper positioning of intranuclear lamin A/C on chromatin to allow efficient myogenic differentiation.

Evolutionarily conserved enhancer-associated features within the MYEOV locus suggest a regulatory role for this non-coding DNA region in cancer.

The myeloma overexpressed gene (MYEOV) has been proposed to be a proto-oncogene due to high RNA transcript levels found in multiple cancers, including myeloma, breast, lung, pancreas and esophageal cancer. The presence of an open reading frame (ORF) in humans and other primates suggests protein-coding potential. Yet, we still lack evidence of a functional MYEOV protein. It remains undetermined how MYEOV overexpression affects cancerous tissues. In this work, we show that MYEOV has likely originated and may still function as an enhancer, regulating CCND1 and LTO1. Firstly, MYEOV 3' enhancer activity was confirmed in humans using publicly available ATAC-STARR-seq data, performed on B-cell-derived GM12878 cells. We detected enhancer histone marks H3K4me1 and H3K27ac overlapping MYEOV in multiple healthy human tissues, which include B cells, liver and lung tissue. The analysis of 3D genome datasets revealed chromatin interactions between a MYEOV-3'-putative enhancer and the proto-oncogene CCND1. BLAST searches and multi-sequence alignment results showed that DNA sequence from this human enhancer element is conserved from the amphibians/amniotes divergence, with a 273bp conserved region also found in all mammals, and even in chickens, where it is consistently located near the corresponding CCND1 orthologues. Furthermore, we observed conservation of an active enhancer state in the MYEOV orthologues of four non-human primates, dogs, rats, and mice. When studying this homologous region in mice, where the ORF of MYEOV is absent, we not only observed an enhancer chromatin state but also found interactions between the mouse enhancer homolog and Ccnd1 using 3D-genome interaction data. This is similar to the interaction observed in humans and, interestingly, coincides with CTCF binding sites in both species. Taken together, this suggests that MYEOV is a primate-specific gene with a de novo ORF that originated at an evolutionarily older enhancer region. This deeply conserved putative enhancer element could regulate CCND1 in both humans and mice, opening the possibility of studying MYEOV regulatory functions in cancer using non-primate animal models.

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo.

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues. Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~ 50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3). HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

ASH1L guards cis-regulatory elements against cyclobutane pyrimidine dimer induction.

The histone methyltransferase ASH1L, first discovered for its role in transcription, has been shown to accelerate the removal of ultraviolet (UV) light-induced cyclobutane pyrimidine dimers (CPDs) by nucleotide excision repair. Previous reports demonstrated that CPD excision is most efficient at transcriptional regulatory elements, including enhancers, relative to other genomic sites. Therefore, we analyzed DNA damage maps in ASH1L-proficient and ASH1L-deficient cells to understand how ASH1L controls enhancer stability. This comparison showed that ASH1L protects enhancer sequences against the induction of CPDs besides stimulating repair activity. ASH1L reduces CPD formation at C-containing but not at TT dinucleotides, and no protection occurs against pyrimidine-(6,4)-pyrimidone photoproducts or cisplatin crosslinks. The diminished CPD induction extends to gene promoters but excludes retrotransposons. This guardian role against CPDs in regulatory elements is associated with the presence of H3K4me3 and H3K27ac histone marks, which are known to interact with the PHD and BRD motifs of ASH1L, respectively. Molecular dynamics simulations identified a DNA-binding AT hook of ASH1L that alters the distance and dihedral angle between neighboring C nucleotides to disfavor dimerization. The loss of this protection results in a higher frequency of C->T transitions at enhancers of skin cancers carrying ASH1L mutations compared to ASH1L-intact counterparts.

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo.

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues. Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3). HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

KLF2 enhancer variant rs4808485 increases lupus risk by modulating inflammasome machinery and cellular homoeostasis

ObjectiveA recent genome-wide association study linked KLF2 as a novel Asian-specific locus for systemic lupus erythematosus (SLE) susceptibility. However, the underlying causal functional variant(s), cognate target gene(s) and genetic mechanisms...

F.2 Comprehensively mapping transcriptionally relevant histone modifications in aggressive meningioma leads to novel biologic insights and therapeutic vulnerabilities

Background: We recently identified four molecular subgroups of meningioma with distinct biology and outcomes. While two (MG3/MG4) are associated with poor outcome, they display divergent transcriptional profiles (enriched in metabolic and cell cycling pathways, respectively) and therapeutic vulnerabilities (MG3 has no clear treatment target). We sought to understand drivers of these key differences at a chromatin level. Methods: We profiled MG3/MG4 meningiomas for common histone marks H3K27me3, H3K27Ac, H3K4me1, H3K4me3, H3K9me3, and H3K36me3. Multiple computational approaches were used to compare MG3 and MG4 tumours including superenhancer ranking, differential binding analysis, and unsupervised clustering. Results: Our cohort includes 11-20 meningiomas per histone mark. Clustering revealed striking separation of subgroups based on multiple histone marks, particularly H3K36me3. FOXC1, a known driver of the epithelial to mesenchymal transition, was identified as a recurrent superenhancer in both groups, whereas MG3-specific superenhancers mapped to immune regulatory networks. Integrated differential binding analysis confirmed an immune-rich microenvironment in MG3 tumours driven by multiple histone marks, suggesting a role for targeting novel immune checkpoint genes CD84 and CD48. Conclusions: This study is the first to apply integrated analysis of multiple histone modifications to aggressive meningioma. We further characterize MG3 tumours by identifying an epigenetically-driven immune phenotype and propose novel treatment targets.

Human cytomegalovirus infection coopts chromatin organization to diminish TEAD1 transcription factor activity.

Human cytomegalovirus (HCMV) infects up to 80% of the world's population. Here, we show that HCMV infection leads to widespread changes in human chromatin accessibility and chromatin looping, with hundreds of thousands of genomic regions affected 48 hours after infection. Integrative analyses reveal HCMV-induced perturbation of Hippo signaling through drastic reduction of TEAD1 transcription factor activity. We confirm extensive concordant loss of TEAD1 binding, active H3K27ac histone marks, and chromatin looping interactions upon infection. Our data position TEAD1 at the top of a hierarchy involving multiple altered important developmental pathways. HCMV infection reduces TEAD1 activity through four distinct mechanisms: closing of TEAD1-bound chromatin, reduction of YAP1 and phosphorylated YAP1 levels, reduction of TEAD1 transcript and protein levels, and alteration of TEAD1 exon-6 usage. Altered TEAD1-based mechanisms are highly enriched at genetic risk loci associated with eye and ear development, providing mechanistic insight into HCMV's established roles in these processes.

Benchmarking tools for transcription factor prioritization

Spatiotemporal regulation of gene expression is controlled by transcription factor (TF) binding to regulatory elements, resulting in a plethora of cell types and cell states from the same genetic information. Due to the importance of regulatory elements, various sequencing methods have been developed to localise them in genomes, for example using ChIP-seq profiling of the histone mark H3K27ac that marks active regulatory regions. Moreover, multiple tools have been developed to predict TF binding to these regulatory elements based on DNA sequence. As altered gene expression is a hallmark of disease phenotypes, identifying TFs driving such gene expression programs is critical for the identification of novel drug targets.In this study, we curated 84 chromatin profiling experiments (H3K27ac ChIP-seq) where TFs were perturbed through e.g., genetic knockout or overexpression. We ran nine published tools to prioritize TFs using these real-world datasets and evaluated the performance of the methods in identifying the perturbed TFs. This allowed the nomination of three frontrunner tools, namely RcisTarget, MEIRLOP and monaLisa. Our analyses revealed opportunities and commonalities of tools that will help to guide further improvements and developments in the field.

Abstract 7555: Selective inhibition of the CoREST complex as a novel therapeutic target of atypical teratoid rhabdoid tumors

Abstract Atypical teratoid/rhabdoid tumors (AT/RT) are deadly infantile brain tumors, primarily arising from a single recurring biallelic loss-of-function mutation in SMARCB1 which is a member of the SWI-SNF chromatin remodeling complex. The partially functional SWI-SNF complex (without SMARCB1) is able to inhibit EZH2 at most locations of the genome except at gene promoter regions where EZH2 co-localizes with the REST complex. This results in hypermethylation and subsequent repression of genes coding for neuronal differentiation and tumor suppressors. CoREST is an important transcriptional repressor and consists of LSD1 and HDAC1/2 bound together by the scaffolding protein RCOR. Increased activity of EZH2 and the CoREST complex may help AT/RT cells maintain a stem cell-like state, driving its aggressive growth and therapy resistance. We hypothesize that functional inhibition of the CoREST complex by targeting LSD1 and HDAC1/2 might help restore the epigenetic balance lost in AT/RT cells.We selectively inhibited the CoREST complex in AT/RT cells (BT37, CHLA05 and CHLA06) by using a bi-functional LSD1 and HDAC1/2 inhibitor, Corin, which has previously been shown to be effective against melanoma and diffuse midline glioma. Administration of corin to AT/RT cells significantly increased H3K4me1, H3K9Ac and H3K27Ac histone marks (p &amp;lt; 0.001, Western blotting). Corin inhibited AT/RT cell growth and proliferation (MUSE Viacount p&amp;lt;0.05, BrdU immunofluorescence compared to control), and increased apoptosis (cPARP Western blotting, CC3 immunofluorescence, MUSE Annexin V assay, p&amp;lt;0.05 compared to control). Corin appears to be specific in targeting AT/RT as it did not affect brain organoids derived from iPSCs (cPARP Western blotting). Corin induced neuronal differentiation in AT/RT as seen by the increased expression of Beta-3-tubulin, MAP2 and Synaptophysin (Immunofluorescence). Enrichment analysis of differentially expressed genes after whole mRNA sequencing of AT/RT cells treated with corin for 24 h revealed upregulation of pathways related to neuronal differentiation and axonal development and downregulation of pathways related to proliferation and cell cycle (pFDR &amp;lt; 0.05, FC&amp;gt;2). In mice bearing orthotopic CHLA06 tumors, intratumoral delivery of a single dose of corin for 24 h increased H3K27me3 and H3K27ac, demonstrating that corin engaged its target in vivo, reduced growth and induced apoptosis (Ki67, cPARP). Survival studies show a trend toward single agent activity against the very aggressive CHLA06 orthotopic xenograft (increased in median survival from 14 to 35 days compared to control). These studies demonstrate that targeting the CoREST complex disrupts AT/RT epigenetic abnormalities, reversing their stem cell-like state, slowing tumor cell growth, and inducing apoptosis, that may translate into the clinical setting to help improve AT/RT’s dismal survival. Citation Format: Anupa Geethadevi, Marianne Collard, Tyler Findlay, Yiming Deng, Robert Fisher, Nikhil Vaidya, Charles Eberhart, Philip A. Cole, Rhoda Alani, Jeffrey Rubens, Eric Raabe. Selective inhibition of the CoREST complex as a novel therapeutic target of atypical teratoid rhabdoid tumors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7555.

Abstract 2861: Rhabdomyosarcoma fusion-oncogene alters the chromatin landscape to initially drive a neural signature in vivo

Abstract Fusion-positive rhabdomyosarcoma is an aggressive pediatric cancer lacking curative therapies, and outcomes for children with this disease have not improved in decades. The chimeric transcription factor PAX3::FOXO1 is the most common and lethal driver of fusion-positive rhabdomyosarcoma, and it consists of the DNA binding domains of PAX3 fused to the transactivation domain of FOXO1. Despite its well-established essentiality for tumorigenesis, how PAX3::FOXO1 initially generates a tumorigenic cell state has been challenging to study due to the lack of in vivo models expressing the fusion-oncogene at accessible early developmental time points. Here, we developed a novel zebrafish mRNA injection model that ubiquitously expresses PAX3::FOXO1 during gastrulation representing all three germ layers. Using this high-throughput model, we used ChIP-seq to identify initial in vivo PAX3::FOXO1 binding sites and targets. PAX3::FOXO1 has pioneering activity, where it can utilize partial-motif recognition to bind to inaccessible chromatin and alter chromatin structure. With biochemical and 2D chromatin sequencing approaches in our model, we find evidence for in vivo PAX3::FOXO1 pioneering activity. PAX3::FOXO1 consists of two DNA binding domains, a paired and homeobox domain, and we demonstrate a new mode of pioneering activity driven by its homeobox domain given a strong enrichment of homeobox-related motifs at PAX3::FOXO1 binding sites. PAX3::FOXO1 activity in developing zebrafish embryos results in an arrested development phenotype and transcriptional signatures enriched for pathways such as an inhibition of segmentation and myogenic development. Critically, we find that PAX3::FOXO1 directly activates neural-related gene targets by increasing chromatin accessibility and re-distribution of the active histone mark H3K27Ac to these loci. These discoveries are particularly striking because rhabdomyosarcoma is traditionally associated with skeletal muscle characteristics. With our findings, we hypothesize neural gene activation is a critical mechanism for fusion-positive rhabdomyosarcoma tumor initiation and could be a critical pathway across various stages of rhabdomyosarcoma tumorigenesis. This hypothesis agrees with our previous discovery of a novel PAX3::FOXO1 cooperating neural transcription factor, HES3, which we find is directly activated by PAX3::FOXO1 in our mRNA injection model. Our long-term goal is to utilize this versatile model to functionally evaluate cooperation between PAX3::FOXO1 and both canonical rhabdomyosarcoma and neural transcription factors. Our approach will employ a cross-species comparative analysis across this injection model, our zebrafish tumor models, other model systems, and patient data. This work will identify conserved targets and mechanisms of PAX3::FOXO1 tumorigenesis that could become new therapeutic vulnerabilities. Citation Format: Jack Kucinski, Cenny Taslim, Alexi Tallan, Matthew Cannon, Katherine Silvius, Benjamin Stanton, Genevieve Kendall. Rhabdomyosarcoma fusion-oncogene alters the chromatin landscape to initially drive a neural signature in vivo [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2861.

Deubiquitylase ubiquitin-specific protease 7 plays a crucial role in the lineage differentiation of preimplantation blastocysts†.

Preimplantation embryos undergo a series of important biological events, including epigenetic reprogramming and lineage differentiation, and the key genes and specific mechanisms that regulate these events are critical to reproductive success. Ubiquitin-specific protease 7 (USP7) is a deubiquitinase involved in the regulation of a variety of cellular functions, yet its precise function and mechanism in preimplantation embryonic development remain unknown. Our results showed that RNAi-mediated silencing of USP7 in mouse embryos or treatment with P5091, a small molecule inhibitor of USP7, significantly reduced blastocyst rate and blastocyst quality, and decreased total and trophectoderm cell numbers per blastocyst, as well as destroyed normal lineage differentiation. The results of single-cell RNA-seq, reverse transcription-quantitative polymerase chain reaction, western blot, and immunofluorescence staining indicated that interference with USP7 caused failure of the morula-to-blastocyst transition and was accompanied by abnormal expression of key genes (Cdx2, Oct4, Nanog, Sox2) for lineage differentiation, decreased transcript levels, increased global DNA methylation, elevated repressive histone marks (H3K27me3), and decreased active histone marks (H3K4me3 and H3K27ac). Notably, USP7 may regulate the transition from the morula to blastocyst by stabilizing the target protein YAP through the ubiquitin-proteasome pathway. In conclusion, our results suggest that USP7 may play a crucial role in preimplantation embryonic development by regulating lineage differentiation and key epigenetic modifications.

Lithocholic Acid Alleviates Deoxynivalenol-Induced Inflammation and Oxidative Stress via PPARγ-Mediated Epigenetically Transcriptional Reprogramming in Porcine Intestinal Epithelial Cells.

Deoxynivalenol (DON) is a common mycotoxin that induces intestinal inflammation and oxidative damage in humans and animals. Given that lithocholic acid (LCA) has been suggested to inhibit intestinal inflammation, we aimed to investigate the protective effects of LCA on DON-exposed porcine intestinal epithelial IPI-2I cells and the underlying mechanisms. Indeed, LCA rescued DON-induced cell death in IPI-2I cells and reduced DON-stimulated inflammatory cytokine levels and oxidative stress. Importantly, the nuclear receptor PPARγ was identified as a key transcriptional factor involved in the DON-induced inflammation and oxidative stress processes in IPI-2I cells. The PPARγ function was found compromised, likely due to the hyperphosphorylation of the p38 and ERK signaling pathways. In contrast, the DON-induced inflammatory responses and oxidative stress were restrained by LCA via PPARγ-mediated reprogramming of the core inflammatory and antioxidant genes. Notably, the PPARγ-modulated transcriptional regulations could be attributed to the altered recruitments of coactivator SRC-1/3 and corepressor NCOR1/2, along with the modified histone marks H3K27ac and H3K18la. This study emphasizes the protective actions of LCA on DON-induced inflammatory damage and oxidative stress in intestinal epithelial cells via PPARγ-mediated epigenetically transcriptional reprogramming, including histone acetylation and lactylation.

Immune-related transcriptomic and epigenetic reconfiguration in BV2 cells after lipopolysaccharide exposure: an in vitro omics integrative study.

Molecular alterations affecting microglia have been consistently associated with the inflammatory response. These cells can have pro- or anti-inflammatory activity, phenotypes thought to be regulated by epigenetic mechanisms. Still, little is known about the details on how epigenetic marks regulate the expression of genes in the context of an inflammatory response. Through CUT&RUN, we profiled four genome-wide histone marks (HM) (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in lipopolysaccharide-exposed cells and compared their distributions to control cells. Transcriptomic profiles were determined through RNA-seq and differentially expressed genes were identified and contrasted with the epigenetic landscapes. Other downstream analyses were also included in this study. Our results illustrate an effectively induced M1 phenotype in microglial cells derived from LPS exposure. We observed differential bound regions associated with the genes classically involved in the inflammatory response in the expected direction according to each histone modification. Consistently, our transcriptomic analysis yielded a conspicuous illustration of the LPS-induced immune activity showing the up-regulation of Nf-κB-induced mRNAs (TNF-α, nfκbiz, nfκbia) and other important genes (Marco, Il-6, etc.). Furthermore, we integrated both omics profiles and identified an important reconfiguration of the genome induced by LPS. The latter was depicted by 8 different chromatin states that changed between conditions and that associated with unique clusters of differentially expressed genes, which not only represented regulatory elements, but also underlined distinct biological functions (inhibition of morphogenesis; changes in metabolism, homeostasis, and cytokine regulation; activation of the inflammatory response). This study exhibits important differences in the distribution of histone modifications in treated and control BV2 cells, constituting an epigenetic reconfiguration that leads to the inflammatory response. Also, it highlights the importance of these marks' regulatory role in gene expression and provides possible targets for further studies in the context of inflammation.

EPCO-46. CHARACTERIZING ATYPICAL TERATOID RHABDOID TUMOR SUBTYPE ENHANCER INTERACTION LANDSCAPE

Abstract Atypical Teratoid Rhabdoid Tumor (ATRT) is one of the most aggressive pediatric tumors of the central nervous system (CNS). The majority of ATRTs are characterized by the recurrent bi-allelic loss of the SWI/SNF nucleosome remodeling complex core component gene, SMARCB1, or SMARCA4 in rare cases. Despite its bland genome, three ATRT subgroups “SHH”, “TYR”, and “MYC” were discovered with unique molecular profiles. Also, distinct activating histone mark H3K27ac profiles between molecular subgroups were identified, suggesting subgroup-specific regulatory networks through enhancers. However, comprehensive cis-regulatory network descriptions in ATRT are lacking. Here, we aim to elucidate ATRT subgroup pathogenesis through its dysregulated enhancer-promoter interactions. We have created a regulatory mutational profile of the three ATRT subgroups using whole-genome sequencing method of 60 ATRT tumors and matching genome-wide H3K27ac profiles on 16 tumors. To investigate mutation burden on cis-regulations, we have employed high-throughput Hi-C chromatin immunoprecipitation (HiChIP) technique to directly map enhancer sequences involved in DNA-looping with high specificity and sensitivity. Recent findings from our lab suggest that activating interaction PDGFRB confer to therapeutics venerability in the MYC subgroup. Here, our preliminary analysis shows a hub of interactions near PDGFRB reveals an extensive interaction hub with long-range enhancer interactions extending as far as 1500Kb, suggesting a much more complicated regulatory picture than previously expected. Interestingly, we have identified long-range DNA interactions between PDGFRB and H3K27ac marks upstream of its gene body and overstretching into other gene elements. This result suggests a potential alternative activation mechanism for MYC cell lines. Our next steps are to perform performing additional HiChIP assays on other ATRT cell lines from other subgroups to capture the full range of ATRT’s heterogeneity. Comprehensive mapping of the non-coding genome and enhancer regulatory network would help identify critical pathogenic elements underlying ATRT biology and provide insights into therapeutic targets.

53 Establishing a Pan-Epigenome for Cattle and Sheep

Abstract Genetic and epigenetic factors contribute towards phenotypic variation in complex traits. Epigenetic modifications can affect transcription and translation and have been shown to influence variation in economically important phenotypes in agricultural species. The onset of pangenomics in livestock species is disconnected from pan-epigenomics and hence lacks a fully comprehensive view of diversity in livestock species. To establish a pan-epigenome for cattle and sheep we are comprehensively characterizing epigenomic conservation and diversity through examination of 5-methylcytosine, 5-hydroxymethylcytosine, histone marks H3K4me3, H3K27ac, H3K4me1 and H3K27me3, open chromatin, chromatin accessibility and n6-methyladenosine in the liver, muscle and prefrontal cortex of Brahman, Angus and Holstein cattle as well as Suffolk, Rambouillet and Katahdin sheep. Finally, comparative epigenomics is being used to identify orthologous regions of the sheep and bovine epigenomes that harbor conserved epigenetic modifications and variants. The establishment of a cattle and sheep pan-epigenome will increase our understanding of how different genomes, epigenomes, and gene products from diverse breeds of cattle and sheep affect a variety of important biological phenotypes, supporting more accurate prediction of traits and improving breeding strategies.

STAG2 loss Induces HSC Programs By Modulating Accessibility of AP-1 Bound Enhancers

Myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML) are blood disorders characterized by clonal expansion of mutant hematopoietic stem and progenitor cells (HSPCs). Mutations in genes encoding the cohesin complex, including STAG2, RAD21, SMC3 and SMC1A, are frequent genetic drivers in MDS and AML. However, the mechanisms by which cohesin mutations lead to clonal expansion of HSPCs are not well understood. Cohesin is essential for sister chromatid cohesion, DNA damage repair, maintenance of chromatin architecture, and gene regulation. To investigate the impact of cohesin mutations on chromatin accessibility, we performed ATAC-Seq in a panel of six isogenic STAG2 wild type (WT) and knockout (KO) U937 leukemia cell lines. Among the total of 97,798 ATAC-Seq peaks, we identified significantly more STAG2 KO (8,859) versus WT (1,622) uniquely accessible chromatin sites, which were strongly enriched for intergenic regions. Integration of ChIP-Seq profiling of cohesin binding and histone modifications demonstrated cohesin binding and an increase in the active enhancer and promoter histone mark H3K27Ac in the STAG2 KO newly accessible chromatin, suggesting a direct role for cohesin in regulating the activity of such enhancers. Together, our data demonstrate that STAG2 KO-induced open chromatin may be enriched for novel STAG2 KO-specific enhancers. To further infer the regulatory function of these putative enhancers, we predicted binding of transcription factors (TFs) in ATAC-Seq-defined regions by performing TF footprinting using the Transcription factor Occupancy Prediction By Investigation of ATAC-seq Signal (TOBIAS) algorithm. We observed that binding of members of the AP-1 transcription family, including JUN FOS andATF, previously implicated in oncogenic transformation, was significantly enriched in the STAG2 KO cells (Figure 1A). Conversely, STAG2 KO cells lost CEBP TF binding in regions that lost chromatin accessibility. To identify potential regulatory targets of these TFs, we mapped TF footprints to annotated promoter regions and to enhancer-gene pairs determined by the Activity by Contact (ABC) model. We observed a strong correlation between predicted AP-1 binding and gene expression of putative regulated genes in STAG2 KO cells, further supporting the role of these newly accessible AP-1 bound sites as regulatory elements. To validate our findings using in vivo models of cohesin-mutant disease, we performed RNA-Seq and ATAC-Seq in HSPC isolated from a Tet2/Stag2 mouse model of MDS previously developed in our lab. We confirmed an increase in accessible chromatin and demonstrated a strong enrichment of predicted Ap-1 binding in intergenic regions, in agreement with our observations in the U937 model. Gene set enrichment analysis in the Tet2/Stag2-mutant versus WT HSPC revealed that putative targets of Ap-1 binding were enriched for hematopoietic stem cell gene signatures (Figure 1B) . This is consistent with our observation of clonal expansion of HSPCs in this model, and suggests a mechanism by which novel Ap-1 bound enhancers may facilitate expression of key genes involved in HSPC expansion and development of MDS/AML. To explore whether our findings were applicable to other published models of Stag2 mutant myeloid neoplasms, we performed TF footprinting using publicly available ATAC-Seq datasets generated in HSPC isolated from Stag2 and Runx1/Stag2 conditional knockout mice. We recapitulated increased chromatin accessibility and significant enrichment of cohesin-bound Ap-1 footprints in Stag2 mutant over WT conditions. Together, these data further validate conservation of STAG2 KO-induced AP-1 bound regulatory elements in human and mouse models of MDS and AML. In summary, we show that STAG2 loss leads to the establishment of a novel set of accessible intergenic regions with features of enhancer activity in vitro and in vivo. We identify multiple members of the AP-1 complex with enriched predicted binding at these putative enhancers, and demonstrate that AP-1 may drive expression of HSPC signatures in models of Stag2 mutant MDS in vivo. We illustrate here a potential role for the AP-1 complex in cohesin mutant myeloid malignancies, which we are currently functionally interrogating.

Multitissue H3K27ac profiling of GTEx samples links epigenomic variation to disease

Genetic variants associated with complex traits are primarily noncoding, and their effects on gene-regulatory activity remain largely uncharacterized. To address this, we profile epigenomic variation of histone mark H3K27ac across 387 brain, heart, muscle and lung samples from Genotype-Tissue Expression (GTEx). We annotate 282 k active regulatory elements (AREs) with tissue-specific activity patterns. We identify 2,436 sex-biased AREs and 5,397 genetically influenced AREs associated with 130 k genetic variants (haQTLs) across tissues. We integrate genetic and epigenomic variation to provide mechanistic insights for disease-associated loci from 55 genome-wide association studies (GWAS), by revealing candidate tissues of action, driver SNPs and impacted AREs. Lastly, we build ARE–gene linking scores based on genetics (gLink scores) and demonstrate their unique ability to prioritize SNP–ARE–gene circuits. Overall, our epigenomic datasets, computational integration and mechanistic predictions provide valuable resources and important insights for understanding the molecular basis of human diseases/traits such as schizophrenia.

Novel method of isolating nuclei of human oligodendrocyte precursor cells reveals substantial developmental changes in gene expression and H3K27ac histone modification.

Oligodendrocyte precursor cells (OPCs) generate differentiated mature oligodendrocytes (MOs) during development. In adult brain, OPCs replenish MOs in adaptive plasticity, neurodegenerative disorders, and after trauma. The ability of OPCs to differentiate to MOs decreases with age and is compromised in disease. Here we explored the cell specific and age-dependent differences in gene expression and H3K27ac histone mark in these two cell types. H3K27ac is indicative of active promoters and enhancers. We developed a novel flow-cytometry-based approach to isolate OPC and MO nuclei from human postmortem brain and profiled gene expression and H3K27ac in adult and infant OPCs and MOs genome-wide. In adult brain, we detected extensive H3K27ac differences between the two cell types with high concordance between gene expression and epigenetic changes. Notably, the expression of genes that distinguish MOs from OPCs appears to be under a strong regulatory control by the H3K27ac modification in MOs but not in OPCs. Comparison of gene expression and H3K27ac between infants and adults uncovered numerous developmental changes in each cell type, which were linked to several biological processes, including cell proliferation and glutamate signaling. A striking example was a subset of histone genes that were highly active in infant samples but fully lost activity in adult brain. Our findings demonstrate a considerable rearrangement of the H3K27ac landscape that occurs during the differentiation of OPCs to MOs and during postnatal development of these cell types, which aligned with changes in gene expression. The uncovered regulatory changes justify further in-depth epigenetic studies of OPCs and MOs in development and disease.