Figeno: multi-region genomic figures with long-read support.

SUMMARY: The vast amount of publicly available genomic data requires analysis and visualization tools. Here, we present figeno, an application for generating publication-quality FIgures for GENOmics. Figeno particularly focuses on multi-region views across genomic breakpoints and on long reads with base modifications. In addition, we support epigenomic data including ATAC-seq, ChIP-seq or HiC, as well as whole genome sequencing data with copy numbers and structural variants. AVAILABILITY AND IMPLEMENTATION: Figeno is available as a python package with both a command line and graphical user interface. It can be installed via PyPI and the source code is available at https://github.com/CompEpigen/figeno.

ARID1A regulates DNA repair through chromatin organization and its deficiency triggers DNA damage-mediated anti-tumor immune response.

AT-rich interaction domain protein 1A (ARID1A), a SWI/SNF chromatin remodeling complex subunit, is frequently mutated across various cancer entities. Loss of ARID1A leads to DNA repair defects. Here, we show that ARID1A plays epigenetic roles to promote both DNA double-strand breaks (DSBs) repair pathways, non-homologous end-joining (NHEJ) and homologous recombination (HR). ARID1A is accumulated at DSBs after DNA damage and regulates chromatin loops formation by recruiting RAD21 and CTCF to DSBs. Simultaneously, ARID1A facilitates transcription silencing at DSBs in transcriptionally active chromatin by recruiting HDAC1 and RSF1 to control the distribution of activating histone marks, chromatin accessibility, and eviction of RNAPII. ARID1A depletion resulted in enhanced accumulation of micronuclei, activation of cGAS-STING pathway, and an increased expression of immunomodulatory cytokines upon ionizing radiation. Furthermore, low ARID1A expression in cancer patients receiving radiotherapy was associated with higher infiltration of several immune cells. The high mutation rate of ARID1A in various cancer types highlights its clinical relevance as a promising biomarker that correlates with the level of immune regulatory cytokines and estimates the levels of tumor-infiltrating immune cells, which can predict the response to the combination of radio- and immunotherapy.

Dnmt3a1 regulates hippocampus-dependent memory via the downstream target Nrp1.

Epigenetic factors are well-established players in memory formation. Specifically, DNA methylation is necessary for the formation of long-term memory in multiple brain regions including the hippocampus. Despite the demonstrated role of DNA methyltransferases (Dnmts) in memory formation, it is unclear whether individual Dnmts have unique or redundant functions in long-term memory formation. Furthermore, the downstream processes controlled by Dnmts during memory consolidation have not been investigated. In this study, we demonstrated that Dnmt3a1, the predominant Dnmt in the adult brain, is required for long-term spatial object recognition and contextual fear memory. Using RNA sequencing, we identified an activity-regulated Dnmt3a1-dependent genomic program in which several genes were associated with functional and structural plasticity. Furthermore, we found that some of the identified genes are selectively dependent on Dnmt3a1, but not its isoform Dnmt3a2. Specifically, we identified Neuropilin 1 (Nrp1) as a downstream target of Dnmt3a1 and further demonstrated the involvement of Nrp1 in hippocampus-dependent memory formation. Importantly, we found that Dnmt3a1 regulates hippocampus-dependent memory via Nrp1. In contrast, Nrp1 overexpression did not rescue memory impairments triggered by reduced Dnmt3a2 levels. Taken together, our study uncovered a Dnmt3a-isoform-specific mechanism in memory formation, identified a novel regulator of memory, and further highlighted the complex and highly regulated functions of distinct epigenetic regulators in brain function.

Deposition of onco-histone H3.3-G34W leads to DNA repair deficiency and activates cGAS/STING-mediated immune responses.

Mutations in histone H3.3-encoding genes causing mutant histone tails are associated with specific cancers such as pediatric glioblastomas (H3.3-G34R/V) and giant cell tumor of the bone (H3.3-G34W). The mechanisms by which these mutations promote malignancy are not completely understood. Here we show that cells expressing H3.3-G34W exhibit DNA double-strand breaks (DSBs) repair defects and increased cellular sensitivity to ionizing radiation (IR). Mechanistically, H3.3-G34W can be deposited to damaged chromatin, but in contrast to wild-type H3.3, does not interact with non-homologous end-joining (NHEJ) key effectors KU70/80 and XRCC4 leading to NHEJ deficiency. Together with defective cell cycle checkpoints reported previously, this DNA repair deficiency in H3.3-G34W cells led to accumulation of micronuclei and cytosolic DNA following IR, which subsequently led to activation of the cyclic GMP-AMP synthase/stimulator of interferon genes (cGAS/STING) pathway, thereby inducing release of immune-stimulatory cytokines. These findings suggest a potential for radiotherapy for tumors expressing H3.3-G34W, which can be further improved by combination with STING agonists to induce immune-mediated therapeutic efficacy.

