EpiFactors 2022: expansion and enhancement of a curated database of human epigenetic factors and complexes.

We present an update of EpiFactors, a manually curated database providing information about epigenetic regulators, their complexes, targets, and products which is openly accessible at http://epifactors.autosome.org. An updated version of the EpiFactors contains information on 902 proteins, including 101 histones and protamines, and, as a main update, a newly curated collection of 124 lncRNAs involved in epigenetic regulation. The amount of publications concerning the role of lncRNA in epigenetics is rapidly growing. Yet, the resource that compiles, integrates, organizes, and presents curated information on lncRNAs in epigenetics is missing. EpiFactors fills this gap and provides data on epigenetic regulators in an accessible and user-friendly form. For 820 of the genes in EpiFactors, we include expression estimates across multiple cell types assessed by CAGE-Seq in the FANTOM5 project. In addition, the updated EpiFactors contains information on 73 protein complexes involved in epigenetic regulation. Our resource is practical for a wide range of users, including biologists, bioinformaticians and molecular/systems biologists.

AML displays increased CTCF occupancy associated with aberrant gene expression and transcription factor binding.

CCTC-binding factor (CTCF) is a key regulator of gene expression through organization of the chromatin structure. Still, it is unclear how CTCF binding is perturbed in leukemia or in cancer in general. We studied CTCF binding by chromatin immunoprecipitation sequencing in cells from patients with acute myeloid leukemia (AML) and in normal bone marrow (NBM) in the context of gene expression, DNA methylation, and azacitidine exposure. CTCF binding was increased in AML compared with NBM. Aberrant CTCF binding was enriched for motifs for key myeloid transcription factors such as CEBPA, PU.1, and RUNX1. AML with TET2 mutations was characterized by a particularly strong gain of CTCF binding, highly enriched for gain in promoter regions, while AML in general was enriched for changes at enhancers. There was a strong anticorrelation between CTCF binding and DNA methylation. Gain of CTCF occupancy was associated with increased gene expression; however, the genomic location (promoter vs distal regions) and enrichment of motifs (for repressing vs activating cofactors) were decisive for the gene expression pattern. Knockdown of CTCF in K562 cells caused loss of CTCF binding and transcriptional repression of genes with changed CTCF binding in AML, as well as loss of RUNX1 binding at RUNX1/CTCF-binding sites. In addition, CTCF knockdown caused increased differentiation. Azacitidine exposure caused major changes in CTCF occupancy in AML patient cells, partly by restoring a CTCF-binding pattern similar to NBM. We conclude that AML displays an aberrant increase in CTCF occupancy that targets key genes for AML development and impacts gene expression.

RUNX1/CEBPA Mutation in Acute Myeloid Leukemia Promotes Hypermethylation and Indicates for Demethylation Therapy.

Acute myeloid leukemia (AML) is a rapidly progressing heterogeneous disease with a high mortality rate, which is characterized by hyperproliferation of atypical immature myeloid cells. The number of AML patients is expected to increase in the near future, due to the old-age-associated nature of AML and increased longevity in the human population. RUNX1 and CEBPA, key transcription factors (TFs) of hematopoiesis, are frequently and independently mutated in AML. RUNX1 and CEBPA can bind TET2 demethylase and attract it to their binding sites (TFBS) in cell lines, leading to DNA demethylation of the regions nearby. Since TET2 does not have a DNA-binding domain, TFs are crucial for its guidance to target genomic locations. In this paper, we show that RUNX1 and CEBPA mutations in AML patients affect the methylation of important regulatory sites that resulted in the silencing of several RUNX1 and CEBPA target genes, most likely in a TET2-dependent manner. We demonstrated that hypermethylation of TFBS in AML cells with RUNX1 mutations was associated with resistance to anticancer chemotherapy. Demethylation therapy restored expression of the RUNX1 target gene, BIK, and increased sensitivity of AML cells to chemotherapy. If our results are confirmed, mutations in RUNX1 could be an indication for prescribing the combination of cytotoxic and demethylation therapies.

Sirtuins and neuropeptide y downregulation in Flinders Sensitive Line rat model of depression.

