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.

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.

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.

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.

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.

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, subsetting 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.

Altered enhancer-promoter interaction leads to MNX1 expression in pediatric acute myeloid leukemia with t(7;12)(q36;p13).

ABSTRACT: Acute myeloid leukemia (AML) with the t(7;12)(q36;p13) translocation occurs only in very young children and has a poor clinical outcome. The expected oncofusion between break point partners (motor neuron and pancreas homeobox 1 [MNX1] and ETS variant transcription factor 6 [ETV6]) has only been reported in a subset of cases. However, a universal feature is the strong transcript and protein expression of MNX1, a homeobox transcription factor that is normally not expressed in hematopoietic cells. Here, we map the translocation break points on chromosomes 7 and 12 in affected patients to a region proximal to MNX1 and either introns 1 or 2 of ETV6. The frequency of MNX1 overexpression in pediatric AML is 2.4% and occurs predominantly in t(7;12)(q36;p13) AML. Chromatin interaction assays in a t(7;12)(q36;p13) induced pluripotent stem cell line model unravel an enhancer-hijacking event that explains MNX1 overexpression in hematopoietic cells. Our data suggest that enhancer hijacking may be a more widespread consequence of translocations in which no oncofusion product was identified, including t(1;3) or t(4;12) AML.

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.

Club cells employ regeneration mechanisms during lung tumorigenesis.

The high plasticity of lung epithelial cells, has for many years, confounded the correct identification of the cell-of-origin of lung adenocarcinoma (LUAD), one of the deadliest malignancies worldwide. Here, we employ lineage-tracing mouse models to investigate the cell of origin of Eml4-Alk LUAD, and show that Club and Alveolar type 2 (AT2) cells give rise to tumours. We focus on Club cell originated tumours and find that Club cells experience an epigenetic switch by which they lose their lineage fidelity and gain an AT2-like phenotype after oncogenic transformation. Single-cell transcriptomic analyses identified two trajectories of Club cell evolution which are similar to the ones used during lung regeneration, suggesting that lung epithelial cells leverage on their plasticity and intrinsic regeneration mechanisms to give rise to a tumour. Together, this study highlights the role of Club cells in LUAD initiation, identifies the mechanism of Club cell lineage infidelity, confirms the presence of these features in human tumours, and unveils key mechanisms conferring LUAD heterogeneity.

Global DNA Methylation Analysis of Cancer-Associated Fibroblasts Reveals Extensive Epigenetic Rewiring Linked with RUNX1 Upregulation in Breast Cancer Stroma.

Cancer cells recruit and rewire normal fibroblasts in their microenvironment to become protumorigenic cancer-associated fibroblasts (CAF). These CAFs are genomically stable, yet their transcriptional programs are distinct from those of their normal counterparts. Transcriptional regulation plays a major role in this reprogramming, but the extent to which epigenetic modifications of DNA also contribute to the rewiring of CAF transcription is not clear. Here we address this question by dissecting the epigenetic landscape of breast CAFs. Applying tagmentation-based whole-genome bisulfite sequencing in a mouse model of breast cancer, we found that fibroblasts undergo massive DNA methylation changes as they transition into CAFs. Transcriptional and epigenetic analyses revealed RUNX1 as a potential mediator of this process and identified a RUNX1-dependent stromal gene signature. Coculture and mouse models showed that both RUNX1 and its stromal signature are induced as normal fibroblasts transition into CAFs. In breast cancer patients, RUNX1 was upregulated in CAFs, and expression of the RUNX1 signature was associated with poor disease outcome, highlighting the relevance of these findings to human disease. This work presents a comprehensive genome-wide map of DNA methylation in CAFs and reveals a previously unknown facet of the dynamic plasticity of the stroma. SIGNIFICANCE: The first genome-wide map of DNA methylation in breast cancer-associated fibroblasts unravels a previously unknown facet of the dynamic plasticity of the stroma, with far-reaching therapeutic implications.

The antileukemic activity of decitabine upon PML/RARA-negative AML blasts is supported by all-trans retinoic acid: in vitro and in vivo evidence for cooperation.

