DNA methylation landscape of tumor-associated macrophages reveals pathways, transcription factors and prognostic value relevant to triple-negative breast cancer patients

The accumulation of myeloid cells, particularly tumor-associated macrophages (TAMs), characterizes the tumor microenvironment (TME) of many solid cancers, including breast cancer. Compared to healthy tissue-resident macrophages, TAMs acquire distinct transcriptomes and tumor-promoting functions by largely unknown mechanisms. Here, we hypothesize the involvement of TME signaling and subsequent epigenetic reprogramming of TAMs. Using the 4T1 mouse model of triple-negative breast cancer, we demonstrate that the presence of cancer cells significantly alters the DNA methylation landscape of macrophages and, to a lesser extent, bone marrow-derived monocytes (BMDMs). TAM methylomes, dissected into BMDM-originating and TAM-specific epigenetic programs, implicated transcription factors (TFs) and signaling pathways involved in TAM reprogramming, correlated with cancer-specific gene expression patterns. Utilizing published single-cell gene expression data, we linked microenvironmentally-derived signals to the cancer-specific DNA methylation landscape of TAMs. These integrative analyses highlighted the role of altered cytokine production in the TME (eg, TGF-β, IFN-γ and CSF1) on the induction of specific TFs (eg, FOSL2, STAT1 and RUNX3) responsible for the epigenetic reprogramming of TAMs. DNA methylation deconvolution identified a TAM-specific signature associated with the identified signaling pathways and TFs, corresponding with severe tumor grade and poor prognosis of breast cancer patients. Similarly, immunosuppressive TAM functions were identified, such as induction of the immune inhibitory receptor-ligand PD-L1 by DNA hypomethylation of Cd274. Collectively, these results provide strong evidence that the epigenetic landscapes of macrophages and monocytes are perturbed by the presence of breast cancer, pointing to molecular mechanisms of TAM reprogramming, impacting patient outcomes.

DNA methylation biomarkers for noninvasive detection of triple-negative breast cancer using liquid biopsy

Noninvasive detection of aberrant DNA methylation could provide invaluable biomarkers for earlier detection of triple-negative breast cancer (TNBC) which could help clinicians with easier and more efficient treatment options. We evaluated genome-wide DNA methylation data derived from TNBC and normal breast tissues, peripheral blood of TNBC cases and controls and reference samples of sorted blood and mammary cells. Differentially methylated regions (DMRs) between TNBC and normal breast tissues were stringently selected, verified and externally validated. A machine-learning algorithm was applied to select the top DMRs, which then were evaluated on plasma-derived circulating cell-free DNA (cfDNA) samples of TNBC patients and healthy controls. We identified 23 DMRs accounting for the methylation profile of blood cells and reference mammary cells and then selected six top DMRs for cfDNA analysis. We quantified un-/methylated copies of these DMRs by droplet digital PCR analysis in a plasma test set from TNBC patients and healthy controls and confirmed our findings obtained on tissues. Differential cfDNA methylation was confirmed in an independent validation set of plasma samples. A methylation score combining signatures of the top three DMRs overlapping with the SPAG6, LINC10606 and TBCD/ZNF750 genes had the best capability to discriminate TNBC patients from controls (AUC = 0.78 in the test set and AUC = 0.74 in validation set). Our findings demonstrate the usefulness of cfDNA-based methylation signatures as noninvasive liquid biopsy markers for the diagnosis of TNBC.

Dynamic DNA methylation reveals novel cis-regulatory elements in mouse hematopoiesis

Differentiation of hematopoietic stem and progenitor cells to terminally differentiated immune cells is accompanied by large-scale remodeling of the DNA methylation landscape. Although significant insights into the molecular mechanisms of hematopoietic tissue regeneration were derived from mouse models, profiling of DNA methylation has been hampered by high cost or low resolution using available methods. The recent development of the Infinium Mouse Methylation BeadChip (MMBC) array facilitates methylation profiling of the mouse genome at a single CpG resolution at affordable cost. We extended the RnBeads package to provide a computational framework for the analysis of MMBC data. This framework was applied to a newly generated reference map of mouse hematopoiesis encompassing nine different cell types. Analysis of dynamically regulated CpG sites showed progressive and unidirectional DNA methylation changes from hematopoietic stem and progenitor cells to differentiated hematopoietic cells and allowed the identification of lineage- and cell type-specific DNA methylation programs. Comparison with previously published catalogs of cis-regulatory elements (CREs) revealed 12,856 novel putative CREs that were dynamically regulated by DNA methylation (mdCREs). These mdCREs were predominantly associated with patterns of cell type-specific DNA hypomethylation and could be identified as epigenetic control regions regulating the expression of key hematopoietic genes during differentiation. In summary, we established an analysis pipeline for MMBC data sets and provide a DNA methylation atlas of mouse hematopoiesis. This resource allowed us to identify novel putative CREs involved in hematopoiesis and will serve as a platform to study epigenetic regulation of normal and malignant hematopoiesis.