Since the NAD+-dependent histone deacetylases sirtuin-1 (SIRT1) and sirtuin-2 (SIRT2) are critically involved in epigenetics, endocrinology and immunology and affect the longevity in model organisms, we investigated their expression in brains of 3-month-old and 14-15 months old rat model of depression Flinders Sensitive Line (FSL) and control Flinders Resistant Line (FRL) rats. In view of the dysregulated NPY system in depression, we also studied NPY in young and old FSL to explore the temporal trajectory of depressive-like-ageing interaction. Sirt1, Sirt2 and Npy mRNA were determined using qRT-PCR in prefrontal cortex (PFC) from young and old FSL and FRL, and in hippocampi from young FSL and FRL. PFC: Sirt1 expression was decreased in FSL (p = 0.001). An interaction between age and genotype was found (p = 0.032); young FSL had lower Sirt1 with respect to both age (p = 0.026) and genotype (p = 0.001). Sirt2 was lower in FSL (p = 0.003). Npy mRNA was downregulated in FSL (p = 0.001) but did not differ between the young and old rat groups. Hippocampus: Sirt1 was reduced in young FSL compared to young FRL (p = 0.005). There was no difference in Sirt2 between FSL and FRL. Npy levels were decreased in hippocampus of young FSL compared to young FRL (p = 0.003). Effects of ageing could not be investigated due to loss of samples. To conclude, i this is the first demonstration that SIRT1 and SIRT2 are changed in brain of FSL, a rat model of depression; ii the changes are age-dependent; iii sirtuins are potential targets for treatment of age-related neurodegenerative diseases.

The Role of Non-Catalytic Domains of Hrp3 in Nucleosome Remodeling.

The Helicase-related protein 3 (Hrp3), an ATP-dependent chromatin remodeling enzyme from the CHD family, is crucial for maintaining global nucleosome occupancy in Schizosaccharomyces pombe (S. pombe). Although the ATPase domain of Hrp3 is essential for chromatin remodeling, the contribution of non-ATPase domains of Hrp3 is still unclear. Here, we investigated the role of non-ATPase domains using in vitro methods. In our study, we expressed and purified recombinant S. pombe histone proteins, reconstituted them into histone octamers, and assembled nucleosome core particles. Using reconstituted nucleosomes and affinity-purified wild type and mutant Hrp3 from S. pombe we created a homogeneous in vitro system to evaluate the ATP hydrolyzing capacity of truncated Hrp3 proteins. We found that all non-ATPase domain deletions (∆chromo, ∆SANT, ∆SLIDE, and ∆coupling region) lead to reduced ATP hydrolyzing activities in vitro with DNA or nucleosome substrates. Only the coupling region deletion showed moderate stimulation of ATPase activity with the nucleosome. Interestingly, affinity-purified Hrp3 showed co-purification with all core histones suggesting a strong association with the nucleosomes in vivo. However, affinity-purified Hrp3 mutant with SANT and coupling regions deletion showed complete loss of interactions with the nucleosomes, while SLIDE and chromodomain deletions reduced Hrp3 interactions with the nucleosomes. Taken together, nucleosome association and ATPase stimulation by DNA or nucleosomes substrate suggest that the enzymatic activity of Hrp3 is fine-tuned by unique contributions of all four non-catalytic domains.

An atlas of active enhancers across human cell types and tissues.

Enhancers control the correct temporal and cell-type-specific activation of gene expression in multicellular eukaryotes. Knowing their properties, regulatory activity and targets is crucial to understand the regulation of differentiation and homeostasis. Here we use the FANTOM5 panel of samples, covering the majority of human tissues and cell types, to produce an atlas of active, in vivo-transcribed enhancers. We show that enhancers share properties with CpG-poor messenger RNA promoters but produce bidirectional, exosome-sensitive, relatively short unspliced RNAs, the generation of which is strongly related to enhancer activity. The atlas is used to compare regulatory programs between different cells at unprecedented depth, to identify disease-associated regulatory single nucleotide polymorphisms, and to classify cell-type-specific and ubiquitous enhancers. We further explore the utility of enhancer redundancy, which explains gene expression strength rather than expression patterns. The online FANTOM5 enhancer atlas represents a unique resource for studies on cell-type-specific enhancers and gene regulation.

Cancer-specific changes in DNA methylation reveal aberrant silencing and activation of enhancers in leukemia.

