Ördög, Tamás

Introduction: Truncating mutations in Additional Sex Combs-Like 1 (ASXL1) are associated with a high-risk disease phenotype in myeloid malignancies. In chronic myelomonocytic leukemia (CMML), truncating ASXL1 mutations are known to increase transcriptional activity of leukemic driver genes and have been associated with gene body hypermethylation. We interrogated the transcriptome and methylome of patients with ASXL1-mutant (MT) and -wildtype (WT) CMML using a multi-omics approach to test the hypothesis that gene expression is mediated through gene body methylation. Methods: Bone marrow mononuclear cells from patients with ASXL1 WT (n=8) and MT (n=8) CMML were subjected to targeted NGS of DNA, whole transcriptome shotgun sequencing (RNA-seq), and immunoprecipitation of DNA methyl residues (DIP-seq). After quality control all samples were sequenced on an Illumina HiSeq 4000 before further processing and data analysis. Differential gene expression analysis was performed to identify genes up-regulated in MT CMML. The samples in the two groups were treated as biological replicates and subjected to a consensus peak calling strategy requiring an overlap of at least 30% between samples and an adjusted p-value < 5x10-5 for a methylation peak to be considered statistically significant. For validation purposes methylation analysis was performed on 3 ASXL1 MT and 3 WT CMML patients using Illumina Infinium MethylationEPIC microarrays and differentially methylated regions were identified using a bump hunting strategy. Gene body methylation was defined as methylation in gene bodies (outside the promoter region, i.e. transcription start site ±2kb). Gene body methylation was compared between WT and MT CMML for the up-regulated genes and correlated with expression of all genes in MT CMML. Results: Sixteen WHO-defined CMML patients were included, median age 69 years (48-77), 63% male, 50% had truncating ASXL1 frame shift mutations. Abnormal karyotypes were observed in the same number of patients and the burden of co-mutations was similar between the two groups (median number per group 3 vs. 3, p=0.508). This included several modulators of DNA methylation including TET2, DNMT3A, and IDH2 (median number per group 1 vs. 1, p=0.699). There was a predominant up-regulation of gene expression in MT CMML: 707 genes up- and 124 down-regulated (FDR<0.050) without evidence of differential methylation in promoter regions of the differentially expressed genes. There was no difference in the number of genes with consensus methylation peaks in the gene body or the extent of gene body methylation (area under the consensus peaks) between MT and WT CMML for the differentially expressed genes (Figure 1a). We further validated these results using methylation microarrays (Figure 1b). Among the MT CMML patients, there was a curvilinear relationship between gene body methylation and gene expression across all genes with intermediate gene body methylation being associated with the highest gene expression (Figure 1c). As an alternative unbiased approach, we identified 1595 differentially methylated regions (DMR, regions with FDR<0.05) in the validation data set. We mapped all hypermethylated regions to gene bodies requiring that there was no concurrent hypermethylation of the promoter region of the same gene. With this approach we identified one of the 707 up-regulated genes (0.14%) with evidence of isolated gene body hypermethylation (HBZ, log2-fold change in gene expression 4.25, FDR=0.0008, DMR area=1.61, DMR FDR=0.018). Representative examples of up-regulated mitotic kinases without evidence of differential gene body methylation are shown in Figure 1d (pie charts represent the mean β-values of microarray probes, p-values represent the comparisons of mean β-values between MT and WT CMML per gene). Conclusions: Gene body methylation was positively associated with gene expression in MT CMML. However, the lack of differential gene body methylation between WT and MT CMML for the up-regulated genes make it an unlikely explanation for the observed increase in transcriptional activity among patients with MT CMML. Figure 1 Disclosures Ordog: Millipore Sigma: Patents & Royalties.