The prognosis of AML patients with adverse genetics, such as a complex, monosomal karyotype and TP53 lesions, is still dismal even with standard chemotherapy. DNA-hypomethylating agent monotherapy induces an encouraging response rate in these patients. When combined with decitabine (DAC), all-trans retinoic acid (ATRA) resulted in an improved response rate and longer overall survival in a randomized phase II trial (DECIDER; NCT00867672). The molecular mechanisms governing this in vivo synergism are unclear. We now demonstrate cooperative antileukemic effects of DAC and ATRA on AML cell lines U937 and MOLM-13. By RNA-sequencing, derepression of >1200 commonly regulated transcripts following the dual treatment was observed. Overall chromatin accessibility (interrogated by ATAC-seq) and, in particular, at motifs of retinoic acid response elements were affected by both single-agent DAC and ATRA, and enhanced by the dual treatment. Cooperativity regarding transcriptional induction and chromatin remodeling was demonstrated by interrogating the HIC1, CYP26A1, GBP4, and LYZ genes, in vivo gene derepression by expression studies on peripheral blood blasts from AML patients receiving DAC + ATRA. The two drugs also cooperated in derepression of transposable elements, more effectively in U937 (mutated TP53) than MOLM-13 (intact TP53), resulting in a "viral mimicry" response. In conclusion, we demonstrate that in vitro and in vivo, the antileukemic and gene-derepressive epigenetic activity of DAC is enhanced by ATRA.

Transcriptional regulation of the thymus master regulator Foxn1.

FOXN1 is a transcription factor critical for the development of both thymic epithelial cell (TEC) and hair follicle cell (HFC) compartments. However, mechanisms controlling its expression remain poorly understood. To address this question, we performed thorough analyses of the evolutionary conservation and chromatin status of the Foxn1 locus in different tissues and states and identified several putative cis-regulatory regions unique to TECs versus HFCs. Furthermore, experiments using genetically modified mice with specific deletions in the Foxn1 locus and additional bioinformatic analyses helped us identify key regions and transcription factors involved in either positive or negative regulation of Foxn1 in both TECs and HFCs. Specifically, we identified SIX1 and FOXN1 itself as key factors inducing Foxn1 expression in embryonic and neonatal TECs. Together, our data provide important mechanistic insights into the transcriptional regulation of the Foxn1 gene in TEC versus HFC and highlight the role of FOXN1 in its autoregulation.

Epigenetic reprogramming of airway macrophages promotes polarization and inflammation in muco-obstructive lung disease.

Lung diseases, such as cystic fibrosis and COPD, are characterized by mucus obstruction and chronic airway inflammation, but their mechanistic link remains poorly understood. Here, we focus on the function of the mucostatic airway microenvironment on epigenetic reprogramming of airway macrophages (AM) and resulting transcriptomic and phenotypical changes. Using a mouse model of muco-obstructive lung disease (Scnn1b-transgenic), we identify epigenetically controlled, differentially regulated pathways and transcription factors involved in inflammatory responses and macrophage polarization. Functionally, AMs from Scnn1b-transgenic mice have reduced efferocytosis and phagocytosis, and excessive inflammatory responses upon lipopolysaccharide challenge, mediated through enhanced Irf1 function and expression. Ex vivo stimulation of wild-type AMs with native mucus impairs efferocytosis and phagocytosis capacities. In addition, mucus induces gene expression changes, comparable with those observed in AMs from Scnn1b-transgenic mice. Our data show that mucostasis induces epigenetic reprogramming of AMs, leading to changes favoring tissue damage and disease progression. Targeting these altered AMs may support therapeutic approaches in patients with muco-obstructive lung diseases.

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.

Epigenetic analysis of patients with T-ALL identifies poor outcomes and a hypomethylating agent-responsive subgroup.

Adult "T cell" acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy that is associated with poor outcomes, requiring additional therapeutic options. The DNA methylation landscapes of adult T-ALL remain undercharacterized. Here, we systematically analyzed the DNA methylation profiles of normal thymic-sorted T cell subpopulations and 143 primary adult T-ALLs as part of the French GRAALL 2003-2005 trial. Our results indicated that T-ALL is epigenetically heterogeneous consisting of five subtypes (C1-C5), which were either associated with co-occurring DNA methyltransferase 3 alpha (DNMT3A)/isocitrate dehydrogenase [NADP(+)] 2 (IDH2) mutations (C1), TAL bHLH transcription factor 1, erythroid differentiation factor (TAL1) deregulation (C2), T cell leukemia homeobox 3 (TLX3) (C3), TLX1/in cis-homeobox A9 (HOXA9) (C4), or in trans-HOXA9 overexpression (C5). Integrative analysis of DNA methylation and gene expression identified potential cluster-specific oncogenes and tumor suppressor genes. In addition to an aggressive hypomethylated subgroup (C1), our data identified an unexpected subset of hypermethylated T-ALL (C5) associated with poor outcome and primary therapeutic response. Using mouse xenografts, we demonstrated that hypermethylated T-ALL samples exhibited therapeutic responses to the DNA hypomethylating agent 5-azacytidine, which significantly (survival probability; P = 0.001 for C3, 0.01 for C4, and 0.0253 for C5) delayed tumor progression. These findings suggest that epigenetic-based therapies may provide an alternative treatment option in hypermethylated T-ALL.

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.