Acute myeloid leukemia (AML) is characterized by an impaired differentiation process leading to an accumulation of immature blasts in the blood. One feature of cytogenetically normal AML is alterations to the DNA methylome. We analyzed 57 AML patients with normal karyotype by using Illumina's 450k array and showed that aberrant DNA methylation is significantly altered at enhancer regions and that the methylation levels at specific enhancers predict overall survival of AML patients. The majority of sites that become differentially methylated in AML occur in regulatory elements of the human genome. Hypermethylation associates with enhancer silencing. In addition, chromatin immunoprecipitation sequencing analyses showed that a subset of hypomethylated sites correlate with enhancer activation, indicated by increased H3K27 acetylation. DNA hypomethylation is therefore not sufficient for enhancer activation. Some sites of hypomethylation occur at weak/poised enhancers marked with H3K4 monomethylation in hematopoietic progenitor cells. Other hypomethylated regions occur at sites inactive in progenitors and reflect the de novo acquisition of AML-specific enhancers. Altered enhancer dynamics are reflected in the gene expression of enhancer target genes, including genes involved in oncogenesis and blood cell development. This study demonstrates that histone variants and different histone modifications interact with aberrant DNA methylation and cause perturbed enhancer activity in cytogenetically normal AML that contributes to a leukemic transcriptome.

Update of the FANTOM web resource: high resolution transcriptome of diverse cell types in mammals.