Multi-omic and functional analysis for classification and treatment of sarcomas with FUS-TFCP2 or EWSR1-TFCP2 fusions.

Linking clinical multi-omics with mechanistic studies may improve the understanding of rare cancers. We leverage two precision oncology programs to investigate rhabdomyosarcoma with FUS/EWSR1-TFCP2 fusions, an orphan malignancy without effective therapies. All tumors exhibit outlier ALK expression, partly accompanied by intragenic deletions and aberrant splicing resulting in ALK variants that are oncogenic and sensitive to ALK inhibitors. Additionally, recurrent CKDN2A/MTAP co-deletions provide a rationale for PRMT5-targeted therapies. Functional studies show that FUS-TFCP2 blocks myogenic differentiation, induces transcription of ALK and truncated TERT, and inhibits DNA repair. Unlike other fusion-driven sarcomas, TFCP2-rearranged tumors exhibit genomic instability and signs of defective homologous recombination. DNA methylation profiling demonstrates a close relationship with undifferentiated sarcomas. In two patients, sarcoma was preceded by benign lesions carrying FUS-TFCP2, indicating stepwise sarcomagenesis. This study illustrates the potential of linking precision oncology with preclinical research to gain insight into the classification, pathogenesis, and therapeutic vulnerabilities of rare cancers.

Aberrant MNX1 expression associated with t(7;12)(q36;p13) pediatric acute myeloid leukemia induces the disease through altering histone methylation.

Certain subtypes of acute myeloid leukemia (AML) in children have inferior outcome, such as AML with translocation t(7;12)(q36;p13) leading to an MNX1::ETV6 fusion along with high expression of MNX1. We have identified the transforming event in this AML and possible ways of treatment. Retroviral expression of MNX1 was able to induce AML in mice, with similar gene expression and pathway enrichment to t(7;12) AML patient data. Importantly, this leukemia was only induced in immune incompetent mice using fetal but not adult hematopoietic stem and progenitor cells. The restriction in transforming capacity to cells from fetal liver is in alignment with t(7;12)(q36;p13) AML being mostly seen in infants. Expression of MNX1 led to increased histone 3 lysine 4 mono-, di- and trimethylation, reduction in H3K27me3, accompanied with changes in genome-wide chromatin accessibility and genome expression, likely mediated through MNX1 interaction with the methionine cycle and methyltransferases. MNX1 expression increased DNA damage, depletion of the Lin-/Sca1+/c-Kit+ population and skewing toward the myeloid lineage. These effects, together with leukemia development, were prevented by pre-treatment with the S-adenosylmethionine analog Sinefungin. In conclusion, we have shown the importance of MNX1 in development of AML with t(7;12), supporting a rationale for targeting MNX1 and downstream pathways.

Predictive value of DNA methylation patterns in AML patients treated with an azacytidine containing induction regimen.

BACKGROUND: Acute myeloid leukemia (AML) is a heterogeneous disease with a poor prognosis. Dysregulation of the epigenetic machinery is a significant contributor to disease development. Some AML patients benefit from treatment with hypomethylating agents (HMAs), but no predictive biomarkers for therapy response exist. Here, we investigated whether unbiased genome-wide assessment of pre-treatment DNA-methylation profiles in AML bone marrow blasts can help to identify patients who will achieve a remission after an azacytidine-containing induction regimen. RESULTS: A total of n = 155 patients with newly diagnosed AML treated in the AMLSG 12-09 trial were randomly assigned to a screening and a refinement and validation cohort. The cohorts were divided according to azacytidine-containing induction regimens and response status. Methylation status was assessed for 664,227 500-bp-regions using methyl-CpG immunoprecipitation-seq, resulting in 1755 differentially methylated regions (DMRs). Top regions were distilled and included genes such as WNT10A and GATA3. 80% of regions identified as a hit were represented on HumanMethlyation 450k Bead Chips. Quantitative methylation analysis confirmed 90% of these regions (36 of 40 DMRs). A classifier was trained using penalized logistic regression and fivefold cross validation containing 17 CpGs. Validation based on mass spectra generated by MALDI-TOF failed (AUC 0.59). However, discriminative ability was maintained by adding neighboring CpGs. A recomposed classifier with 12 CpGs resulted in an AUC of 0.77. When evaluated in the non-azacytidine containing group, the AUC was 0.76. CONCLUSIONS: Our analysis evaluated the value of a whole genome methyl-CpG screening assay for the identification of informative methylation changes. We also compared the informative content and discriminatory power of regions and single CpGs for predicting response to therapy. The relevance of the identified DMRs is supported by their association with key regulatory processes of oncogenic transformation and support the idea of relevant DMRs being enriched at distinct loci rather than evenly distribution across the genome. Predictive response to therapy could be established but lacked specificity for treatment with azacytidine. Our results suggest that a predictive epigenotype carries its methylation information at a complex, genome-wide level, that is confined to regions, rather than to single CpGs. With increasing application of combinatorial regimens, response prediction may become even more complicated.

