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.

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.

Integrative analysis reveals early and distinct genetic and epigenetic changes in intraductal papillary and tubulopapillary cholangiocarcinogenesis.

OBJECTIVE: A detailed understanding of the molecular alterations in different forms of cholangiocarcinogenesis is crucial for a better understanding of cholangiocarcinoma (CCA) and may pave the way to early diagnosis and better treatment options. DESIGN: We analysed a clinicopathologically well-characterised patient cohort (n=54) with high-grade intraductal papillary (IPNB) or tubulopapillary (ITPN) neoplastic precursor lesions of the biliary tract and correlated the results with an independent non-IPNB/ITPN associated CCA cohort (n=294). The triplet sample set of non-neoplastic biliary epithelium, precursor and invasive CCA was analysed by next generation sequencing, DNA copy number and genome-wide methylation profiling. RESULTS: Patients with invasive CCA arising from IPNB/ITPN had better prognosis than patients with CCA not associated with IPNB/ITPN. ITPN was localised mostly intrahepatic, whereas IPNB was mostly of extrahepatic origin. IPNB/ITPN were equally associated with small-duct and large-duct type intrahepatic CCA. IPNB exhibited mutational profiles of extrahepatic CCA, while ITPN had significantly fewer mutations. Most mutations were shared between precursor lesions and corresponding invasive CCA but ROBO2 mutations occurred exclusively in invasive CCA and CTNNB1 mutations were mainly present in precursor lesions. In addition, IPNB and ITPN differed in their DNA methylation profiles and analyses of latent methylation components suggested that IPNB and ITPN may have different cells-of-origin. CONCLUSION: Integrative analysis revealed that IPNB and ITPN harbour distinct early genetic alterations, IPNB are enriched in mutations typical for extrahepatic CCA, whereas ITPN exhibited few genetic alterations and showed distinct epigenetic profiles. In conclusion, IPNB/ITPN may represent a distinctive, intermediate form of intrahepatic and extrahepatic cholangiocarcinogenesis.

Identification of NHLRC1 as a Novel AKT Activator from a Lung Cancer Epigenome-Wide Association Study (EWAS).

Changes in DNA methylation identified by epigenome-wide association studies (EWAS) have been recently linked to increased lung cancer risk. However, the cellular effects of these differentially methylated positions (DMPs) are often unclear. Therefore, we investigated top differentially methylated positions identified from an EWAS study. This included a putative regulatory region of NHLRC1. Hypomethylation of this gene was recently linked with decreased survival rates in lung cancer patients. HumanMethylation450 BeadChip array (450K) analysis was performed on 66 lung cancer case-control pairs from the European Prospective Investigation into Cancer and Nutrition Heidelberg lung cancer EWAS (EPIC HD) cohort. DMPs identified in these pre-diagnostic blood samples were then investigated for differential DNA methylation in lung tumor versus adjacent normal lung tissue from The Cancer Genome Atlas (TCGA) and replicated in two independent lung tumor versus adjacent normal tissue replication sets with MassARRAY. The EPIC HD top hypermethylated DMP cg06646708 was found to be a hypomethylated region in multiple data sets of lung tumor versus adjacent normal tissue. Hypomethylation within this region caused increased mRNA transcription of the closest gene NHLRC1 in lung tumors. In functional assays, we demonstrate attenuated proliferation, viability, migration, and invasion upon NHLRC1 knock-down in lung cancer cells. Furthermore, diminished AKT phosphorylation at serine 473 causing expression of pro-apoptotic AKT-repressed genes was detected in these knock-down experiments. In conclusion, this study demonstrates the powerful potential for discovery of novel functional mechanisms in oncogenesis based on EWAS DNA methylation data. NHLRC1 holds promise as a new prognostic biomarker for lung cancer survival and prognosis, as well as a target for novel treatment strategies in lung cancer patients.

Maftools: efficient and comprehensive analysis of somatic variants in cancer.

Numerous large-scale genomic studies of matched tumor-normal samples have established the somatic landscapes of most cancer types. However, the downstream analysis of data from somatic mutations entails a number of computational and statistical approaches, requiring usage of independent software and numerous tools. Here, we describe an R Bioconductor package, Maftools, which offers a multitude of analysis and visualization modules that are commonly used in cancer genomic studies, including driver gene identification, pathway, signature, enrichment, and association analyses. Maftools only requires somatic variants in Mutation Annotation Format (MAF) and is independent of larger alignment files. With the implementation of well-established statistical and computational methods, Maftools facilitates data-driven research and comparative analysis to discover novel results from publicly available data sets. In the present study, using three of the well-annotated cohorts from The Cancer Genome Atlas (TCGA), we describe the application of Maftools to reproduce known results. More importantly, we show that Maftools can also be used to uncover novel findings through integrative analysis.