Introduction: Additional Sex Combs-Like 1 (ASXL1) is a chromatin modifier frequently affected by truncating mutations in myeloid malignancies. In chronic myelomonocytic leukemia (CMML), truncating ASXL1 mutations are associated with the overexpression of leukemogenic driver genes, however the molecular mechanisms underlying this transcriptional activation remain controversial. We performed single-cell chromatin accessibility studies of ASXL1-mutant (MT) and -wildtype (WT) CMML to identify cis-regulatory elements (CREs) such as distal enhancer elements that may explain the observed transcriptional activity in MT CMML. Methods: Bone marrow mononuclear cells from a patient with WT and a patient with MT CMML were subjected to single-cell DNA transposase accessibility studies (scATAC-seq). The cryopreserved cells were thawed and resuspended and approximately 100000 viable mononuclear cells per sample were subjected to transposase assays before proceeding to single-cell partitioning into gel beads in emulsion, barcoding, library construction, and sequencing following an established 10X Genomics workflow. The target cell recovery was approximately 2000 cells per sample. Genomic libraries were sequenced on an Illumina HiSeq 4000 before demultiplexing, alignment to the reference genome, and post-alignment quality control. The 10X Genomics Cell Ranger ATAC software was used for demultiplexing, alignment of the reads to the GRCh38 reference genome, filtering and quality control, counting of barcodes and unique molecular identifiers, identification of transposase cut sites, detection of accessible chromatin peaks, count matrix generation for peaks and transcription factors, clustering, and differential accessibility analysis. The two-sample Kolmogorov-Smirnov test for equality of distribution functions was used to compare the distributions of single-cell entropies. Single-cell cis-co-accessibility scores were estimated using Cicero. Results: To mitigate potential effects on the analysis introduced by unequal cell recovery due to experimental conditions, we randomly sampled 1500 single cells from each genotype. Visualizing the first two principal components of the chromatin accessibility data revealed the differences in chromatin conformation between genotypes and also the increased heterogeneity within the MT sample (Figure 1a). To quantify this increase in intratumoral heterogeneity of chromatin accessibility in MT CMML, we calculated the single-cell entropies for both genotypes. There was an increase in the dispersion of single-cell entropy values (mean 0.793±0.026 versus 0.811±0.012, p=6.19x10-151) in MT compared to WT disease (Figure 1b). To leverage the power of the single-cell chromatin accessibility data, we estimated cis-co-accessibility scores (two-way chromatin interactions). There were 631 genes involved in two-way chromatin interactions that were either specific to MT CMML or were shared between MT and WT CMML and showed an increase in co-accessibility in MT CMML. Co-accessibility was higher in MT compared to WT CMML, regardless whether the chromatin interaction was MT-specific or shared between the genotypes (Figure 1c). Gene ontology analysis revealed increased chromatin accessibility in genes involved in the regulation of cell cycle, nuclear division, chromosome organization, cell cycle phase transition, DNA conformation change, and DNA replication (FDR<1.00x10-5). The construction of cis-co-accessibility networks demonstrated a remodeling of the chromatin landscape with MT-specific chromatin interactions involving key mitotic kinases such as MAP2K3 (Figure 1d). Conclusions: MT CMML is characterized by an increase in intratumoral heterogeneity and increased chromatin accessibility. The increased chromatin co-accessibility involved key mitotic kinases and may serve as a plausible explanation for the increased transcriptional activity observed in MT CMML. Figure 1 Disclosures Ordog: Millipore Sigma: Patents & Royalties.

Abstract Adenosquamous cancer of the pancreas (ASCP) is a subtype of pancreatic cancer that has a worse prognosis and greater metastatic potential than the more common pancreatic ductal adenocarcinoma (PDAC) subtype. To distinguish the genomic landscape of ASCP and identify actionable targets for this lethal cancer, we applied DNA content flow cytometry to a series of 15 tumor samples including five patient-derived xenografts (PDX). We interrogated purified sorted tumor fractions from these samples with whole-genome copy-number variant (CNV), whole-exome sequencing, and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) analyses. These identified a variety of somatic genomic lesions targeting chromatin regulators in ASCP genomes that were superimposed on well-characterized genomic lesions including mutations in TP53 (87%) and KRAS (73%), amplification of MYC (47%), and homozygous deletion of CDKN2A (40%) that are common in PDACs. Furthermore, a comparison of ATAC-seq profiles of three ASCP and three PDAC genomes using flow-sorted PDX models identified genes with accessible chromatin unique to the ASCP genomes, including the lysine methyltransferase SMYD2 and the pancreatic cancer stem cell regulator RORC in all three ASCPs, and a FGFR1-ERLIN2 fusion associated with focal CNVs in both genes in a single ASCP. Finally, we demonstrate significant activity of a pan FGFR inhibitor against organoids derived from the FGFR1-ERLIN2 fusion–positive ASCP PDX model. Our results suggest that the genomic and epigenomic landscape of ASCP provide new strategies for targeting this aggressive subtype of pancreatic cancer. Significance: These data provide a unique description of the ASCP genomic and epigenomic landscape and identify candidate therapeutic targets for this dismal cancer.