DNMT and HDAC inhibition induces immunogenic neoantigens from human endogenous retroviral element-derived transcripts.

Immunotherapies targeting cancer-specific neoantigens have revolutionized the treatment of cancer patients. Recent evidence suggests that epigenetic therapies synergize with immunotherapies, mediated by the de-repression of endogenous retroviral element (ERV)-encoded promoters, and the initiation of transcription. Here, we use deep RNA sequencing from cancer cell lines treated with DNA methyltransferase inhibitor (DNMTi) and/or Histone deacetylase inhibitor (HDACi), to assemble a de novo transcriptome and identify several thousand ERV-derived, treatment-induced novel polyadenylated transcripts (TINPATs). Using immunopeptidomics, we demonstrate the human leukocyte antigen (HLA) presentation of 45 spectra-validated treatment-induced neopeptides (t-neopeptides) arising from TINPATs. We illustrate the potential of the identified t-neopeptides to elicit a T-cell response to effectively target cancer cells. We further verify the presence of t-neopeptides in AML patient samples after in vivo treatment with the DNMT inhibitor Decitabine. Our findings highlight the potential of ERV-derived neoantigens in epigenetic and immune therapies.

Clinical impact of the genomic landscape and leukemogenic trajectories in non-intensively treated elderly acute myeloid leukemia patients.

To characterize the genomic landscape and leukemogenic pathways of older, newly diagnosed, non-intensively treated patients with AML and to study the clinical implications, comprehensive genetics analyses were performed including targeted DNA sequencing of 263 genes in 604 patients treated in a prospective Phase III clinical trial. Leukemic trajectories were delineated using oncogenetic tree modeling and hierarchical clustering, and prognostic groups were derived from multivariable Cox regression models. Clonal hematopoiesis-related genes (ASXL1, TET2, SRSF2, DNMT3A) were most frequently mutated. The oncogenetic modeling algorithm produced a tree with five branches with ASXL1, DDX41, DNMT3A, TET2, and TP53 emanating from the root suggesting leukemia-initiating events which gave rise to further subbranches with distinct subclones. Unsupervised clustering mirrored the genetic groups identified by the tree model. Multivariable analysis identified FLT3 internal tandem duplications (ITD), SRSF2, and TP53 mutations as poor prognostic factors, while DDX41 mutations exerted an exceptionally favorable effect. Subsequent backwards elimination based on the Akaike information criterion delineated three genetic risk groups: DDX41 mutations (favorable-risk), DDX41wildtype/FLT3-ITDneg/TP53wildtype (intermediate-risk), and FLT3-ITD or TP53 mutations (high-risk). Our data identified distinct trajectories of leukemia development in older AML patients and provide a basis for a clinically meaningful genetic outcome stratification for patients receiving less intensive therapies.

Embryonic stem cell factor FOXD3 (Genesis) defects in gastrointestinal stromal tumors.