Control of APOBEC3B induction and cccDNA decay by NF-κB and miR-138-5p

Background & Aims

Immune-mediated induction of cytidine deaminase APOBEC3B (A3B) expression leads to HBV covalently closed circular DNA (cccDNA) decay. Here, we aimed to decipher the signalling pathway(s) and regulatory mechanism(s) involved in A3B induction and related HBV control.

Methods

Differentiated HepaRG cells (dHepaRG) knocked-down for NF-κB signalling components, transfected with siRNA or micro RNAs (miRNA), and primary human hepatocytes ± HBV or HBVΔX or HBV-RFP, were treated with lymphotoxin beta receptor (LTβR)-agonist (BS1). The biological outcomes were analysed by reverse transcriptase-qPCR, immunoblotting, luciferase activity, chromatin immune precipitation, electrophoretic mobility-shift assay, targeted-bisulfite-, miRNA-, RNA-, genome-sequencing, and mass-spectrometry.

Results

We found that canonical and non-canonical NF-κB signalling pathways are mandatory for A3B induction and anti-HBV effects. The degree of immune-mediated A3B production is independent of A3B promoter demethylation but is controlled post-transcriptionally by the miRNA 138-5p expression (hsa-miR-138-5p), promoting A3B mRNA decay. Hsa-miR-138-5p over-expression reduced A3B levels and its antiviral effects. Of note, established infection inhibited BS1-induced A3B expression through epigenetic modulation of A3B promoter. Twelve days of treatment with a LTβR-specific agonist BS1 is sufficient to reduce the cccDNA pool by 80% without inducing significant damages to a subset of cancer-related host genes. Interestingly, the A3B-mediated effect on HBV is independent of the transcriptional activity of cccDNA as well as on rcDNA synthesis.

Conclusions

Altogether, A3B represents the only described enzyme to target both transcriptionally active and inactive cccDNA. Thus, inhibiting hsa-miR-138-5p expression should be considered in the combinatorial design of new therapies against HBV, especially in the context of immune-mediated A3B induction.

Lay summary

Immune-mediated induction of cytidine deaminase APOBEC3B is transcriptionally regulated by NF-κB signalling and post-transcriptionally downregulated by hsa-miR-138-5p expression, leading to cccDNA decay. Timely controlled APOBEC3B-mediated cccDNA decay occurs independently of cccDNA transcriptional activity and without damage to a subset of cancer-related genes. Thus, APOBEC3B-mediated cccDNA decay could offer an efficient therapeutic alternative to target hepatitis B virus chronic infection.

•

Impairment of NF-κB signalling prevents APOBEC3B induction and cccDNA decay.

•

APOBEC3B is post-transcriptionally regulated by the hsa-miR-138-5p.

•

Over-expression of the hsa-miR-138-5p inhibits APOBEC3B expression and cccDNA decay.

•

A3B timely induces cccDNA decay without damage to cancer-related genes.

•

APOBEC3B-mediated cccDNA decay is independent of cccDNA transcriptional activity.

Hotspot DNMT3A mutations in clonal hematopoiesis and acute myeloid leukemia sensitize cells to azacytidine via viral mimicry response

Somatic mutations in DNA methyltransferase 3A (DNMT3A) are among the most frequent alterations in clonal hematopoiesis (CH) and acute myeloid leukemia (AML), with a hotspot in exon 23 at arginine 882 (DNMT3A^R882). Here, we demonstrate that DNMT3A^R882H-dependent CH and AML cells are specifically susceptible to the hypomethylating agent azacytidine (AZA). Addition of AZA to chemotherapy prolonged AML survival solely in individuals with DNMT3A^R882 mutations, suggesting its potential as a predictive marker for AZA response. AML and CH mouse models confirmed AZA susceptibility specifically in DNMT3A^R882H-expressing cells. Hematopoietic stem cells (HSCs) and progenitor cells expressing DNMT3A^R882H exhibited cell autonomous viral mimicry response as a result of focal DNA hypomethylation at retrotransposon sequences. Administration of AZA boosted hypomethylation of retrotransposons specifically in DNMT3A^R882H-expressing cells and maintained elevated levels of canonical interferon-stimulated genes (ISGs), thus leading to suppressed protein translation and increased apoptosis.

Methrix: an R/bioconductor package for systematic aggregation and analysis of bisulfite sequencing data

MOTIVATION

Whole genome bisulfite sequencing (WGBS), measures DNA methylation at base pair resolution resulting in large bedGraph like coverage files. Current options for processing such files are hindered by discrepancies in file format specification, speed and memory requirements.