Integrative Analysis Defines Distinct Prognostic Subgroups of Intrahepatic Cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (iCCA) is the second most common primary liver cancer. It is defined by cholangiocytic differentiation and has poor prognosis. Recently, epigenetic processes have been shown to play an important role in cholangiocarcinogenesis. We performed an integrative analysis on 52 iCCAs using both genetic and epigenetic data with a specific focus on DNA methylation components. We found recurrent isocitrate dehydrogenase 1 (IDH1) and IDH2 (28%) gene mutations, recurrent arm-length copy number alterations (CNAs), and focal alterations such as deletion of 3p21 or amplification of 12q15, which affect BRCA1 Associated Protein 1, polybromo 1, and mouse double minute 2 homolog. DNA methylome analysis revealed excessive hypermethylation of iCCA, affecting primarily the bivalent genomic regions marked with both active and repressive histone modifications. Integrative clustering of genetic and epigenetic data identified four iCCA subgroups with prognostic relevance further designated as IDH, high (H), medium (M), and low (L) alteration groups. The IDH group consisted of all samples with IDH1 or IDH2 mutations and showed, together with the H group, a highly disrupted genome, characterized by frequent deletions of chromosome arms 3p and 6q. Both groups showed excessive hypermethylation with distinct patterns. The M group showed intermediate characteristics regarding both genetic and epigenetic marks, whereas the L group exhibited few methylation changes and mutations and a lack of CNAs. Methylation-based latent component analysis of cell-type composition identified differences among these four groups. Prognosis of the H and M groups was significantly worse than that of the L group. Conclusion: Using an integrative genomic and epigenomic analysis approach, we identified four major iCCA subgroups with widespread genomic and epigenomic differences and prognostic implications. Furthermore, our data suggest differences in the cell-of-origin of the iCCA subtypes.

Intratumor heterogeneity in epigenetic patterns.

Analogous to life on earth, tumor cells evolve through space and time and adapt to different micro-environmental conditions. As a result, tumors are composed of millions of genetically diversified cells at the time of diagnosis. Profiling these variants contributes to understanding tumors' clonal origins and might help to better understand response to therapy. However, even genetically homogenous cell populations show remarkable diversity in their response to different environmental stimuli, suggesting that genetic heterogeneity does not explain the full spectrum of tumor plasticity. Understanding epigenetic diversity across cancer cells provides important additional information about the functional state of subclones and therefore allows better understanding of tumor evolution and resistance to current therapies.

Histone methyltransferase SETDB1 contributes to melanoma tumorigenesis and serves as a new potential therapeutic target.

Alterations in histone modifications play a crucial role in the progression of various types of cancer. The histone methyltransferase SETDB1 catalyzes the addition of methyl groups to histone H3 at lysine 9. Here, we describe how overexpression of SETDB1 contributes to melanoma tumorigenesis. SETDB1 is highly amplified in melanoma cells and in the patient tumors. Increased expression of SETDB1, which correlates with SETDB1 amplification, is associated with a more aggressive phenotype in in vitro and in vivo studies. Mechanistically, SETDB1 implements its effects via regulation of thrombospondin 1, and the SET-domain of SETDB1 is essential for the maintenance of its tumorigenic activity. Inhibition of SETDB1 reduces cell growth in melanomas resistant to targeted treatments. Our results indicate that SETDB1 is a major driver of melanoma development and may serve as a potential future target for the treatment of this disease.

Enrichment analysis with EpiAnnotator.

Motivation: Deciphering relevant biological insights from epigenomic data can be a challenging task. One commonly used approach is to perform enrichment analysis. However, finding, downloading and using the publicly available functional annotations require time, programming skills and IT infrastructure. Here we describe the online tool EpiAnnotator for performing enrichment analyses on epigenomic data in a fast and user-friendly way. Results: EpiAnnotator is an R Package accompanied by a web interface. It contains regularly updated annotations from 4 public databases: Blueprint, RoadMap, GENCODE and the UCSC Genome Browser. Annotations are hosted locally or in a server environment and automatically updated by scripts of our own design. Thousands of tracks are available, reflecting data on a variety of tissues, cell types and cell lines from the human and mouse genomes. Users need to upload sets of selected and background regions. Results are displayed in customizable and easily interpretable figures. Availability and implementation: The R package and Shiny app are open source and available under the GPL v3 license. EpiAnnotator's web interface is accessible at http://computational-epigenomics.com/en/epiannotator. Contact: epiannotator@computational-epigenomics.com.