Upon the first publication of the fifth iteration of the Functional Annotation of Mammalian Genomes collaborative project, FANTOM5, we gathered a series of primary data and database systems into the FANTOM web resource (http://fantom.gsc.riken.jp) to facilitate researchers to explore transcriptional regulation and cellular states. In the course of the collaboration, primary data and analysis results have been expanded, and functionalities of the database systems enhanced. We believe that our data and web systems are invaluable resources, and we think the scientific community will benefit for this recent update to deepen their understanding of mammalian cellular organization. We introduce the contents of FANTOM5 here, report recent updates in the web resource and provide future perspectives.

Distinct global binding patterns of the Wilms tumor gene 1 (WT1) -KTS and +KTS isoforms in leukemic cells.

The zinc finger transcription factor Wilms tumor gene 1 (WT1) acts as an oncogene in acute myeloid leukemia. A naturally occurring alternative splice event between zinc fingers three and four, removing or retaining three amino acids (±KTS), is believed to change the DNA binding affinity of WT1, although there are conflicting data regarding the binding affinity and motifs of the different isoforms. Increased expression of the WT1 -KTS isoform at the expense of the WT1 +KTS isoform is associated with poor prognosis in acute myeloid leukemia. We determined the genome-wide binding pattern of WT1 -KTS and WT1 +KTS in leukemic K562 cells by chromatin immunoprecipitation and deep sequencing. We discovered that the WT1 -KTS isoform predominantly binds close to transcription start sites and to enhancers, in a similar fashion to other transcription factors, whereas WT1 +KTS binding is enriched within gene bodies. We observed a significant overlap between WT1 -KTS and WT1 +KTS target genes, despite the binding sites being distinct. Motif discovery revealed distinct binding motifs for the isoforms, some of which have been previously reported as WT1 binding sites. Additional analyses showed that both WT1 -KTS and WT1 +KTS target genes are more likely to be transcribed than non-targets, and are involved in cell proliferation, cell death, and development. Our study provides evidence that WT1 -KTS and WT1 +KTS share target genes yet still bind distinct locations, indicating isoform-specific regulation in transcription of genes related to cell proliferation and differentiation, consistent with the involvement of WT1 in acute myeloid leukemia.

High-throughput transcription profiling identifies putative epigenetic regulators of hematopoiesis.

Hematopoietic differentiation is governed by a complex regulatory program controlling the generation of different lineages of blood cells from multipotent hematopoietic stem cells. The transcriptional program that dictates hematopoietic cell fate and differentiation requires an epigenetic memory function provided by a network of epigenetic factors regulating DNA methylation, posttranslational histone modifications, and chromatin structure. Aberrant interactions between epigenetic factors and transcription factors cause perturbations in the blood cell differentiation program that result in various types of hematopoietic disorders. To elucidate the contributions of different epigenetic factors in human hematopoiesis, high-throughput cap analysis of gene expression was used to build transcription profiles of 199 epigenetic factors in a wide range of blood cells. Our epigenetic transcriptome analysis revealed cell type- (eg, HELLS and ACTL6A), lineage- (eg, MLL), and/or leukemia- (eg, CHD2, CBX8, and EPC1) specific expression of several epigenetic factors. In addition, we show that several epigenetic factors use alternative transcription start sites in different cell types. This analysis could serve as a resource for the scientific community for further characterization of the role of these epigenetic factors in blood development.

Analysis of the DNA methylome and transcriptome in granulopoiesis reveals timed changes and dynamic enhancer methylation.

In development, epigenetic mechanisms such as DNA methylation have been suggested to provide a cellular memory to maintain multipotency but also stabilize cell fate decisions and direct lineage restriction. In this study, we set out to characterize changes in DNA methylation and gene expression during granulopoiesis using 4 distinct cell populations ranging from the oligopotent common myeloid progenitor stage to terminally differentiated neutrophils. We observed that differentially methylated sites (DMSs) generally show decreased methylation during granulopoiesis. Methylation appears to change at specific differentiation stages and overlap with changes in transcription and activity of key hematopoietic transcription factors. DMSs were preferentially located in areas distal to CpG islands and shores. Also, DMSs were overrepresented in enhancer elements and enriched in enhancers that become active during differentiation. Overall, this study depicts in detail the epigenetic and transcriptional changes that occur during granulopoiesis and supports the role of DNA methylation as a regulatory mechanism in blood cell differentiation.

The DEK oncoprotein binds to highly and ubiquitously expressed genes with a dual role in their transcriptional regulation.

BACKGROUND: The DEK gene is highly expressed in a wide range of cancer cells, and a recurrent translocation partner in acute myeloid leukemia. While DEK has been identified as one of the most abundant proteins in human chromatin, its function and binding properties are not fully understood. METHODS: We performed ChIP-seq analysis in the myeloid cell line U937 and coupled it with epigenetic and gene expression analysis to explore the genome-wide binding pattern of DEK and its role in gene regulation. RESULTS: We show that DEK preferentially binds to open chromatin, with a low degree of DNA methylation and scarce in the heterochromatin marker H3K9me(3) but rich in the euchromatin marks H3K4me(2/3), H3K27ac and H3K9ac. More specifically, DEK binding is predominantly located at the transcription start sites of highly transcribed genes and a comparative analysis with previously established transcription factor binding patterns shows a similarity with that of RNA polymerase II. Further bioinformatic analysis demonstrates that DEK mainly binds to genes that are ubiquitously expressed across tissues. The functional significance of DEK binding was demonstrated by knockdown of DEK by shRNA, resulting in both significant upregulation and downregulation of DEK-bound genes. CONCLUSIONS: We find that DEK binds to transcription start sites with a dual role in activation and repression of highly and ubiquitously expressed genes.

Prognostic DNA methylation patterns in cytogenetically normal acute myeloid leukemia are predefined by stem cell chromatin marks.

Cytogenetically normal acute myeloid leukemia (CN-AML) compose between 40% and 50% of all adult acute myeloid leukemia (AML) cases. In this clinically diverse group, molecular aberrations, such as FLT3-ITD, NPM1, and CEBPA mutations, recently have added to the prognostic accuracy. Aberrant DNA methylation is a hallmark of cancer, including AML. We investigated in total 118 CN-AML samples in a test and a validation cohort for genome-wide promoter DNA methylation with Illumina Methylation Bead arrays and compared them with normal myeloid precursors and global gene expression. IDH and NPM1 mutations were associated with different methylation patterns (P = .0004 and .04, respectively). Genome-wide methylation levels were elevated in IDH-mutated samples (P = .006). We observed a negative impact of DNA methylation on transcription. Genes targeted by Polycomb group (PcG) proteins and genes associated with bivalent histone marks in stem cells showed increased aberrant methylation in AML (P < .0001). Furthermore, high methylation levels of PcG target genes were independently associated with better progression-free survival (odds ratio = 0.47, P = .01) and overall survival (odds ratio = 0.36, P = .001). In summary, genome-wide methylation patterns show preferential methylation of PcG targets with prognostic impact in CN-AML.

Forced expression of the DEK-NUP214 fusion protein promotes proliferation dependent on upregulation of mTOR.

BACKGROUND: The t(6;9)(p23;q34) chromosomal translocation is found in 1% of acute myeloid leukemia and encodes the fusion protein DEK-NUP214 (formerly DEK-CAN) with largely uncharacterized functions. METHODS: We expressed DEK-NUP214 in the myeloid cell lines U937 and PL-21 and studied the effects on cellular functions. RESULTS: In this study, we demonstrate that expression of DEK-NUP214 increases cellular proliferation. Western blot analysis revealed elevated levels of one of the key proteins regulating proliferation, the mechanistic target of rapamycin, mTOR. This conferred increased mTORC1 but not mTORC2 activity, as determined by the phosphorylation of their substrates, p70 S6 kinase and Akt. The functional importance of the mTOR upregulation was determined by assaying the downstream cellular processes; protein synthesis and glucose metabolism. A global translation assay revealed a substantial increase in the translation rate and a metabolic assay detected a shift from glycolysis to oxidative phosphorylation, as determined by a reduction in lactate production without a concomitant decrease in glucose consumption. Both these effects are in concordance with increased mTORC1 activity. Treatment with the mTORC1 inhibitor everolimus (RAD001) selectively reversed the DEK-NUP214-induced proliferation, demonstrating that the effect is mTOR-dependent. CONCLUSIONS: Our study shows that the DEK-NUP214 fusion gene increases proliferation by upregulation of mTOR, suggesting that patients with leukemias carrying DEK-NUP214 may benefit from treatment with mTOR inhibitors.

Perturbed epigenetic transcriptional regulation in AML with IDH mutations causes increased susceptibility to NK cells.

Isocitrate dehydrogenase (IDH) mutations are found in 20% of acute myeloid leukemia (AML) patients. However, only 30-40% of the patients respond to IDH inhibitors (IDHi). We aimed to identify a molecular vulnerability to tailor novel therapies for AML patients with IDH mutations. We characterized the transcriptional and epigenetic landscape with the IDH2i AG-221, using an IDH2 mutated AML cell line model and AML patient cohorts, and discovered a perturbed transcriptional regulatory network involving myeloid transcription factors that were partly restored after AG-221 treatment. In addition, hypermethylation of the HLA cluster caused a down-regulation of HLA class I genes, triggering an enhanced natural killer (NK) cell activation and an increased susceptibility to NK cell-mediated responses. Finally, analyses of DNA methylation data from IDHi-treated patients showed that non-responders still harbored hypermethylation in HLA class I genes. In conclusion, this study provides new insights suggesting that IDH mutated AML is particularly sensitive to NK cell-based personalized immunotherapy.

Antidepressant treatment is associated with epigenetic alterations in the promoter of P11 in a genetic model of depression.

P11 (S100A10) has been associated with the pathophysiology of depression both in human and rodent models. Different types of antidepressants have been shown to increase P11 levels in distinct brain regions and P11 gene therapy was recently proven effective in reversing depressive-like behaviours in mice. However, the molecular mechanisms that govern P11 gene expression in response to antidepressants still remain elusive. In this study we report decreased levels of P11, associated with higher DNA methylation in the promoter region, in the prefrontal cortex of the Flinders Sensitive Line (FSL) genetic rodent model of depression. This hypermethylated pattern was reversed to normal, as indicated by the control line, after chronic administration of escitalopram (a selective serotonin reuptake inhibitor; SSRI). The escitalopram-induced hypomethylation was associated with both an increase in P11 gene expression and a reduction in mRNA levels of two DNA methyltransferases that have been shown to maintain DNA methylation in adult forebrain neurons (Dnmt1 and Dnmt3a). In conclusion, our data further support a role for P11 in depression-like states and suggest that this gene is controlled by epigenetic mechanisms that can be affected by antidepressant treatment.

One Omics Approach Does Not Rule Them All: The Metabolome and the Epigenome Join Forces in Haematological Malignancies.

Aberrant DNA methylation, dysregulation of chromatin-modifying enzymes, and microRNAs (miRNAs) play a crucial role in haematological malignancies. These epimutations, with an impact on chromatin accessibility and transcriptional output, are often associated with genomic instability and the emergence of drug resistance, disease progression, and poor survival. In order to exert their functions, epigenetic enzymes utilize cellular metabolites as co-factors and are highly dependent on their availability. By affecting the expression of metabolic enzymes, epigenetic modifiers may aid the generation of metabolite signatures that could be utilized as targets and biomarkers in cancer. This interdependency remains often neglected and poorly represented in studies, despite well-established methods to study the cellular metabolome. This review critically summarizes the current knowledge in the field to provide an integral picture of the interplay between epigenomic alterations and the cellular metabolome in haematological malignancies. Our recent findings defining a distinct metabolic signature upon response to enhancer of zeste homolog 2 (EZH2) inhibition in multiple myeloma (MM) highlight how a shift of preferred metabolic pathways may potentiate novel treatments. The suggested link between the epigenome and the metabolome in haematopoietic tumours holds promise for the use of metabolic signatures as possible biomarkers of response to treatment.

Challenging conventional karyotyping by next-generation karyotyping in 281 intensively treated patients with AML.

Although copy number alterations (CNAs) and translocations constitute the backbone of the diagnosis and prognostication of acute myeloid leukemia (AML), techniques used for their assessment in routine diagnostics have not been reconsidered for decades. We used a combination of 2 next-generation sequencing-based techniques to challenge the currently recommended conventional cytogenetic analysis (CCA), comparing the approaches in a series of 281 intensively treated patients with AML. Shallow whole-genome sequencing (sWGS) outperformed CCA in detecting European Leukemia Net (ELN)-defining CNAs and showed that CCA overestimated monosomies and suboptimally reported karyotype complexity. Still, the concordance between CCA and sWGS for all ELN CNA-related criteria was 94%. Moreover, using in silico dilution, we showed that 1 million reads per patient would be enough to accurately assess ELN-defining CNAs. Total genomic loss, defined as a total loss ≥200 Mb by sWGS, was found to be a better marker for genetic complexity and poor prognosis compared with the CCA-based definition of complex karyotype. For fusion detection, the concordance between CCA and whole-transcriptome sequencing (WTS) was 99%. WTS had better sensitivity in identifying inv(16) and KMT2A rearrangements while showing limitations in detecting lowly expressed PML-RARA fusions. Ligation-dependent reverse transcription polymerase chain reaction was used for validation and was shown to be a fast and reliable method for fusion detection. We conclude that a next-generation sequencing-based approach can replace conventional CCA for karyotyping, provided that efforts are made to cover lowly expressed fusion transcripts.

Histone modification patterns and epigenetic codes.

The eukaryotic DNA is wrapped around histone octamers, which consist of four different histones, H2A, H2B, H3 and H4. The N-terminal tail of each histone is post-transcriptionally modified. The modification patterns constitute codes that regulate chromatin organisation and DNA utilization processes, including transcription. Recent progress in technology development has made it possible to perform systematic genome-wide studies of histone modifications. This helps immensely in deciphering the histone codes and their biological influence. In this review, we discuss the histone modification patterns found in genome-wide studies in different biological models and how they influence cell differentiation and carcinogenesis.

Abo1 is required for the H3K9me2 to H3K9me3 transition in heterochromatin.

Heterochromatin regulation is critical for genomic stability. Different H3K9 methylation states have been discovered, with distinct roles in heterochromatin formation and silencing. However, how the transition from H3K9me2 to H3K9me3 is controlled is still unclear. Here, we investigate the role of the conserved bromodomain AAA-ATPase, Abo1, involved in maintaining global nucleosome organisation in fission yeast. We identified several key factors involved in heterochromatin silencing that interact genetically with Abo1: histone deacetylase Clr3, H3K9 methyltransferase Clr4, and HP1 homolog Swi6. Cells lacking Abo1 cultivated at 30 °C exhibit an imbalance of H3K9me2 and H3K9me3 in heterochromatin. In abo1∆ cells, the centromeric constitutive heterochromatin has increased H3K9me2 but decreased H3K9me3 levels compared to wild-type. In contrast, facultative heterochromatin regions exhibit reduced H3K9me2 and H3K9me3 levels in abo1∆. Genome-wide analysis showed that abo1∆ cells have silencing defects in both the centromeres and subtelomeres, but not in a subset of heterochromatin islands in our condition. Thus, our work uncovers a role of Abo1 in stabilising directly or indirectly Clr4 recruitment to allow the H3K9me2 to H3K9me3 transition in heterochromatin.