Systematic evaluation of cell-type deconvolution pipelines for sequencing-based bulk DNA methylomes.

DNA methylation analysis by sequencing is becoming increasingly popular, yielding methylomes at single-base pair and single-molecule resolution. It has tremendous potential for cell-type heterogeneity analysis using intrinsic read-level information. Although diverse deconvolution methods were developed to infer cell-type composition based on bulk sequencing-based methylomes, systematic evaluation has not been performed yet. Here, we thoroughly benchmark six previously published methods: Bayesian epiallele detection, DXM, PRISM, csmFinder+coMethy, ClubCpG and MethylPurify, together with two array-based methods, MeDeCom and Houseman, as a comparison group. Sequencing-based deconvolution methods consist of two main steps, informative region selection and cell-type composition estimation, thus each was individually assessed. With this elaborate evaluation, we aimed to establish which method achieves the highest performance in different scenarios of synthetic bulk samples. We found that cell-type deconvolution performance is influenced by different factors depending on the number of cell types within the mixture. Finally, we propose a best-practice deconvolution strategy for sequencing data and point out limitations that need to be handled. Array-based methods-both reference-based and reference-free-generally outperformed sequencing-based methods, despite the absence of read-level information. This implies that the current sequencing-based methods still struggle with correctly identifying cell-type-specific signals and eliminating confounding methylation patterns, which needs to be handled in future studies.

The age of the bone marrow microenvironment influences B-cell acute lymphoblastic leukemia progression via CXCR5-CXCL13.

B-cell acute lymphoblastic leukemia (B-ALL) occurs most commonly in children, whereas chronic myeloid leukemia is more frequent in adults. The myeloid bias of hematopoiesis in elderly individuals has been considered causative, but the age of the bone marrow microenvironment (BMM) may be contributory. Using various murine models of B-ALL in young vs old mice, we recapitulated B-ALL preponderance in children vs adults. We showed differential effects of young vs old BM macrophages on B-ALL cell function. Molecular profiling using RNA- and ATAC-sequencing revealed pronounced differences in young vs old BMM-derived macrophages and enrichment for gene sets associated with inflammation. In concordance with the role of C-X-C motif chemokine (CXCL) 13 for disease-associated B-cell chemoattraction, we found CXCL13 to be highly expressed in young macrophages on a translational compared with a transcriptional level. Inhibition of CXCL13 in BM macrophages impaired leukemia cell migration and decreased the proliferation of cocultured B-ALL cells, whereas recombinant CXCL13 increased pAKT and B-ALL cell expansion. Pretreatment of B-ALL-initiating cells with CXCL13 accelerated B-ALL progression. Deficiency of Cxcr5, the receptor for CXCL13, on B-ALL-initiating cells prolonged murine survival, whereas high expression of CXCR5 in pediatric B-ALL may predict central nervous system relapse. CXCL13 staining was increased in bone sections from pediatric compared with adult patients with B-ALL. Taken together, our study shows that the age of the BMM and, in particular, BM macrophages influence the leukemia phenotype. The CXCR5-CXCL13 axis may act as prognostic marker and an attractive novel target for the treatment of B-ALL.

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.

ID3 promotes homologous recombination via non-transcriptional and transcriptional mechanisms and its loss confers sensitivity to PARP inhibition.

The inhibitor of DNA-binding 3 (ID3) is a transcriptional regulator that limits interaction of basic helix-loop-helix transcription factors with their target DNA sequences. We previously reported that ID3 loss is associated with mutational signatures linked to DNA repair defects. Here we demonstrate that ID3 exhibits a dual role to promote DNA double-strand break (DSB) repair, particularly homologous recombination (HR). ID3 interacts with the MRN complex and RECQL helicase to activate DSB repair and it facilitates RAD51 loading and downstream steps of HR. In addition, ID3 promotes the expression of HR genes in response to ionizing radiation by regulating both chromatin accessibility and activity of the transcription factor E2F1. Consistently, analyses of TCGA cancer patient data demonstrate that low ID3 expression is associated with impaired HR. The loss of ID3 leads to sensitivity of tumor cells to PARP inhibition, offering new therapeutic opportunities in ID3-deficient tumors.

Histone H3K27 demethylase KDM6A is an epigenetic gatekeeper of mTORC1 signalling in cancer.