Circulating cytokines and small molecules follow distinct expression patterns in acute myeloid leukaemia.

BACKGROUND: Although aberrant expression of cytokines and small molecules (analytes) is well documented in acute myeloid leukaemia (AML), their co-expression patterns are not yet identified. In addition, plasma baselines for some analytes that are biomarkers for other cancers have not been previously reported in AML. METHODS: We used multiplex array technology to simultaneously detect and quantify 32 plasma analyte (22 reported analytes and 10 novel analytes) levels in 38 patients. RESULTS: In our study, 16 analytes are found to be significantly deregulated (13 higher, 3 lower, Mann-Whitney U-test, P-value <0.005), where 5 of them have never been reported before in AML. We predicted a seven-analyte-containing multiplex panel for diagnosis of AML and, among them, MIF could be a possible therapeutic target. In addition, we observed that circulating analytes show five co-expression signatures. CONCLUSIONS: Circulating analyte expression in AML significantly differs from normal, and follow distinct expression patterns.

Comprehensive analysis of DNA methylation data with RnBeads.

RnBeads is a software tool for large-scale analysis and interpretation of DNA methylation data, providing a user-friendly analysis workflow that yields detailed hypertext reports (http://rnbeads.mpi-inf.mpg.de/). Supported assays include whole-genome bisulfite sequencing, reduced representation bisulfite sequencing, Infinium microarrays and any other protocol that produces high-resolution DNA methylation data. Notable applications of RnBeads include the analysis of epigenome-wide association studies and epigenetic biomarker discovery in cancer cohorts.

Epigenetic profiling to classify cancer of unknown primary: a multicentre, retrospective analysis.

BACKGROUND: Cancer of unknown primary ranks in the top ten cancer presentations and has an extremely poor prognosis. Identification of the primary tumour and development of a tailored site-specific therapy could improve the survival of these patients. We examined the feasability of using DNA methylation profiles to determine the occult original cancer in cases of cancer of unknown primary. METHODS: We established a classifier of cancer type based on the microarray DNA methylation signatures (EPICUP) in a training set of 2790 tumour samples of known origin representing 38 tumour types and including 85 metastases. To validate the classifier, we used an independent set of 7691 known tumour samples from the same tumour types that included 534 metastases. We applied the developed diagnostic test to predict the tumour type of 216 well-characterised cases of cancer of unknown primary. We validated the accuracy of the predictions from the EPICUP assay using autopsy examination, follow-up for subsequent clinical detection of the primary sites months after the initial presentation, light microscopy, and comprehensive immunohistochemistry profiling. FINDINGS: The tumour type classifier based on the DNA methylation profiles showed a 99·6% specificity (95% CI 99·5-99·7), 97·7% sensitivity (96·1-99·2), 88·6% positive predictive value (85·8-91·3), and 99·9% negative predictive value (99·9-100·0) in the validation set of 7691 tumours. DNA methylation profiling predicted a primary cancer of origin in 188 (87%) of 216 patients with cancer with unknown primary. Patients with EPICUP diagnoses who received a tumour type-specific therapy showed improved overall survival compared with that in patients who received empiric therapy (hazard ratio [HR] 3·24, p=0·0051 [95% CI 1·42-7·38]; log-rank p=0·0029). INTERPRETATION: We show that the development of a DNA methylation based assay can significantly improve diagnoses of cancer of unknown primary and guide more precise therapies associated with better outcomes. Epigenetic profiling could be a useful approach to unmask the original primary tumour site of cancer of unknown primary cases and a step towards the improvement of the clinical management of these patients. FUNDING: European Research Council (ERC), Cellex Foundation, the Institute of Health Carlos III (ISCIII), Cancer Australia, Victorian Cancer Agency, Samuel Waxman Cancer Research Foundation, the Health and Science Departments of the Generalitat de Catalunya, and Ferrer.

A DNA methylation fingerprint of 1628 human samples.

Most of the studies characterizing DNA methylation patterns have been restricted to particular genomic loci in a limited number of human samples and pathological conditions. Herein, we present a compromise between an extremely comprehensive study of a human sample population with an intermediate level of resolution of CpGs at the genomic level. We obtained a DNA methylation fingerprint of 1628 human samples in which we interrogated 1505 CpG sites. The DNA methylation patterns revealed show this epigenetic mark to be critical in tissue-type definition and stemness, particularly around transcription start sites that are not within a CpG island. For disease, the generated DNA methylation fingerprints show that, during tumorigenesis, human cancer cells underwent a progressive gain of promoter CpG-island hypermethylation and a loss of CpG methylation in non-CpG-island promoters. Although transformed cells are those in which DNA methylation disruption is more obvious, we observed that other common human diseases, such as neurological and autoimmune disorders, had their own distinct DNA methylation profiles. Most importantly, we provide proof of principle that the DNA methylation fingerprints obtained might be useful for translational purposes by showing that we are able to identify the tumor type origin of cancers of unknown primary origin (CUPs). Thus, the DNA methylation patterns identified across the largest spectrum of samples, tissues, and diseases reported to date constitute a baseline for developing higher-resolution DNA methylation maps and provide important clues concerning the contribution of CpG methylation to tissue identity and its changes in the most prevalent human diseases.

Pediatric T-cell lymphoblastic leukemia evolves into relapse by clonal selection, acquisition of mutations and promoter hypomethylation.

Relapsed precursor T-cell acute lymphoblastic leukemia is characterized by resistance against chemotherapy and is frequently fatal. We aimed at understanding the molecular mechanisms resulting in relapse of T-cell acute lymphoblastic leukemia and analyzed 13 patients at first diagnosis, remission and relapse by whole exome sequencing, targeted ultra-deep sequencing, multiplex ligation dependent probe amplification and DNA methylation array. Compared to primary T-cell acute lymphoblastic leukemia, in relapse the number of single nucleotide variants and small insertions and deletions approximately doubled from 11.5 to 26. Targeted ultra-deep sequencing sensitively detected subclones that were selected for in relapse. The mutational pattern defined two types of relapses. While both are characterized by selection of subclones and acquisition of novel mutations, 'type 1' relapse derives from the primary leukemia whereas 'type 2' relapse originates from a common pre-leukemic ancestor. Relapse-specific changes included activation of the nucleotidase NT5C2 resulting in resistance to chemotherapy and mutations of epigenetic modulators, exemplified by SUZ12, WHSC1 and SMARCA4. While mutations present in primary leukemia and in relapse were enriched for known drivers of leukemia, relapse-specific changes revealed an association with general cancer-promoting mechanisms. This study thus identifies mechanisms that drive progression of pediatric T-cell acute lymphoblastic leukemia to relapse and may explain the characteristic treatment resistance of this condition.

A promoter DNA demethylation landscape of human hematopoietic differentiation.

Global mechanisms defining the gene expression programs specific for hematopoiesis are still not fully understood. Here, we show that promoter DNA demethylation is associated with the activation of hematopoietic-specific genes. Using genome-wide promoter methylation arrays, we identified 694 hematopoietic-specific genes repressed by promoter DNA methylation in human embryonic stem cells and whose loss of methylation in hematopoietic can be associated with gene expression. The association between promoter methylation and gene expression was studied for many hematopoietic-specific genes including CD45, CD34, CD28, CD19, the T cell receptor (TCR), the MHC class II gene HLA-DR, perforin 1 and the phosphoinositide 3-kinase (PI3K) and results indicated that DNA demethylation was not always sufficient for gene activation. Promoter demethylation occurred either early during embryonic development or later on during hematopoietic differentiation. Analysis of the genome-wide promoter methylation status of induced pluripotent stem cells (iPSCs) generated from somatic CD34(+) HSPCs and differentiated derivatives from CD34(+) HSPCs confirmed the role of DNA methylation in regulating the expression of genes of the hemato-immune system, and indicated that promoter methylation of these genes may be associated to stemness. Together, these data suggest that promoter DNA demethylation might play a role in the tissue/cell-specific genome-wide gene regulation within the hematopoietic compartment.

Disease-related blood-based differential methylation in cystic fibrosis and its representation in lung cancer revealed a regulatory locus in PKP3 in lung epithelial cells.

Cystic fibrosis (CF) is a monogenic disease, characterized by massive chronic lung inflammation. The observed variability in clinical phenotypes in monozygotic CF twins is likely associated with the extent of inflammation. This study sought to investigate inflammation-related aberrant DNA methylation in CF twins and to determine to what extent acquired methylation changes may be associated with lung cancer.Blood-based genome-wide DNA methylation analysis was performed to compare the DNA methylomes of monozygotic twins, from the European CF Twin and Sibling Study with various degrees of disease severity. Putatively inflammation-related and differentially methylated positions were selected from a large lung cancer case-control study and investigated in blood by targeted bisulphite next-generation-sequencing. An inflammation-related locus located in the Plakophilin-3 (PKP3) gene was functionally analysed regarding promoter and enhancer activity in presence and absence of methylation using luciferase reporter assays.We confirmed in a unique cohort that monozygotic twins, even if clinically discordant, have only minor differences in global DNA methylation patterns and blood cell composition. Further, we determined the most differentially methylated positions, a high proportion of which are blood cell-type-specific, whereas others may be acquired and thus have potential relevance in the context of inflammation as lung cancer risk factors. We identified a sequence in the gene body of PKP3 which is hypermethylated in blood from CF twins with severe phenotype and highly variably methylated in lung cancer patients and controls, independent of known clinical parameters, and showed that this region exhibits methylation-dependent promoter activity in lung epithelial cells.

Pediatric T-ALL type-1 and type-2 relapses develop along distinct pathways of clonal evolution.

The mechanisms underlying T-ALL relapse remain essentially unknown. Multilevel-omics in 38 matched pairs of initial and relapsed T-ALL revealed 18 (47%) type-1 (defined by being derived from the major ancestral clone) and 20 (53%) type-2 relapses (derived from a minor ancestral clone). In both types of relapse, we observed known and novel drivers of multidrug resistance including MDR1 and MVP, NT5C2 and JAK-STAT activators. Patients with type-1 relapses were specifically characterized by IL7R upregulation. In remarkable contrast, type-2 relapses demonstrated (1) enrichment of constitutional cancer predisposition gene mutations, (2) divergent genetic and epigenetic remodeling, and (3) enrichment of somatic hypermutator phenotypes, related to BLM, BUB1B/PMS2 and TP53 mutations. T-ALLs that later progressed to type-2 relapses exhibited a complex subclonal architecture, unexpectedly, already at the time of initial diagnosis. Deconvolution analysis of ATAC-Seq profiles showed that T-ALLs later developing into type-1 relapses resembled a predominant immature thymic T-cell population, whereas T-ALLs developing into type-2 relapses resembled a mixture of normal T-cell precursors. In sum, our analyses revealed fundamentally different mechanisms driving either type-1 or type-2 T-ALL relapse and indicate that differential capacities of disease evolution are already inherent to the molecular setup of the initial leukemia.

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.

Immunotyping of clinically divergent p.Phe508del homozygous monozygous cystic fibrosis twins.

Blood of the three clinically most concordant and most discordant p.Phe508del homozygous monozygous twin pairs of the European Cystic Fibrosis Twin and Sibling Study was examined in two postzygotic attributes that generate diversity between monozygous twins, i.e. the repertoire of the CDR3 region of the T-cell receptor ß chains and the DNA methylation at 450,000 genomic CpG sites. Methylation patterns in peripheral blood of twins changed at selected cell-type-independent positions and the immune cells of the twins showed individual profiles of the T cell receptor repertoire reflecting the plasticity of the immune system of genetically identical humans with cystic fibrosis to cope with the environment.

Methylome-based cell-of-origin modeling (Methyl-COOM) identifies aberrant expression of immune regulatory molecules in CLL.

BACKGROUND: In cancer, normal epigenetic patterns are disturbed and contribute to gene expression changes, disease onset, and progression. The cancer epigenome is composed of the epigenetic patterns present in the tumor-initiating cell at the time of transformation, and the tumor-specific epigenetic alterations that are acquired during tumor initiation and progression. The precise dissection of these two components of the tumor epigenome will facilitate a better understanding of the biological mechanisms underlying malignant transformation. Chronic lymphocytic leukemia (CLL) originates from differentiating B cells, which undergo extensive epigenetic programming. This poses the challenge to precisely determine the epigenomic ground state of the cell-of-origin in order to identify CLL-specific epigenetic aberrations. METHODS: We developed a linear regression model, methylome-based cell-of-origin modeling (Methyl-COOM), to map the cell-of-origin for individual CLL patients based on the continuum of epigenomic changes during normal B cell differentiation. RESULTS: Methyl-COOM accurately maps the cell-of-origin of CLL and identifies CLL-specific aberrant DNA methylation events that are not confounded by physiologic epigenetic B cell programming. Furthermore, Methyl-COOM unmasks abnormal action of transcription factors, altered super-enhancer activities, and aberrant transcript expression in CLL. Among the aberrantly regulated transcripts were many genes that have previously been implicated in T cell biology. Flow cytometry analysis of these markers confirmed their aberrant expression on malignant B cells at the protein level. CONCLUSIONS: Methyl-COOM analysis of CLL identified disease-specific aberrant gene regulation. The aberrantly expressed genes identified in this study might play a role in immune-evasion in CLL and might serve as novel targets for immunotherapy approaches. In summary, we propose a novel framework for in silico modeling of reference DNA methylomes and for the identification of cancer-specific epigenetic changes, a concept that can be broadly applied to other human malignancies.