Gastrointestinal stromal tumors (GISTs) are mesenchymal neoplasms, believed to originate from the interstitial cells of Cajal (ICC), often caused by overexpression of tyrosine kinase receptors (TKR) KIT or PDGFRA. Here, we present evidence that the embryonic stem cell factor FOXD3, first identified as 'Genesis' and involved in both gastrointestinal and neural crest cell development, is implicated in GIST pathogenesis; its involvement is investigated both in vitro and in zebrafish and a mouse model of FOXD3 deficiency. Samples from a total of 58 patients with wild-type GISTs were used for molecular analyses, including Sanger sequencing, comparative genomic hybridization, and methylation analysis. Immunohistochemistry and western blot evaluation were used to assess FOXD3 expression. Additionally, we conducted in vitro functional studies in tissue samples and in transfected cells to confirm the pathogenicity of the identified genetic variants. Germline partially inactivating FOXD3 sequence variants (p.R54H and p.Ala88_Gly91del) were found in patients with isolated GISTs. Chromosome 1p loss was the most frequent chromosomal abnormality identified in tumors. In vitro experiments demonstrate the impairment of FOXD3 in the presence of those variants. Animal studies showed disruption of the GI neural network and changes in the number and distribution in the ICC. FOXD3 suppresses KIT expression in human cells; its inactivation led to an increase in ICC in zebrafish, as well as mice, providing evidence for a functional link between FOXD3 defects and KIT overexpression leading to GIST formation.

Buffy coat signatures of breast cancer risk in a prospective cohort study.

BACKGROUND: Epigenetic alterations are a near-universal feature of human malignancy and have been detected in malignant cells as well as in easily accessible specimens such as blood and urine. These findings offer promising applications in cancer detection, subtyping, and treatment monitoring. However, much of the current evidence is based on findings in retrospective studies and may reflect epigenetic patterns that have already been influenced by the onset of the disease. METHODS: Studying breast cancer, we established genome-scale DNA methylation profiles of prospectively collected buffy coat samples (n = 702) from a case-control study nested within the EPIC-Heidelberg cohort using reduced representation bisulphite sequencing (RRBS). RESULTS: We observed cancer-specific DNA methylation events in buffy coat samples. Increased DNA methylation in genomic regions associated with SURF6 and REXO1/CTB31O20.3 was linked to the length of time to diagnosis in the prospectively collected buffy coat DNA from individuals who subsequently developed breast cancer. Using machine learning methods, we piloted a DNA methylation-based classifier that predicted case-control status in a held-out validation set with 76.5% accuracy, in some cases up to 15 years before clinical diagnosis of the disease. CONCLUSIONS: Taken together, our findings suggest a model of gradual accumulation of cancer-associated DNA methylation patterns in peripheral blood, which may be detected long before clinical manifestation of cancer. Such changes may provide useful markers for risk stratification and, ultimately, personalized cancer prevention.

High-resolution transcriptomic and epigenetic profiling identifies novel regulators of COPD.

Patients with chronic obstructive pulmonary disease (COPD) are still waiting for curative treatments. Considering its environmental cause, we hypothesized that COPD will be associated with altered epigenetic signaling in lung cells. We generated genome-wide DNA methylation maps at single CpG resolution of primary human lung fibroblasts (HLFs) across COPD stages. We show that the epigenetic landscape is changed early in COPD, with DNA methylation changes occurring predominantly in regulatory regions. RNA sequencing of matched fibroblasts demonstrated dysregulation of genes involved in proliferation, DNA repair, and extracellular matrix organization. Data integration identified 110 candidate regulators of disease phenotypes that were linked to fibroblast repair processes using phenotypic screens. Our study provides high-resolution multi-omic maps of HLFs across COPD stages. We reveal novel transcriptomic and epigenetic signatures associated with COPD onset and progression and identify new candidate regulators involved in the pathogenesis of chronic lung diseases. The presence of various epigenetic factors among the candidates demonstrates that epigenetic regulation in COPD is an exciting research field that holds promise for novel therapeutic avenues for patients.

Indirect epigenetic testing identifies a diagnostic signature of cardiomyocyte DNA methylation in heart failure.

Precision-based molecular phenotyping of heart failure must overcome limited access to cardiac tissue. Although epigenetic alterations have been found to underlie pathological cardiac gene dysregulation, the clinical utility of myocardial epigenomics remains narrow owing to limited clinical access to tissue. Therefore, the current study determined whether patient plasma confers indirect phenotypic, transcriptional, and/or epigenetic alterations to ex vivo cardiomyocytes to mirror the failing human myocardium. Neonatal rat ventricular myocytes (NRVMs) and single-origin human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) and were treated with blood plasma samples from patients with dilated cardiomyopathy (DCM) and donor subjects lacking history of cardiovascular disease. Following plasma treatments, NRVMs and hiPSC-CMs underwent significant hypertrophy relative to non-failing controls, as determined via automated high-content screening. Array-based DNA methylation analysis of plasma-treated hiPSC-CMs and cardiac biopsies uncovered robust, and conserved, alterations in cardiac DNA methylation, from which 100 sites were validated using an independent cohort. Among the CpG sites identified, hypo-methylation of the ATG promoter was identified as a diagnostic marker of HF, wherein cg03800765 methylation (AUC = 0.986, P < 0.0001) was found to out-perform circulating NT-proBNP levels in differentiating heart failure. Taken together, these findings support a novel approach of indirect epigenetic testing in human HF.