INTRODUCTION: Epigenetic modifications have been implicated to mediate several complications of diabetes mellitus (DM), especially nephropathy and retinopathy. Our aim was to ascertain whether epigenetic alterations in whole blood discriminate among patients with DM with normal, delayed, and rapid gastric emptying (GE). METHODS: Using the ChIP-seq (chromatin immunoprecipitation combined with next-generation sequencing) assays, we compared the genome-wide enrichment of 3 histone modifications (i.e., H3K4me3, H3K9ac, and H3K27ac) in buffy coats from 20 diabetic patients with gastrointestinal symptoms and normal (n = 6), delayed (n = 8), or rapid (n = 6) GE. RESULTS: Between patients with DM with delayed vs normal GE, there were 108 and 54 genes that were differentially bound (false discovery rate < 0.05) with H3K27ac and H3K9ac, respectively; 100 genes were differentially bound with H3K9ac in patients with rapid vs normal GE. The differentially bound genes with H3K27ac were functionally linked to the type 2 immune response, particularly Th2 cell activation and function (e.g., CCR3 , CRLF2 , CXCR4 , IL5RA , and IL1RL1 ) and glucose homeostasis ( FBP-1 , PDE4A , and CMKLR1 ). For H3K9ac, the differentially occupied genes were related to T-cell development and function (e.g., ICOS and CCR3 ) and innate immunity ( RELB , CD300LB , and CLEC2D ). Compared with normal GE, rapid GE had differential H3K9ac peaks at the promoter site of diverse immunity-related genes (e.g., TNFRSF25 and CXCR4 ) and genes related to insulin resistance and glucose metabolism. Motif analysis disclosed enrichment of binding sites for transcription factors relevant to the pathogenesis and complications of DM. DISCUSSION: GE disturbances in DM are associated with epigenetic alterations that pertain to dysimmunity, glucose metabolism, and other complications of DM.

Objective This study was designed to evaluate the roles of microRNAs (miRNAs) in slow transit constipation (STC). Design All human tissue samples were from the muscularis externa of the colon. Expression of 372 miRNAs was examined in a discovery cohort of four patients with STC versus three age/sex-matched controls by a quantitative PCR array. Upregulated miRNAs were examined by quantitative reverse transcription qPCR (RT-qPCR) in a validation cohort of seven patients with STC and age/sex-matched controls. The effect of a highly differentially expressed miRNA on a custom human smooth muscle cell line was examined in vitro by RT-qPCR, electrophysiology, traction force microscopy, and ex vivo by lentiviral transduction in rat muscularis externa organotypic cultures. Results The expression of 13 miRNAs was increased in STC samples. Of those miRNAs, four were predicted to target SCN5A , the gene that encodes the Na + channel Na V 1.5. The expression of SCN5A mRNA was decreased in STC samples. Let-7f significantly decreased Na + current density in vitro in human smooth muscle cells. In rat muscularis externa organotypic cultures, overexpression of let-7f resulted in reduced frequency and amplitude of contraction. Conclusions A small group of miRNAs is upregulated in STC, and many of these miRNAs target the SCN5A-encoded Na + channel Na V 1.5. Within this set, a novel Na V 1.5 regulator, let-7f, resulted in decreased Na V 1.5 expression, current density and reduced motility of GI smooth muscle. These results suggest Na V 1.5 and miRNAs as novel diagnostic and potential therapeutic targets in STC.

Abstract Establishment and subsequent maintenance of distinct chromatin domains during embryonic stem cell (ESC) differentiation are crucial for lineage specification and cell fate determination. Here we show that the histone chaperone Chromatin Assembly Factor 1 (CAF-1), which is recruited to DNA replication forks through its interaction with proliferating cell nuclear antigen (PCNA) for nucleosome assembly, participates in the establishment of H3K27me3-mediated silencing during differentiation. Deletion of CAF-1 p150 subunit impairs the silencing of many genes including Oct4, Sox2 and Nanog as well as the establishment of H3K27me3 at these gene promoters during ESC differentiation. Mutations of PCNA residues involved in recruiting CAF-1 to the chromatin also result in defects in differentiation in vitro and impair early embryonic development as p150 deletion. Together, these results reveal that the CAF-1-PCNA nucleosome assembly pathway plays an important role in the establishment of H3K27me3-mediated silencing during cell fate determination.

ABSTRACT Tissue fibrosis is a chronic disease driven by persistent fibroblast activation that has recently been linked to epigenetic modifications. Here, we screened a small library of epigenetic small-molecule modulators to identify compounds capable of inhibiting or reversing TGFβ-mediated fibroblast activation. We identified pracinostat, an HDAC inhibitor, as a potent attenuator of lung fibroblast activation and confirmed its efficacy in patient-derived fibroblasts isolated from fibrotic lung tissue. Mechanistically, we found that HDAC-dependent transcriptional repression was an early and essential event in TGFβ-mediated fibroblast activation. Treatment of lung fibroblasts with pracinostat broadly attenuated TGFβ-mediated epigenetic repression and promoted fibroblast quiescence. We confirmed a specific role for HDAC-dependent histone deacetylation in the promoter region of the anti-fibrotic gene PPARGC1A (PGC1α) in response to TGFβ stimulation. Finally, we identified HDAC7 as a key factor whose siRNA-mediated knockdown attenuates fibroblast activation without altering global histone acetylation. Together, these results provide novel mechanistic insight into the essential role HDACs play in TGFβ-mediated fibroblast activation via targeted gene repression.