RESULTS

We developed methrix, an R package, which provides a toolset for systematic analysis of large datasets. Core functionality of the package includes a comprehensive bedGraph or similar tab-separated text file reader - which summarizes methylation calls based on annotated reference indices, infers and collapses strands, and handles uncovered reference CpG sites while facilitating a flexible input file format specification. Additional optimized functions for quality control filtering, sub-setting, and visualization allow user-friendly and effective processing of WGBS results. Easy integration with tools for differentially methylated region (DMR) calling and annotation further eases the analysis of genome-wide methylation data. Overall, methrix enriches established WGBS workflows by bringing together computational efficiency and versatile functionality.

AVAILABILITY AND IMPLEMENTATION

Methrix is implemented as an R package, made available under MIT license at https://github.com/CompEpigen/methrix and can be installed from the Bioconductor repository.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

International Consensus Definition of DNA Methylation Subgroups in Juvenile Myelomonocytic Leukemia

PURPOSE

Known clinical and genetic markers have limitations in predicting disease course and outcome in juvenile myelomonocytic leukemia (JMML). DNA methylation patterns in JMML have correlated with outcome across multiple studies, suggesting it as a biomarker to improve patient stratification. However, standardized approaches to classify JMML on the basis of DNA methylation patterns are lacking. We, therefore, sought to define an international consensus for DNA methylation subgroups in JMML and develop classification methods for clinical implementation.

EXPERIMENTAL DESIGN

Published DNA methylation data from 255 patients with JMML were used to develop and internally validate a classifier model. Accuracy across platforms (EPIC-arrays and MethylSeq) was tested using a technical validation cohort (32 patients). The suitability of both methods for single-patient classification was demonstrated using an independent cohort (47 patients).

RESULTS

Analysis of pooled, published data established three DNA methylation subgroups as a de facto standard. Unfavorable prognostic parameters (PTPN11 mutation, elevated fetal hemoglobin, and older age) were significantly enriched in the high methylation (HM) subgroup. A classifier was then developed that predicted subgroups with 98% accuracy across different technological platforms. Applying the classifier to an independent validation cohort confirmed an association of HM with secondary mutations, high relapse incidence, and inferior overall survival (OS), while the low methylation subgroup was associated with a favorable disease course. Multivariable analysis established DNA methylation subgroups as the only significant factor predicting OS.

CONCLUSIONS

This study provides an international consensus definition for DNA methylation subgroups in JMML. We developed and validated methods which will facilitate the design of risk-stratified clinical trials in JMML.

AmpliconDesign - an interactive web server for the design of high-throughput targeted DNA methylation assays

Targeted analysis of DNA methylation patterns based on bisulfite-treated genomic DNA (BT-DNA) is considered as a gold-standard for epigenetic biomarker development. Existing software tools facilitate primer design, primer quality control or visualization of primer localization. However, high-throughput design of primers for BT-DNA amplification is hampered by limits in throughput and functionality of existing tools, requiring users to repeatedly perform specific tasks manually. Consequently, the design of PCR primers for BT-DNA remains a tedious and time-consuming process. To bridge this gap, we developed AmpliconDesign, a webserver providing a scalable and user-friendly platform for the design and analysis of targeted DNA methylation studies based on BT-DNA, e.g. deep amplicon bisulfite sequencing (ampBS-seq) or EpiTYPER MassArray. Core functionality of the web server includes high-throughput primer design and binding site validation based on in silico bisulfite-converted DNA sequences, prediction of fragmentation patterns for EpiTYPER MassArray, an interactive quality control as well as a streamlined analysis workflow for ampBS-seq.

RAS-pathway mutation patterns define epigenetic subclasses in juvenile myelomonocytic leukemia

Juvenile myelomonocytic leukemia (JMML) is an aggressive myeloproliferative disorder of early childhood characterized by mutations activating RAS signaling. Established clinical and genetic markers fail to fully recapitulate the clinical and biological heterogeneity of this disease. Here we report DNA methylome analysis and mutation profiling of 167 JMML samples. We identify three JMML subgroups with unique molecular and clinical characteristics. The high methylation group (HM) is characterized by somatic PTPN11 mutations and poor clinical outcome. The low methylation group is enriched for somatic NRAS and CBL mutations, as well as for Noonan patients, and has a good prognosis. The intermediate methylation group (IM) shows enrichment for monosomy 7 and somatic KRAS mutations. Hypermethylation is associated with repressed chromatin, genes regulated by RAS signaling, frequent co-occurrence of RAS pathway mutations and upregulation of DNMT1 and DNMT3B, suggesting a link between activation of the DNA methylation machinery and mutational patterns in JMML.