Harnessing the MYB-dependent TAL1 5'super-enhancer for targeted therapy in T-ALL.

The acquisition of genetic abnormalities engendering oncogene dysregulation underpins cancer development. Certain proto-oncogenes possess several dysregulation mechanisms, yet how each mechanism impacts clinical outcome is unclear. Using T-cell acute lymphoblastic leukemia (T-ALL) as an example, we show that patients harboring 5'super-enhancer (5'SE) mutations of the TAL1 oncogene identifies a specific patient subgroup with poor prognosis irrespective of the level of oncogene dysregulation. Remarkably, the MYB dependent oncogenic 5'SE can be targeted using Mebendazole to induce MYB protein degradation and T-ALL cell death. Of note Mebendazole treatment demonstrated efficacy in vivo in T-ALL preclinical models. Our work provides proof of concept that within a specific oncogene driven cancer, the mechanism of oncogene dysregulation rather than the oncogene itself can identify clinically distinct patient subgroups and pave the way for future super-enhancer targeting therapy.

Global hypomethylation in childhood asthma identified by genome-wide DNA-methylation sequencing preferentially affects enhancer regions.

BACKGROUND: Childhood asthma is a result of a complex interaction of genetic and environmental components causing epigenetic and immune dysregulation, airway inflammation and impaired lung function. Although different microarray based EWAS studies have been conducted, the impact of epigenetic regulation in asthma development is still widely unknown. We have therefore applied unbiased whole genome bisulfite sequencing (WGBS) to characterize global DNA-methylation profiles of asthmatic children compared to healthy controls. METHODS: Peripheral blood samples of 40 asthmatic and 42 control children aged 5-15 years from three birth cohorts were sequenced together with paired cord blood samples. Identified differentially methylated regions (DMRs) were categorized in genotype-associated, cell-type-dependent, or prenatally primed. Network analysis and subsequent natural language processing of DMR-associated genes was complemented by targeted analysis of functional translation of epigenetic regulation on the transcriptional and protein level. RESULTS: In total, 158 DMRs were identified in asthmatic children compared to controls of which 37% were related to the eosinophil content. A global hypomethylation was identified affecting predominantly enhancer regions and regulating key immune genes such as IL4, IL5RA, and EPX. These DMRs were confirmed in n = 267 samples and could be linked to aberrant gene expression. Out of the 158 DMRs identified in the established phenotype, 56 were perturbed already at birth and linked, at least in part, to prenatal influences such as tobacco smoke exposure or phthalate exposure. CONCLUSION: This is the first epigenetic study based on whole genome sequencing to identify marked dysregulation of enhancer regions as a hallmark of childhood asthma.

TAL1 activation in T-cell acute lymphoblastic leukemia: a novel oncogenic 3' neo-enhancer.

T-cell acute lymphocytic leukemia protein 1 (TAL1) is one of the most frequently deregulated oncogenes in T-cell acute lymphoblastic leukemia (T-ALL). Its deregulation can occur through diverse cis-alterations, including SIL-TAL1 microdeletions, translocations with T-cell Receptor loci, and more recently described upstream intergenic non-coding mutations. These mutations consist of recurrent focal microinsertions that create an oncogenic neo-enhancer accompanied by activating epigenetic marks. This observation laid the groundwork for an innovative paradigm concerning the activation of proto-oncogenes via genomic alterations of non-coding intergenic regions. However, for the majority of T-ALL expressing TAL1 (TAL1+), the deregulation mechanism remains 'unresolved'. We took advantage of H3K27ac and H3K4me3 chromatin immunoprecipitation sequencing data of eight cases of T-ALL, including five TAL1+ cases. We identified a putative novel oncogenic neo-enhancer downstream of TAL1 in an unresolved monoallelic TAL1+ case. A rare but recurrent somatic heterozygous microinsertion within this region creates a de novo binding site for MYB transcription factor. Here we demonstrate that this mutation leads to increased enhancer activity, gain of active epigenetic marks, and TAL1 activation via recruitment of MYB. These results highlight the diversity of non-coding mutations that can drive oncogene activation.