Abstract Background: Epigenetic dysregulation is involved in the etiology and progression of various human diseases. Formalin-fixed paraffin-embedded (FFPE) samples represent the gold standard for archiving pathology samples, and thus FFPE samples are a major resource of samples in clinical research. However, chromatin-based epigenetic assays in the clinical settings are limited to fresh or frozen samples and are hampered by low chromatin yield in FFPE samples due to the lack of reliable and efficient chromatin preparation methods. Here, we introduce a new chromatin extraction method from FFPE tissues (Chrom-EX PE) for chromatin-based epigenetic assays. Results: During rehydration of FFPE tissues, applying tissue-level cross-link reversal to the deparaffinized tissue at 65°C dramatically increased chromatin yield in the soluble fraction. We found the resulting chromatin to be compatible with ChIP-qPCR and ChIP-seq analysis of histone marks (H3K4me3, H3K27me3, and H3K27Ac) and RNA polymerase II. The chromatin prepared by Chrom-EX PE showed a gradual fragmentation pattern with varying incubation temperature. At temperatures below 37°C, the majority of soluble chromatin was over 1 kb in size. The soluble chromatin prepared at 45-60°C showed a typical nucleosomal pattern shown in cell lines and frozen tissues. The majority of the chromatin prepared at 65°C is close to mononucleosomal size. These observations indicate that chromatin preparation from FFPE samples can be controlled for downstream chromatin-based epigenetic assays. Conclusions: Our method achieves efficient extraction of high-quality chromatin from FFPE samples suitable for targeted and genome-wide epigenetic assays while minimizing degradation and damages. This approach provides a new platform for developing chromatin-based assays applicable to archived FFPE samples to circumvent the over-fixed nature of FFPE tissues. This work was supported, in part, by National Institutes of Health grant P01DK068055 and the Mayo Clinic Center for Individualized Medicine Epigenomics Program. Citation Format: Jian Zhong, Zhenqing Ye, Chad Clark, Samuel Lenz, Justin Nguyen, Huihuang Yan, Keith Robertson, Gianrico Farrugia, Zhiguo Zhang, Tamas Ordog, Jeong-Heon Lee. Enhanced and controlled chromatin extraction for chromatin-based epigenetic assays in FFPE tissues [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 5180.

AbstractAnaplastic large cell lymphomas (ALCLs) represent a relatively common group of T-cell non-Hodgkin lymphomas (T-NHLs) that are unified by similar pathologic features but demonstrate marked genetic heterogeneity. ALCLs are broadly classified as being anaplastic lymphoma kinase (ALK)+ or ALK−, based on the presence or absence of ALK rearrangements. Exome sequencing of 62 T-NHLs identified a previously unreported recurrent mutation in the musculin gene, MSCE116K, exclusively in ALK− ALCLs. Additional sequencing for a total of 238 T-NHLs confirmed the specificity of MSCE116K for ALK− ALCL and further demonstrated that 14 of 15 mutated cases (93%) had coexisting DUSP22 rearrangements. Musculin is a basic helix-loop-helix (bHLH) transcription factor that heterodimerizes with other bHLH proteins to regulate lymphocyte development. The E116K mutation localized to the DNA binding domain of musculin and permitted formation of musculin–bHLH heterodimers but prevented their binding to authentic target sequence. Functional analysis showed MSCE116K acted in a dominant-negative fashion, reversing wild-type musculin-induced repression of MYC and cell cycle inhibition. Chromatin immunoprecipitation–sequencing and transcriptome analysis identified the cell cycle regulatory gene E2F2 as a direct transcriptional target of musculin. MSCE116K reversed E2F2-induced cell cycle arrest and promoted expression of the CD30–IRF4–MYC axis, whereas its expression was reciprocally induced by binding of IRF4 to the MSC promoter. Finally, ALCL cells expressing MSCE116K were preferentially targeted by the BET inhibitor JQ1. These findings identify a novel recurrent MSC mutation as a key driver of the CD30–IRF4–MYC axis and cell cycle progression in a unique subset of ALCLs.