OBJECTIVE: Large-scale genome sequencing efforts of human tumours identified epigenetic modifiers as one of the most frequently mutated gene class in human cancer. However, how these mutations drive tumour development and tumour progression are largely unknown. Here, we investigated the function of the histone demethylase KDM6A in gastrointestinal cancers, such as liver cancer and pancreatic cancer. DESIGN: Genetic alterations as well as expression analyses of KDM6A were performed in patients with liver cancer. Genetic mouse models of liver and pancreatic cancer coupled with Kdm6a-deficiency were investigated, transcriptomic and epigenetic profiling was performed, and in vivo and in vitro drug treatments were conducted. RESULTS: KDM6A expression was lost in 30% of patients with liver cancer. Kdm6a deletion significantly accelerated tumour development in murine liver and pancreatic cancer models. Kdm6a-deficient tumours showed hyperactivation of mTORC1 signalling, whereas endogenous Kdm6a re-expression by inducible RNA-interference in established Kdm6a-deficient tumours diminished mTORC1 activity resulting in attenuated tumour progression. Genome-wide transcriptional and epigenetic profiling revealed direct binding of Kdm6a to crucial negative regulators of mTORC1, such as Deptor, and subsequent transcriptional activation by epigenetic remodelling. Moreover, in vitro and in vivo genetic epistasis experiments illustrated a crucial function of Deptor and mTORC1 in Kdm6a-dependent tumour suppression. Importantly, KDM6A expression in human tumours correlates with mTORC1 activity and KDM6A-deficient tumours exhibit increased sensitivity to mTORC1 inhibition. CONCLUSION: KDM6A is an important tumour suppressor in gastrointestinal cancers and acts as an epigenetic toggle for mTORC1 signalling. Patients with KDM6A-deficient tumours could benefit of targeted therapy focusing on mTORC1 inhibition.

Epigenetic Modulation of Radiation-Induced Diacylglycerol Kinase Alpha Expression Prevents Pro-Fibrotic Fibroblast Response.

Radiotherapy, a common component in cancer treatment, can induce adverse effects including fibrosis in co-irradiated tissues. We previously showed that differential DNA methylation at an enhancer of diacylglycerol kinase alpha (DGKA) in normal dermal fibroblasts is associated with radiation-induced fibrosis. After irradiation, the transcription factor EGR1 is induced and binds to the hypomethylated enhancer, leading to increased DGKA and pro-fibrotic marker expression. We now modulated this DGKA induction by targeted epigenomic and genomic editing of the DGKA enhancer and administering epigenetic drugs. Targeted DNA demethylation of the DGKA enhancer in HEK293T cells resulted in enrichment of enhancer-related histone activation marks and radiation-induced DGKA expression. Mutations of the EGR1-binding motifs decreased radiation-induced DGKA expression in BJ fibroblasts and caused dysregulation of multiple fibrosis-related pathways. EZH2 inhibitors (GSK126, EPZ6438) did not change radiation-induced DGKA increase. Bromodomain inhibitors (CBP30, JQ1) suppressed radiation-induced DGKA and pro-fibrotic marker expression. Similar drug effects were observed in donor-derived fibroblasts with low DNA methylation. Overall, epigenomic manipulation of DGKA expression may offer novel options for a personalized treatment to prevent or attenuate radiotherapy-induced fibrosis.

Globally altered epigenetic landscape and delayed osteogenic differentiation in H3.3-G34W-mutant giant cell tumor of bone.

The neoplastic stromal cells of giant cell tumor of bone (GCTB) carry a mutation in H3F3A, leading to a mutant histone variant, H3.3-G34W, as a sole recurrent genetic alteration. We show that in patient-derived stromal cells H3.3-G34W is incorporated into the chromatin and associates with massive epigenetic alterations on the DNA methylation, chromatin accessibility and histone modification level, that can be partially recapitulated in an orthogonal cell line system by the introduction of H3.3-G34W. These epigenetic alterations affect mainly heterochromatic and bivalent regions and provide possible explanations for the genomic instability, as well as the osteolytic phenotype of GCTB. The mutation occurs in differentiating mesenchymal stem cells and associates with an impaired osteogenic differentiation. We propose that the observed epigenetic alterations reflect distinct differentiation stages of H3.3 WT and H3.3 MUT stromal cells and add to H3.3-G34W-associated changes.