MeDeCom: discovery and quantification of latent components of heterogeneous methylomes

MetDecode: methylation-based deconvolution of cell-free DNA for noninvasive multi-cancer typing

MOTIVATION

Circulating-cell free DNA (cfDNA) is widely explored as a noninvasive biomarker for cancer screening and diagnosis. The ability to decode the cells of origin in cfDNA would provide biological insights into pathophysiological mechanisms, aiding in cancer characterization and directing clinical management and follow-up.

RESULTS

We developed a DNA methylation signature-based deconvolution algorithm, MetDecode, for cancer tissue origin identification. We built a reference atlas exploiting de novo and published whole-genome methylation sequencing data for colorectal, breast, ovarian, and cervical cancer, and blood-cell-derived entities. MetDecode models the contributors absent in the atlas with methylation patterns learnt on-the-fly from the input cfDNA methylation profiles. In addition, our model accounts for the coverage of each marker region to alleviate potential sources of noise. In-silico experiments showed a limit of detection down to 2.88% of tumor tissue contribution in cfDNA. MetDecode produced Pearson correlation coefficients above 0.95 and outperformed other methods in simulations (P < 0.001; T-test; one-sided). In plasma cfDNA profiles from cancer patients, MetDecode assigned the correct tissue-of-origin in 84.2% of cases. In conclusion, MetDecode can unravel alterations in the cfDNA pool components by accurately estimating the contribution of multiple tissues, while supplied with an imperfect reference atlas.

AVAILABILITY AND IMPLEMENTATION

MetDecode is available at https://github.com/JorisVermeeschLab/MetDecode.

Computational deconvolution of DNA methylation data from mixed DNA samples

In this review, we provide a comprehensive overview of the different computational tools that have been published for the deconvolution of bulk DNA methylation (DNAm) data. Here, deconvolution refers to the estimation of cell-type proportions that constitute a mixed sample. The paper reviews and compares 25 deconvolution methods (supervised, unsupervised or hybrid) developed between 2012 and 2023 and compares the strengths and limitations of each approach. Moreover, in this study, we describe the impact of the platform used for the generation of methylation data (including microarrays and sequencing), the applied data pre-processing steps and the used reference dataset on the deconvolution performance. Next to reference-based methods, we also examine methods that require only partial reference datasets or require no reference set at all. In this review, we provide guidelines for the use of specific methods dependent on the DNA methylation data type and data availability.

Identification of region-specific gene isoforms in the human brain using long-read transcriptome sequencing

In neurological and neuropsychiatric diseases, different brain regions are affected, and differences in gene expression patterns could potentially explain this mechanism. However, limited studies have precisely explored gene expression in different regions of the human brain. In this study, we performed long-read RNA sequencing on three different brain regions of the same individuals: the cerebellum, hypothalamus, and temporal cortex. Despite stringent filtering criteria excluding isoforms predicted to be artifacts, over half of the isoforms expressed in multiple samples across multiple regions were found to be unregistered in the GENCODE reference. We then especially focused on genes with different major isoforms in each brain region, even with similar overall expression levels, and identified that many of such genes including GAS7 might have distinct roles in dendritic spine and neuronal formation in each region. We also found that DNA methylation might, in part, drive different isoform expressions in different regions. These findings highlight the significance of analyzing isoforms expressed in disease-relevant sites.

A Novel Type of IDH-wildtype Glioma Characterized by Gliomatosis Cerebri-like Growth Pattern, TERT Promoter Mutation, and Distinct Epigenetic Profile

Diffuse gliomas in adults encompass a heterogenous group of central nervous system neoplasms. In recent years, extensive (epi-)genomic profiling has identified several glioma subgroups characterized by distinct molecular characteristics, most importantly IDH1/2 and histone H3 mutations. A group of 16 diffuse gliomas classified as "adult-type diffuse high-grade glioma, IDH-wildtype, subtype F (HGG-F)" was identified by the DKFZ v12.5 Brain Tumor Classifier . Histopathologic characterization, exome sequencing, and review of clinical data was performed in all cases. Based on unsupervised t -distributed stochastic neighbor embedding and clustering analysis of genome-wide DNA methylation data, HGG-F shows distinct epigenetic profiles separate from established central nervous system tumors. Exome sequencing demonstrated frequent TERT promoter (12/15 cases), PIK3R1 (11/16), and TP53 mutations (5/16). Radiologic characteristics were reminiscent of gliomatosis cerebri in 9/14 cases (64%). Histopathologically, most cases were classified as diffuse gliomas (7/16, 44%) or were suspicious for the infiltration zone of a diffuse glioma (5/16, 31%). None of the cases demonstrated microvascular proliferation or necrosis. Outcome of 14 patients with follow-up data was better compared to IDH-wildtype glioblastomas with a median progression-free survival of 58 months and overall survival of 74 months (both P <0.0001). Our series represents a novel type of adult-type diffuse glioma with distinct molecular and clinical features. Importantly, we provide evidence that TERT promoter mutations in diffuse gliomas without further morphologic or molecular signs of high-grade glioma should be interpreted in the context of the clinicoradiologic presentation as well as epigenetic profile and may not be suitable as a standalone marker for glioblastoma, IDH-wildtype.

EpiMix is an integrative tool for epigenomic subtyping using DNA methylation

DNA methylation (DNAme) is a major epigenetic factor influencing gene expression with alterations leading to cancer and immunological and cardiovascular diseases. Recent technological advances have enabled genome-wide profiling of DNAme in large human cohorts. There is a need for analytical methods that can more sensitively detect differential methylation profiles present in subsets of individuals from these heterogeneous, population-level datasets. We developed an end-to-end analytical framework named "EpiMix" for population-level analysis of DNAme and gene expression. Compared with existing methods, EpiMix showed higher sensitivity in detecting abnormal DNAme that was present in only small patient subsets. We extended the model-based analyses of EpiMix to cis-regulatory elements within protein-coding genes, distal enhancers, and genes encoding microRNAs and long non-coding RNAs (lncRNAs). Using cell-type-specific data from two separate studies, we discover epigenetic mechanisms underlying childhood food allergy and survival-associated, methylation-driven ncRNAs in non-small cell lung cancer.

EpiMix: an integrative tool for epigenomic subtyping using DNA methylation

DNA methylation (DNAme) is a major epigenetic factor influencing gene expression with alterations leading to cancer, immunological, and cardiovascular diseases. Recent technological advances enable genome-wide quantification of DNAme in large human cohorts. So far, existing methods have not been evaluated to identify differential DNAme present in large and heterogeneous patient cohorts. We developed an end-to-end analytical framework named "EpiMix" for population-level analysis of DNAme and gene expression. Compared to existing methods, EpiMix showed higher sensitivity in detecting abnormal DNAme that was present in only small patient subsets. We extended the model-based analyses of EpiMix to cis-regulatory elements within protein-coding genes, distal enhancers, and genes encoding microRNAs and lncRNAs. Using cell-type specific data from two separate studies, we discovered novel epigenetic mechanisms underlying childhood food allergy and survival-associated, methylation-driven non-coding RNAs in non-small cell lung cancer.

Gliomatosis cerebri (GC) growth pattern: A single-center analysis of clinical, histological, and molecular characteristics of GC and non-GC glioblastoma

Background

The biological understanding of glioblastoma (GB) with gliomatosis cerebri (GC) pattern is poor due to the absence of GC-specific studies. Here, we aimed to identify molecular or clinical parameters that drive GC growth.

Methods

From our methylome database of IDH (isocitrate dehydrogenase)-wildtype GB, we identified 158 non-GC and 65 GC cases. GC cases were subdivided into diffuse-infiltrative (subtype 1), multifocal (subtype 2), or tumors with 1 solid mass (subtype 3). We compared clinical, histological, and molecular parameters and conducted a reference-free tumor deconvolution of DNA methylation data based on latent methylation components (LMC).

Results

GC subtype 1 less frequently showed contrast-enhancing tumors, and more frequently lacked morphological GB criteria despite displaying GB DNA methylation profile. However, the tumor deconvolution did not deliver a specific LMC cluster for either of the GC subtypes. Employing the reference-based analysis MethylCIBERSORT, we did not identify significant differences in tumor cell composition. The majority of both GC and non-GC patients received radiochemotherapy as first-line treatment, but there was a major imbalance for resection. The entire GC cohort had significantly shorter overall survival (OS) and time to treatment failure (TTF) than the non-GC cohort. However, when filtering for cases in which only stereotactic biopsy was performed, the comparison of OS and TTF lost statistical significance.

Conclusions

Our study offers clinically relevant information by demonstrating a similar outcome for GB with GC growth pattern in the surgically matched analysis. The limited number of cases in the GC subgroups encourages the validation of our DNA methylation analysis in larger cohorts.

Ageing-associated small RNA cargo of extracellular vesicles

Previous work on murine models and humans demonstrated global as well as tissue-specific molecular ageing trajectories of RNAs. Extracellular vesicles (EVs) are membrane vesicles mediating the horizontal transfer of genetic information between different tissues. We sequenced small regulatory RNAs (sncRNAs) in two mouse plasma fractions at five time points across the lifespan from 2-18 months: (1) sncRNAs that are free-circulating (fc-RNA) and (2) sncRNAs bound outside or inside EVs (EV-RNA). Different sncRNA classes exhibit unique ageing patterns that vary between the fcRNA and EV-RNA fractions. While tRNAs showed the highest correlation with ageing in both fractions, rRNAs exhibited inverse correlation trajectories between the EV- and fc-fractions. For miRNAs, the EV-RNA fraction was exceptionally strongly associated with ageing, especially the miR-29 family in adipose tissues. Sequencing of sncRNAs and coding genes in fat tissue of an independent cohort of aged mice up to 27 months highlighted the pivotal role of miR-29a-3p and miR-29b-3p in ageing-related gene regulation that we validated in a third cohort by RT-qPCR.

Multi-omics research strategies in ischemic stroke: A multidimensional perspective

Ischemic stroke (IS) is a multifactorial and heterogeneous neurological disorder with high rate of death and long-term impairment. Despite years of studies, there are still no stroke biomarkers for clinical practice, and the molecular mechanisms of stroke remain largely unclear. The high-throughput omics approach provides new avenues for discovering biomarkers of IS and explaining its pathological mechanisms. However, single-omics approaches only provide a limited understanding of the biological pathways of diseases. The integration of multiple omics data means the simultaneous analysis of thousands of genes, RNAs, proteins and metabolites, revealing networks of interactions between multiple molecular levels. Integrated analysis of multi-omics approaches will provide helpful insights into stroke pathogenesis, therapeutic target identification and biomarker discovery. Here, we consider advances in genomics, transcriptomics, proteomics and metabolomics and outline their use in discovering the biomarkers and pathological mechanisms of IS. We then delineate strategies for achieving integration at the multi-omics level and discuss how integrative omics and systems biology can contribute to our understanding and management of IS.

Genetic and methylation profiles distinguish benign, malignant and spitzoid melanocytic tumors

Accurate classification of melanocytic tumors is important for prognostic evaluation, treatment and follow-up protocols of patients. The majority of melanocytic proliferations can be classified solely based on clinical and pathological criteria, however in select cases a definitive diagnostic assessment remains challenging and additional diagnostic biomarkers would be advantageous. We analyzed melanomas, nevi, Spitz nevi and atypical spitzoid tumors using parallel sequencing (exons of 611 genes and 507 gene translocation analysis) and methylation arrays (850k Illumina EPIC). By combining detailed genetic and epigenetic analysis with reference-based and reference-free DNA methylome deconvolution we compared Spitz nevi to nevi and melanoma and assessed the potential for these methods in classifying challenging spitzoid tumors. Results were correlated with clinical and histologic features. Spitz nevi were found to cluster independently of nevi and melanoma and demonstrated a different mutation profile. Multiple copy number alterations and TERT promoter mutations were identified only in melanomas. Genome-wide methylation in Spitz nevi was comparable to benign nevi while the Leukocytes UnMethylation for Purity (LUMP) algorithm in Spitz nevi was comparable to melanoma. Histologically difficult to classify Spitz tumor cases were assessed which, based on methylation arrays, clustered between Spitz nevi and melanoma and in terms of genetic profile or copy number variations demonstrated worrisome features suggesting a malignant neoplasm. Comprehensive sequencing and methylation analysis verify Spitz nevi as an independent melanocytic entity distinct from both nevi and melanoma. Combined genetic and methylation assays can offer additional insights in diagnosing difficult to classify Spitzoid tumors.

Deconvolution of tumor composition using partially available DNA methylation data

Background

Deciphering proportions of constitutional cell types in tumor tissues is a crucial step for the analysis of tumor heterogeneity and the prediction of response to immunotherapy. In the process of measuring cell population proportions, traditional experimental methods have been greatly hampered by the cost and extensive dropout events. At present, the public availability of large amounts of DNA methylation data makes it possible to use computational methods to predict proportions.

Results

In this paper, we proposed PRMeth, a method to deconvolve tumor mixtures using partially available DNA methylation data. By adopting an iteratively optimized non-negative matrix factorization framework, PRMeth took DNA methylation profiles of a portion of the cell types in the tissue mixtures (including blood and solid tumors) as input to estimate the proportions of all cell types as well as the methylation profiles of unknown cell types simultaneously. We compared PRMeth with five different methods through three benchmark datasets and the results show that PRMeth could infer the proportions of all cell types and recover the methylation profiles of unknown cell types effectively. Then, applying PRMeth to four types of tumors from The Cancer Genome Atlas (TCGA) database, we found that the immune cell proportions estimated by PRMeth were largely consistent with previous studies and met biological significance.

Conclusions

Our method can circumvent the difficulty of obtaining complete DNA methylation reference data and obtain satisfactory deconvolution accuracy, which will be conducive to exploring the new directions of cancer immunotherapy. PRMeth is implemented in R and is freely available from GitHub (https://github.com/hedingqin/PRMeth).

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04893-7.

Genetic and Methylation Analysis of CTNNB1 in Benign and Malignant Melanocytic Lesions

Simple Summary

Recurrent CTNNB1 exon 3 mutations have been recognized in the distinct group of melanocytic tumors showing deep penetrating nevus-like morphology and in 1–2% of advanced melanoma. We performed a detailed genetic analysis of difficult-to-classify nevi and melanomas with CTNNB1 mutations and found that benign tumors (nevi) show characteristic morphological, genetic and epigenetic traits, which distinguish them from other nevi and melanoma. Malignant CTNNB1-mutant tumors (melanoma) demonstrated a different genetic profile, grouping clearly with other non-CTNNB1 melanomas in methylation assays. To further evaluate the role of CTNNB1 mutations in melanoma, we assessed a large cohort of clinically sequenced melanomas, identifying 38 tumors with CTNNB1 exon 3 mutations, including recurrent S45 (n = 13, 34%), G34 (n = 5, 13%), and S27 (n = 5, 13%) mutations. Locations and histological subtype of CTNNB1-mutated melanoma varied; none were reported as showing deep penetrating nevus-like morphology. The most frequent concurrent activating mutations were BRAF V600 (55%) and NRAS Q61 (34%).

Abstract

Melanocytic neoplasms have been genetically characterized in detail during the last decade. Recurrent CTNNB1 exon 3 mutations have been recognized in the distinct group of melanocytic tumors showing deep penetrating nevus-like morphology. In addition, they have been identified in 1–2% of advanced melanoma. Performing a detailed genetic analysis of difficult-to-classify nevi and melanomas with CTNNB1 mutations, we found that benign tumors (nevi) show characteristic morphological, genetic and epigenetic traits, which distinguish them from other nevi and melanoma. Malignant CTNNB1-mutant tumors (melanomas) demonstrated a different genetic profile, instead grouping clearly with other non-CTNNB1 melanomas in methylation assays. To further evaluate the role of CTNNB1 mutations in melanoma, we assessed a large cohort of clinically sequenced melanomas, identifying 38 tumors with CTNNB1 exon 3 mutations, including recurrent S45 (n = 13, 34%), G34 (n = 5, 13%), and S27 (n = 5, 13%) mutations. Locations and histological subtype of CTNNB1-mutated melanoma varied; none were reported as showing deep penetrating nevus-like morphology. The most frequent concurrent activating mutations were BRAF V600 (n = 21, 55%) and NRAS Q61 (n = 13, 34%). In our cohort, four of seven (58%) and one of nine (11%) patients treated with targeted therapy (BRAF and MEK Inhibitors) or immune-checkpoint therapy, respectively, showed disease control (partial response or stable disease). In summary, CTNNB1 mutations are associated with a unique melanocytic tumor type in benign tumors (nevi), which can be applied in a diagnostic setting. In advanced disease, no clear characteristics distinguishing CTNNB1-mutant from other melanomas were observed; however, studies of larger, optimally prospective, cohorts are warranted.

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.

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.

DNA methylation profiling of central nervous system hemangioblastomas identifies two distinct subgroups

Hemangioblastomas (HBs) of the central nervous system are highly vascular neoplasms that occur sporadically or as a manifestation of von Hippel–Lindau (VHL) disease. Despite their benign nature, HBs are clinically heterogeneous and can be associated with significant morbidity due to mass effects of peritumoral cysts or tumor progression. Underlying molecular factors involved in HB tumor biology remain elusive. We investigated genome‐wide DNA methylation profiles and clinical and histopathological features in a series of 47 HBs from 42 patients, including 28 individuals with VHL disease. Thirty tumors occurred in the cerebellum, 8 in the brainstem and 8 HBs were of spinal location, while 1 HB was located in the cerebrum. Histologically, 12 HBs (26%) belonged to the cellular subtype and exclusively occurred in the cerebellum, whereas 35 HBs were reticular (74%). Unsupervised clustering and dimensionality reduction of DNA methylation profiles revealed two distinct subgroups. Methylation cluster 1 comprised 30 HBs of mainly cerebellar location (29/30, 97%), whereas methylation cluster 2 contained 17 HBs predominantly located in non‐cerebellar compartments (16/17, 94%). The sum of chromosomal regions being affected by copy‐number alterations was significantly higher in methylation cluster 1 compared to cluster 2 (mean 262 vs. 109 Mb, p = 0.001). Of note, loss of chromosome 6 occurred in 9/30 tumors (30%) of methylation cluster 1 and was not observed in cluster 2 tumors (p = 0.01). No relevant methylation differences between sporadic and VHL‐related HBs or cystic and non‐cystic HBs could be detected. Deconvolution of the bulk DNA methylation profiles revealed four methylation components that were associated with the two methylation clusters suggesting cluster‐specific cell‐type compositions. In conclusion, methylation profiling of HBs reveals 2 distinct subgroups that mainly associate with anatomical location, cytogenetic profiles and differences in cell type composition, potentially reflecting different cells of origin.

A pan-tissue DNA methylation atlas enables in silico decomposition of human tissue methylomes at cell-type resolution

Bulk-tissue DNA methylomes represent an average over many different cell types, hampering our understanding of cell-type-specific contributions to disease development. As single-cell methylomics is not scalable to large cohorts of individuals, cost-effective computational solutions are needed, yet current methods are limited to tissues such as blood. Here we leverage the high-resolution nature of tissue-specific single-cell RNA-sequencing datasets to construct a DNA methylation atlas defined for 13 solid tissue types and 40 cell types. We comprehensively validate this atlas in independent bulk and single-nucleus DNA methylation datasets. We demonstrate that it correctly predicts the cell of origin of diverse cancer types and discovers new prognostic associations in olfactory neuroblastoma and stage 2 melanoma. In brain, the atlas predicts a neuronal origin for schizophrenia, with neuron-specific differential DNA methylation enriched for corresponding genome-wide association study risk loci. In summary, the DNA methylation atlas enables the decomposition of 13 different human tissue types at a high cellular resolution, paving the way for an improved interpretation of epigenetic data.

This resource presents an in silico generated DNA methylation atlas that can be used for cell-type deconvolution of human tissues.

Novel epigenetic network biomarkers for early detection of esophageal cancer

BACKGROUND

Early detection of esophageal cancer is critical to improve survival. Whilst studies have identified biomarkers, their interpretation and validity is often confounded by cell-type heterogeneity.

RESULTS

Here we applied systems-epigenomic and cell-type deconvolution algorithms to a discovery set encompassing RNA-Seq and DNA methylation data from esophageal adenocarcinoma (EAC) patients and matched normal-adjacent tissue, in order to identify robust biomarkers, free from the confounding effect posed by cell-type heterogeneity. We identify 12 gene-modules that are epigenetically deregulated in EAC, and are able to validate all 12 modules in 4 independent EAC cohorts. We demonstrate that the epigenetic deregulation is present in the epithelial compartment of EAC-tissue. Using single-cell RNA-Seq data we show that one of these modules, a proto-cadherin module centered around CTNND2, is inactivated in Barrett's Esophagus, a precursor lesion to EAC. By measuring DNA methylation in saliva from EAC cases and controls, we identify a chemokine module centered around CCL20, whose methylation patterns in saliva correlate with EAC status.

CONCLUSIONS

Given our observations that a CCL20 chemokine network is overactivated in EAC tissue and saliva from EAC patients, and that in independent studies CCL20 has been found to be overactivated in EAC tissue infected with the bacterium F. nucleatum, a bacterium that normally inhabits the oral cavity, our results highlight the possibility of using DNAm measurements in saliva as a proxy for changes occurring in the esophageal epithelium. Both the CTNND2/CCL20 modules represent novel promising network biomarkers for EAC that merit further investigation.

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.

Methylation and mutation of the inhibin‑α gene in human melanoma cells and regulation of PTEN expression and AKT/PI3K signaling by a demethylating agent

Inhibin suppresses the pituitary secretion of follicle‑stimulating hormone and has been reported to act as a tumor suppressor gene in the gonad in mice. Epigenetic modifications, mutations, changes in the loss of heterozygosity (LOH) of the inhibin‑α gene and regulation of gene expression in response to a demethylating agent [5‑aza‑2'‑deoxycytidine (5‑Aza‑dC)] in human melanoma cells were assessed. In addition, the association between a mutation in the 5'‑untranslated region (5'‑UTR) of the inhibin‑α subunit and the expression of phosphatidylinositol 3,4,5‑trisphosphate‑dependent Rac exchanger 2 (PREX2) and phosphatase and tensin homolog (PTEN) as well as AKT/PI3K signaling was determined. The methylation status of the CpG sites of the inhibin‑α promoter was analyzed by methylation‑specific PCR in bisulfite‑treated DNA. Cell viability was counted using the trypan blue assay, mRNA expression was examined via reverse transcription‑quantitative PCR, and protein expression was examined via western blot analysis. The inhibin‑α promoter was hypermethylated in G361, SK‑MEL‑3, SK‑MEL‑24 and SK‑MEL‑28 cells and moderately methylated in SK‑MEL‑5 cells. Inhibin‑α gene mutations were observed in the 5'‑UTR exon 1 of G361, SK‑MEL‑5, SK‑MEL‑24 and SK‑MEL‑28 cells as well as in exon 2 of SK‑MEL‑3 cells. Allelic imbalance, including LOH, in the inhibin‑α gene was detected in human melanoma cells. Treatment with 5‑Aza‑dC increased inhibin‑α mRNA and protein levels, inhibited cell proliferation, and delayed the doubling times of surviving melanoma cells. In 5‑Aza‑dC‑treated cells, PREX2 protein expression was slightly increased in G361 and SK‑MEL‑24 cells and decreased in SK‑MEL3, SK‑MEL‑5 and SK‑MEL‑28 cells. However, the protein expression of PTEN was decreased in melanoma cells. In addition, AKT and PI3K protein phosphorylation levels increased in all melanoma cells, except of G361 cells, demonstrating decreased PI3K protein phosphorylation. These data provided evidence that methylation, mutation and LOH are observed in the inhibin α‑subunit gene and gene locus in human melanoma cells. Furthermore, the demethylating agent reactivated inhibin‑α gene expression and regulated PREX2 expression. AKT/PI3K signaling increased as PTEN expression decreased. In addition, mutations in the tumor suppressor inhibin‑α, PTEN and p53 genes were not associated with transcriptional silencing, gene expression and cell growth as analyzed through experiments and literature reviews. These data demonstrated that methylation and mutations were associated with the inhibin‑α gene in human melanoma cells and indicated the regulation of PTEN expression and AKT/PI3K signaling by a demethylating agent.

A transdisciplinary approach to understand the epigenetic basis of race/ethnicity health disparities

Health disparities correspond to differences in disease burden and mortality among socially defined population groups. Such disparities may emerge according to race/ethnicity, socioeconomic status and a variety of other social contexts, and are documented for a wide range of diseases. Here, we provide a transdisciplinary perspective on the contribution of epigenetics to the understanding of health disparities, with a special emphasis on disparities across socially defined racial/ethnic groups. Scientists in the fields of biological anthropology, bioinformatics and molecular epidemiology provide a summary of theoretical, statistical and practical considerations for conducting epigenetic health disparities research, and provide examples of successful applications from cancer research using this approach.

DECONbench: a benchmarking platform dedicated to deconvolution methods for tumor heterogeneity quantification

BACKGROUND

Quantification of tumor heterogeneity is essential to better understand cancer progression and to adapt therapeutic treatments to patient specificities. Bioinformatic tools to assess the different cell populations from single-omic datasets as bulk transcriptome or methylome samples have been recently developed, including reference-based and reference-free methods. Improved methods using multi-omic datasets are yet to be developed in the future and the community would need systematic tools to perform a comparative evaluation of these algorithms on controlled data.

RESULTS

We present DECONbench, a standardized unbiased benchmarking resource, applied to the evaluation of computational methods quantifying cell-type heterogeneity in cancer. DECONbench includes gold standard simulated benchmark datasets, consisting of transcriptome and methylome profiles mimicking pancreatic adenocarcinoma molecular heterogeneity, and a set of baseline deconvolution methods (reference-free algorithms inferring cell-type proportions). DECONbench performs a systematic performance evaluation of each new methodological contribution and provides the possibility to publicly share source code and scoring.

CONCLUSION

DECONbench allows continuous submission of new methods in a user-friendly fashion, each novel contribution being automatically compared to the reference baseline methods, which enables crowdsourced benchmarking. DECONbench is designed to serve as a reference platform for the benchmarking of deconvolution methods in the evaluation of cancer heterogeneity. We believe it will contribute to leverage the benchmarking practices in the biomedical and life science communities. DECONbench is hosted on the open source Codalab competition platform. It is freely available at: https://competitions.codalab.org/competitions/27453 .

Impact of the methylation classifier and ancillary methods on CNS tumor diagnostics

BACKGROUND

Accurate CNS tumor diagnosis can be challenging, and methylation profiling can serve as an adjunct to classify diagnostically difficult cases.

METHODS

An integrated diagnostic approach was employed for a consecutive series of 1,258 surgical neuropathology samples obtained primarily in a consultation practice over 2-year period. DNA methylation profiling and classification using the DKFZ/Heidelberg CNS tumor classifier was performed, as well as unsupervised analyses of methylation data. Ancillary testing, where relevant, was performed.

RESULTS

Among the received cases in consultation, a high confidence methylation classifier score (>0.84) was reached in 66.4% of cases. The classifier impacted the diagnosis in 46.5% of these high-confidence classifier score cases, including a substantially new diagnosis in 26.9% cases. Among the 289 cases received with only a descriptive diagnosis, methylation was able to resolve approximately half (144, 49.8%) with high-confidence scores. Additional methods were able to resolve diagnostic uncertainty in 41.6% of the low-score cases. Tumor purity was significantly associated with classifier score (p = 1.15e-11). Deconvolution demonstrated that suspected GBMs matching as control/inflammatory brain tissue could be resolved into GBM methylation profiles, which provided a proof-of-concept approach to resolve tumor classification in the setting of low tumor purity.

CONCLUSIONS

This work assesses the impact of a methylation classifier and additional methods in a consultative practice by defining the proportions with concordant vs. change in diagnosis in a set of diagnostically challenging CNS tumors. We address approaches to low-confidence scores and confounding issues of low tumor purity.

DNA methylation-based prediction of response to immune checkpoint inhibition in metastatic melanoma

Background

Therapies based on targeting immune checkpoints have revolutionized the treatment of metastatic melanoma in recent years. Still, biomarkers predicting long-term therapy responses are lacking.

Methods

A novel approach of reference-free deconvolution of large-scale DNA methylation data enabled us to develop a machine learning classifier based on CpG sites, specific for latent methylation components (LMC), that allowed for patient allocation to prognostic clusters. DNA methylation data were processed using reference-free analyses (MeDeCom) and reference-based computational tumor deconvolution (MethylCIBERSORT, LUMP).

Results

We provide evidence that DNA methylation signatures of tumor tissue from cutaneous metastases are predictive for therapy response to immune checkpoint inhibition in patients with stage IV metastatic melanoma.

Conclusions

These results demonstrate that LMC-based segregation of large-scale DNA methylation data is a promising tool for classifier development and treatment response estimation in cancer patients under targeted immunotherapy.

Detecting cord blood cell type-specific epigenetic associations with gestational diabetes mellitus and early childhood growth

BACKGROUND

Epigenome-wide association studies (EWAS) have provided opportunities to understand the role of epigenetic mechanisms in development and pathophysiology of many chronic diseases. However, an important limitation of conventional EWAS is that profiles of epigenetic variability are often obtained in samples of mixed cell types. Here, we aim to assess whether changes in cord blood DNA methylation (DNAm) associated with gestational diabetes mellitus (GDM) exposure and early childhood growth markers occur in a cell type-specific manner.

RESULTS

We analyzed 275 cord blood samples collected at delivery from a prospective pre-birth cohort with genome-wide DNAm profiled by the Illumina MethylationEPIC array. We estimated proportions of seven common cell types in each sample using a cord blood-specific DNAm reference panel. Leveraging a recently developed approach named CellDMC, we performed cell type-specific EWAS to identify CpG loci significantly associated with GDM, or 3-year-old body mass index (BMI) z-score. A total of 1410 CpG loci displayed significant cell type-specific differences in methylation level between 23 GDM cases and 252 controls with a false discovery rate < 0.05. Gene Ontology enrichment analysis indicated that LDL transportation emerged from CpG specifically identified from B-cells DNAm analyses and the mitogen-activated protein kinase pathway emerged from CpG specifically identified from natural killer cells DNAm analyses. In addition, we identified four and six loci associated with 3-year-old BMI z-score that were specific to CD8+ T-cells and monocytes, respectively. By performing genome-wide permutation tests, we validated that most of our detected signals had low false positive rates.

CONCLUSION

Compared to conventional EWAS adjusting for the effects of cell type heterogeneity, the proposed approach based on cell type-specific EWAS could provide additional biologically meaningful associations between CpG methylation, prenatal maternal GDM or 3-year-old BMI. With careful validation, these findings may provide new insights into the pathogenesis, programming, and consequences of related childhood metabolic dysregulation. Therefore, we propose that cell type-specific analyses are worth cautious explorations.

Upper Aerodigestive Tract Squamous Cell Carcinomas Show Distinct Overall DNA Methylation Profiles and Different Molecular Mechanisms behind WNT Signaling Disruption

Upper aerodigestive tract (UADT) tumors present different biological behavior and prognosis, suggesting specific molecular mechanisms underlying their development. However, they are rarely considered as single entities (particularly head and neck subsites) and share the most common genetic alterations. Therefore, there is a need for a better understanding of the global DNA methylation differences among UADT tumors. We performed a genome-wide DNA methylation analysis of esophageal (ESCC), laryngeal (LSCC), oral (OSCC) and oropharyngeal (OPSCC) squamous cell carcinomas, and their non-tumor counterparts. The unsupervised analysis showed that non-tumor tissues present markedly distinct DNA methylation profiles, while tumors are highly heterogeneous. Hypomethylation was more frequent in LSCC and OPSCC, while ESCC and OSCC presented mostly hypermethylation, with the latter showing a CpG island overrepresentation. Differentially methylated regions affected genes in 127 signaling pathways, with only 3.1% of these being common among different tumor subsites, but with different genes affected. The WNT signaling pathway, known to be dysregulated in different epithelial tumors, is a frequent hit for DNA methylation and gene expression alterations in ESCC and OPSCC, but mostly for genetic alterations in LSCC and OSCC. UADT tumor subsites present differences in genome-wide methylation regarding their profile, intensity, genomic regions and signaling pathways affected.

Deconvolution of sarcoma methylomes reveals varying degrees of immune cell infiltrates with association to genomic aberrations

BACKGROUND

Soft-tissue sarcomas (STS) are a heterogeneous group of mesenchymal tumors for which response to immunotherapies is not well established. Therefore, it is important to risk-stratify and identify STS patients who will most likely benefit from these treatments.

RESULTS

To reveal shared and distinct methylation signatures present in STS, we performed unsupervised deconvolution of DNA methylation data from the TCGA sarcoma and an independent validation cohort. We showed that leiomyosarcoma can be subclassified into three distinct methylation groups. More importantly, we identified a component associated with tumor-infiltrating leukocytes, which suggests varying degrees of immune cell infiltration in STS subtypes and an association with prognosis. We further investigated the genomic alterations that may influence tumor infiltration by leukocytes including RB1 loss in undifferentiated pleomorphic sarcomas and ELK3 amplification in dedifferentiated liposarcomas.

CONCLUSIONS

In summary, we have leveraged unsupervised methylation-based deconvolution to characterize the immune compartment and molecularly stratify subtypes in STS, which may benefit precision medicine in the future.

Machine learning for deciphering cell heterogeneity and gene regulation

Epigenetics studies inheritable and reversible modifications of DNA that allow cells to control gene expression throughout their development and in response to environmental conditions. In computational epigenomics, machine learning is applied to study various epigenetic mechanisms genome wide. Its aim is to expand our understanding of cell differentiation, that is their specialization, in health and disease. Thus far, most efforts focus on understanding the functional encoding of the genome and on unraveling cell-type heterogeneity. Here, we provide an overview of state-of-the-art computational methods and their underlying statistical concepts, which range from matrix factorization and regularized linear regression to deep learning methods. We further show how the rise of single-cell technology leads to new computational challenges and creates opportunities to further our understanding of epigenetic regulation.

Epigenetic signatures in cancer: proper controls, current challenges and the potential for clinical translation

Epigenetic alterations are associated with normal biological processes such as aging or differentiation. Changes in global epigenetic signatures, together with genetic alterations, are driving events in several diseases including cancer. Comparative studies of cancer and healthy tissues found alterations in patterns of DNA methylation, histone posttranslational modifications, and changes in chromatin accessibility. Driven by sophisticated, next-generation sequencing-based technologies, recent studies discovered cancer epigenomes to be dominated by epigenetic patterns already present in the cell-of-origin, which transformed into a neoplastic cell. Tumor-specific epigenetic changes therefore need to be redefined and factors influencing epigenetic patterns need to be studied to unmask truly disease-specific alterations. The underlying mechanisms inducing cancer-associated epigenetic alterations are poorly understood. Studies of mutated epigenetic modifiers, enzymes that write, read, or edit epigenetic patterns, or mutated chromatin components, for example oncohistones, help to provide functional insights on how cancer epigenomes arise. In this review, we highlight the importance and define challenges of proper control tissues and cell populations to exploit cancer epigenomes. We summarize recent advances describing mechanisms leading to epigenetic changes in tumorigenesis and briefly discuss advances in investigating their translational potential.

Systems Biology Analysis for Ewing Sarcoma

Ewing sarcoma (EwS) is a highly aggressive pediatric bone cancer that is defined by a somatic fusion between the EWSR1 gene and an ETS family member, most frequently the FLI1 gene, leading to expression of a chimeric transcription factor EWSR1-FLI1. Otherwise, EwS is one of the most genetically stable cancers. The situation when the major cancer driver is well known looks like a unique opportunity for applying the systems biology approach in order to understand the EwS mechanisms as well as to uncover some general mechanistic principles of carcinogenesis. A number of studies have been performed revealing the direct and indirect effects of EWSR1-FLI1 on multiple aspects of cellular life. Nevertheless, the emerging picture of the oncogene action appears to be highly complex and systemic, with multiple reciprocal influences between the immediate consequences of the driver mutation and intracellular and intercellular molecular mechanisms, including regulation of transcription, epigenome, and tumoral microenvironment. In this chapter, we present an overview of existing molecular profiling resources available for EwS tumors and cell lines and provide an online comprehensive catalogue of publicly available omics and other datasets. We further highlight the systems biology studies of EwS, involving mathematical modeling of networks and integration of molecular data. We conclude that despite the seeming simplicity, a lot has yet to be understood on the systems-wide mechanisms connecting the driver mutation and the major cellular phenotypes of this pediatric cancer. Overall, this chapter can serve as a guide for a systems biology researcher to start working on EwS.

Deconvolution of heterogeneous tumor samples using partial reference signals

Deconvolution of heterogeneous bulk tumor samples into distinct cellular populations is an important yet challenging problem, particularly when only partial references are available. A common approach to dealing with this problem is to deconvolve the mixed signals using available references and leverage the remaining signal as a new cell component. However, as indicated in our simulation, such an approach tends to over-estimate the proportions of known cell types and fails to detect novel cell types. Here, we propose PREDE, a partial reference-based deconvolution method using an iterative non-negative matrix factorization algorithm. Our method is verified to be effective in estimating cell proportions and expression profiles of unknown cell types based on simulated datasets at a variety of parameter settings. Applying our method to TCGA tumor samples, we found that proportions of pure cancer cells better indicate different subtypes of tumor samples. We also detected several cell types for each cancer type whose proportions successfully predicted patient survival. Our method makes a significant contribution to deconvolution of heterogeneous tumor samples and could be widely applied to varieties of high throughput bulk data. PREDE is implemented in R and is freely available from GitHub (https://xiaoqizheng.github.io/PREDE).

Tumor tissues are mixtures of different cell types. Identification and quantification of constitutional cell types within tumor tissues are important tasks in cancer research. The problem can be readily solved using regression-based methods if reference signals are available. But in most clinical applications, only partial references are available, which significantly reduces the deconvolution accuracy of the existing regression-based methods. In this paper, we propose a partial-reference based deconvolution model, PREDE, integrating the non-negative matrix factorization framework with an iterative optimization strategy. We conducted comprehensive evaluations for PREDE using both simulation and real data analyses, demonstrating better performance of our method than other existing methods.

Reference-free deconvolution, visualization and interpretation of complex DNA methylation data using DecompPipeline, MeDeCom and FactorViz

DNA methylation profiling offers unique insights into human development and diseases. Often the analysis of complex tissues and cell mixtures is the only feasible option to study methylation changes across large patient cohorts. Since DNA methylomes are highly cell type specific, deconvolution methods can be used to recover cell type-specific information in the form of latent methylation components (LMCs) from such 'bulk' samples. Reference-free deconvolution methods retrieve these components without the need for DNA methylation profiles of purified cell types. Currently no integrated and guided procedure is available for data preparation and subsequent interpretation of deconvolution results. Here, we describe a three-stage protocol for reference-free deconvolution of DNA methylation data comprising: (i) data preprocessing, confounder adjustment using independent component analysis (ICA) and feature selection using DecompPipeline, (ii) deconvolution with multiple parameters using MeDeCom, RefFreeCellMix or EDec and (iii) guided biological inference and validation of deconvolution results with the R/Shiny graphical user interface FactorViz. Our protocol simplifies the analysis and guides the initial interpretation of DNA methylation data derived from complex samples. The harmonized approach is particularly useful to dissect and evaluate cell heterogeneity in complex systems such as tumors. We apply the protocol to lung cancer methylomes from The Cancer Genome Atlas (TCGA) and show that our approach identifies the proportions of stromal cells and tumor-infiltrating immune cells, as well as associations of the detected components with clinical parameters. The protocol takes slightly >3 d to complete and requires basic R skills.

Towards multi-omics characterization of tumor heterogeneity: a comprehensive review of statistical and machine learning approaches

The multi-omics molecular characterization of cancer opened a new horizon for our understanding of cancer biology and therapeutic strategies. However, a tumor biopsy comprises diverse types of cells limited not only to cancerous cells but also to tumor microenvironmental cells and adjacent normal cells. This heterogeneity is a major confounding factor that hampers a robust and reproducible bioinformatic analysis for biomarker identification using multi-omics profiles. Besides, the heterogeneity itself has been recognized over the years for its significant prognostic values in some cancer types, thus offering another promising avenue for therapeutic intervention. A number of computational approaches to unravel such heterogeneity from high-throughput molecular profiles of a tumor sample have been proposed, but most of them rely on the data from an individual omics layer. Since the heterogeneity of cells is widely distributed across multi-omics layers, methods based on an individual layer can only partially characterize the heterogeneous admixture of cells. To help facilitate further development of the methodologies that synchronously account for several multi-omics profiles, we wrote a comprehensive review of diverse approaches to characterize tumor heterogeneity based on three different omics layers: genome, epigenome and transcriptome. As a result, this review can be useful for the analysis of multi-omics profiles produced by many large-scale consortia. Contact:sunkim.bioinfo@snu.ac.kr.

Capturing functional epigenomes for insight into metabolic diseases

BACKGROUND

Metabolic diseases such as obesity are known to be driven by both environmental and genetic factors. Although genome-wide association studies of common variants and their impact on complex traits have provided some biological insight into disease etiology, identified genetic variants have been found to contribute only a small proportion to disease heritability, and to map mainly to non-coding regions of the genome. To link variants to function, association studies of cellular traits, such as epigenetic marks, in disease-relevant tissues are commonly applied.

SCOPE OF THE REVIEW

We review large-scale efforts to generate genome-wide maps of coordinated epigenetic marks and their utility in complex disease dissection with a focus on DNA methylation. We contrast DNA methylation profiling methods and discuss the advantages of using targeted methods for single-base resolution assessments of methylation levels across tissue-specific regulatory regions to deepen our understanding of contributing factors leading to complex diseases.

MAJOR CONCLUSIONS

Large-scale assessments of DNA methylation patterns in metabolic disease-linked study cohorts have provided insight into the impact of variable epigenetic variants in disease etiology. In-depth profiling of epigenetic marks at regulatory regions, particularly at tissue-specific elements, will be key to dissect the genetic and environmental components contributing to metabolic disease onset and progression.

Quantitative comparison of within-sample heterogeneity scores for DNA methylation data

DNA methylation is an epigenetic mark with important regulatory roles in cellular identity and can be quantified at base resolution using bisulfite sequencing. Most studies are limited to the average DNA methylation levels of individual CpGs and thus neglect heterogeneity within the profiled cell populations. To assess this within-sample heterogeneity (WSH) several window-based scores that quantify variability in DNA methylation in sequencing reads have been proposed. We performed the first systematic comparison of four published WSH scores based on simulated and publicly available datasets. Moreover, we propose two new scores and provide guidelines for selecting appropriate scores to address cell-type heterogeneity, cellular contamination and allele-specific methylation. Most of the measures were sensitive in detecting DNA methylation heterogeneity in these scenarios, while we detected differences in susceptibility to technical bias. Using recently published DNA methylation profiles of Ewing sarcoma samples, we show that DNA methylation heterogeneity provides information complementary to the DNA methylation level. WSH scores are powerful tools for estimating variance in DNA methylation patterns and have the potential for detecting novel disease-associated genomic loci not captured by established statistics. We provide an R-package implementing the WSH scores for integration into analysis workflows.

Computational methods in tumor immunology

The remarkable success of cancer immunotherapies, especially the checkpoint blocking antibodies, in a subset of patients has reinvigorated the study of tumor-immune crosstalk and its role in heterogeneity of response. High-throughput sequencing and imaging technologies can help recapitulate various aspects of the tumor ecosystem. Computational approaches provide an arsenal of tools to efficiently analyze, quantify and integrate multiple parameters of tumor immunity mined from these diverse but complementary high-throughput datasets. This chapter describes numerous such computational approaches in tumor immunology that leverage high-throughput data from diverse sources (genomic, transcriptomics, epigenomics and digitized histopathology images) to systematically interrogate tumor immunity in context of its microenvironment, and to identify mechanisms that confer resistance or sensitivity to cancer therapies, in particular immunotherapy.

Disease prediction by cell-free DNA methylation

Disease diagnosis using cell-free DNA (cfDNA) has been an active research field recently. Most existing approaches perform diagnosis based on the detection of sequence variants on cfDNA; thus, their applications are limited to diseases associated with high mutation rate such as cancer. Recent developments start to exploit the epigenetic information on cfDNA, which could have substantially wider applications. In this work, we provide thorough reviews and discussions on the statistical method developments and data analysis strategies for using cfDNA epigenetic profiles, in particular DNA methylation, to construct disease diagnostic models. We focus on two important aspects: marker selection and prediction model construction, under different scenarios. We perform simulations and real data analysis to compare different approaches, and provide recommendations for data analysis.

Guidelines for cell-type heterogeneity quantification based on a comparative analysis of reference-free DNA methylation deconvolution software

Background

Cell-type heterogeneity of tumors is a key factor in tumor progression and response to chemotherapy. Tumor cell-type heterogeneity, defined as the proportion of the various cell-types in a tumor, can be inferred from DNA methylation of surgical specimens. However, confounding factors known to associate with methylation values, such as age and sex, complicate accurate inference of cell-type proportions. While reference-free algorithms have been developed to infer cell-type proportions from DNA methylation, a comparative evaluation of the performance of these methods is still lacking.

Results

Here we use simulations to evaluate several computational pipelines based on the software packages MeDeCom, EDec, and RefFreeEWAS. We identify that accounting for confounders, feature selection, and the choice of the number of estimated cell types are critical steps for inferring cell-type proportions. We find that removal of methylation probes which are correlated with confounder variables reduces the error of inference by 30–35%, and that selection of cell-type informative probes has similar effect. We show that Cattell’s rule based on the scree plot is a powerful tool to determine the number of cell-types. Once the pre-processing steps are achieved, the three deconvolution methods provide comparable results. We observe that all the algorithms’ performance improves when inter-sample variation of cell-type proportions is large or when the number of available samples is large. We find that under specific circumstances the methods are sensitive to the initialization method, suggesting that averaging different solutions or optimizing initialization is an avenue for future research.

Conclusion

Based on the lessons learned, to facilitate pipeline validation and catalyze further pipeline improvement by the community, we develop a benchmark pipeline for inference of cell-type proportions and implement it in the R package medepir.

Virtual methylome dissection facilitated by single-cell analyses

Background

Numerous cell types can be identified within plant tissues and animal organs, and the epigenetic modifications underlying such enormous cellular heterogeneity are just beginning to be understood. It remains a challenge to infer cellular composition using DNA methylomes generated for mixed cell populations. Here, we propose a semi-reference-free procedure to perform virtual methylome dissection using the nonnegative matrix factorization (NMF) algorithm.

Results

In the pipeline that we implemented to predict cell-subtype percentages, putative cell-type-specific methylated (pCSM) loci were first determined according to their DNA methylation patterns in bulk methylomes and clustered into groups based on their correlations in methylation profiles. A representative set of pCSM loci was then chosen to decompose target methylomes into multiple latent DNA methylation components (LMCs). To test the performance of this pipeline, we made use of single-cell brain methylomes to create synthetic methylomes of known cell composition. Compared with highly variable CpG sites, pCSM loci achieved a higher prediction accuracy in the virtual methylome dissection of synthetic methylomes. In addition, pCSM loci were shown to be good predictors of the cell type of the sorted brain cells. The software package developed in this study is available in the GitHub repository (https://github.com/Gavin-Yinld).

Conclusions

We anticipate that the pipeline implemented in this study will be an innovative and valuable tool for the decoding of cellular heterogeneity.

Epigenome-wide association study of narcolepsy-affected lateral hypothalamic brains, and overlapping DNA methylation profiles between narcolepsy and multiple sclerosis

Narcolepsy with cataplexy is a sleep disorder caused by a deficiency in hypocretin neurons in the lateral hypothalamus (LH). Here we performed an epigenome-wide association study (EWAS) of DNA methylation for narcolepsy and replication analyses using DNA samples extracted from two brain regions: LH (Cases: N = 4; Controls: N = 4) and temporal cortex (Cases: N = 7; Controls: N = 7). Seventy-seven differentially methylated regions (DMRs) were identified in the LH analysis, with the top association of a DMR in the myelin basic protein (MBP) region. Only five DMRs were detected in the temporal cortex analysis. Genes annotated to LH DMRs were significantly associated with pathways related to fatty acid response or metabolism. Two additional analyses applying the EWAS data were performed: (1) investigation of methylation profiles shared between narcolepsy and other disorders and (2) an integrative analysis of DNA methylation data and a genome-wide association study for narcolepsy. The results of the two approaches, which included significant overlap of methylated positions associated with narcolepsy and multiple sclerosis, indicated that the two diseases may partly share their pathogenesis. In conclusion, DNA methylation in LH where loss of orexin-producing neurons occurs may play a role in the pathophysiology of the disease.

TOAST: improving reference-free cell composition estimation by cross-cell type differential analysis

In the analysis of high-throughput data from complex samples, cell composition is an important factor that needs to be accounted for. Except for a limited number of tissues with known pure cell type profiles, a majority of genomics and epigenetics data relies on the "reference-free deconvolution" methods to estimate cell composition. We develop a novel computational method to improve reference-free deconvolution, which iteratively searches for cell type-specific features and performs composition estimation. Simulation studies and applications to six real datasets including both DNA methylation and gene expression data demonstrate favorable performance of the proposed method. TOAST is available at https://bioconductor.org/packages/TOAST .

Computational deconvolution of transcriptomics data from mixed cell populations

Summary

Gene expression analyses of bulk tissues often ignore cell type composition as an important confounding factor, resulting in a loss of signal from lowly abundant cell types. In this review, we highlight the importance and value of computational deconvolution methods to infer the abundance of different cell types and/or cell type-specific expression profiles in heterogeneous samples without performing physical cell sorting. We also explain the various deconvolution scenarios, the mathematical approaches used to solve them and the effect of data processing and different confounding factors on the accuracy of the deconvolution results.

Contact

katleen.depreter@ugent.be.

Supplementary information

Supplementary data are available at Bioinformatics online.

CONFINED: distinguishing biological from technical sources of variation by leveraging multiple methylation datasets

Methylation datasets are affected by innumerable sources of variability, both biological (cell-type composition, genetics) and technical (batch effects). Here, we propose a reference-free method based on sparse canonical correlation analysis to separate the biological from technical sources of variability. We show through simulations and real data that our method, CONFINED, is not only more accurate than the state-of-the-art reference-free methods for capturing known, replicable biological variability, but it is also considerably more robust to dataset-specific technical variability than previous approaches. CONFINED is available as an R package as detailed at https://github.com/cozygene/CONFINED.

The role of epigenetics for understanding mental health difficulties and its implications for psychotherapy research

Many mental health difficulties have developmental origins. Understanding the mechanisms for how psychosocial experiences are biologically embedded and influence lifelong development is a key challenge for the mental health disciplines. In recent years, epigenetic processes have emerged as a potential mechanism mediating the long-lasting vulnerability following the experience of adversity. Animal models provide evidence that early-life adversity can produce enduring epigenetic modifications in the brain, which mediate disorder-like behaviours, and there is emerging evidence to support that environmental factors influence epigenetic processes in humans. The investigation of DNA methylation, a chemical modification of the DNA with a role in gene regulatory processes, is becoming increasingly popular in psychological studies. A particular interest for the psychotherapy field lies in the potential for psychological interventions to influence epigenetic processes. Hence, the focus of this review will be on studies that have investigated intervention-associated changes in DNA methylation. Results of the first few studies will be critically reviewed, and a model of how therapy-associated changes of DNA methylation in peripheral, non-brain tissue might be useful as epigenetic biomarkers of treatment outcome will be presented. PRACTITIONER POINTS: Many mental health difficulties have substantial developmental origin. Epigenetic processes have emerged as a potential mechanism mediating the long-term effects of early adversity Epigenetic refers to cellular mechanisms that control gene expression states, independent of changes to the underlying DNA sequence. The epigenome can be highly dynamic and potentially influenced by external factors A particular interest for the psychotherapy field lies in the potential for psychological interventions to influence epigenetic processes.

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.

Identification of differentially methylated cell types in epigenome-wide association studies

An outstanding challenge of epigenome-wide association studies (EWASs) performed in complex tissues is the identification of the specific cell type(s) responsible for the observed differential DNA methylation. Here we present a statistical algorithm called CellDMC ( https://github.com/sjczheng/EpiDISH ), which can identify differentially methylated positions and the specific cell type(s) driving the differential methylation. We validated CellDMC on in silico mixtures of DNA methylation data generated with different technologies, as well as on real mixtures from epigenome-wide association and cancer epigenome studies. CellDMC achieved over 90% sensitivity and specificity in scenarios where current state-of-the-art methods did not identify differential methylation. By applying CellDMC to an EWAS performed in buccal swabs, we identified smoking-associated differentially methylated positions occurring in the epithelial compartment, which we validated in smoking-related lung cancer. CellDMC may be useful in the identification of causal DNA-methylation alterations in disease.

BayesCCE: a Bayesian framework for estimating cell-type composition from DNA methylation without the need for methylation reference

We introduce a Bayesian semi-supervised method for estimating cell counts from DNA methylation by leveraging an easily obtainable prior knowledge on the cell-type composition distribution of the studied tissue. We show mathematically and empirically that alternative methods which attempt to infer cell counts without methylation reference only capture linear combinations of cell counts rather than provide one component per cell type. Our approach allows the construction of components such that each component corresponds to a single cell type, and provides a new opportunity to investigate cell compositions in genomic studies of tissues for which it was not possible before.

Epigenome-Wide Analyses Identify Two Novel Associations With Recurrent Stroke in the Vitamin Intervention for Stroke Prevention Clinical Trial

DNA methylation, a well-characterized epigenetic modification that is influenced by both environment and genetic variation, has previously been implicated in a number of complex diseases, including cardiovascular disease and stroke. The goal of this study was to evaluate epigenome-wide associations with recurrent stroke and the folate one-carbon metabolism-related trait, plasma homocysteine (hcy). Differential methylation analyses were performed on 473,864 autosomal CpG loci, using Illumina HumanMethylation 450K array data in 180 ischemic stroke cases from the Vitamin Intervention for Stroke Prevention (VISP) clinical trial. Linear regression was used to assess associations between number of strokes prior to VISP enrollment and measures of hcy with degree of methylation (β-values), while logistic regression was used to evaluate recurrent stroke status and incident recurrent stroke associations. All regression analyses were stratified by race. Two differentially methylated CpG sites exceeded epigenome-wide significance (p ≤ 1.055 × 10⁻⁷) for prior number of strokes (PNS) in European Americans. The top locus, cg22812874, was located in the ankyrin repeat and SOCS box containing 10 gene (ASB10; p = 3.4 × 10⁻⁹; β = −0.0308; 95% CI = −0.040, −0.002). Methylation locus cg00340919, located in an intron of the tetratricopeptide repeat domain 37 gene, was also statistically significant (TTC37; p = 8.74 × 10⁻⁸; β = −0.0517; 95% CI = −0.069, −0.034). An additional 138 CpG sites met our threshold for suggestive significance (p ≤ 5 × 10⁻⁵). We evaluated DNA methylation associated with recurrent stroke and hcy phenotypes across the epigenome. Hypermethylation at two CpG sites located in ASB10 and TTC37 was associated with fewer strokes prior to VISP enrollment. Our findings present a foundation for additional epigenome-wide studies, as well as mechanistic studies into epigenetic marks that influence recurrent stroke risk.

DNA methylation analysis on purified neurons and glia dissects age and Alzheimer's disease-specific changes in the human cortex

Background

Epigenome-wide association studies (EWAS) based on human brain samples allow a deep and direct understanding of epigenetic dysregulation in Alzheimer’s disease (AD). However, strong variation of cell-type proportions across brain tissue samples represents a significant source of data noise. Here, we report the first EWAS based on sorted neuronal and non-neuronal (mostly glia) nuclei from postmortem human brain tissues.

Results

We show that cell sorting strongly enhances the robust detection of disease-related DNA methylation changes even in a relatively small cohort. We identify numerous genes with cell-type-specific methylation signatures and document differential methylation dynamics associated with aging specifically in neurons such as CLU, SYNJ2 and NCOR2 or in glia RAI1,CXXC5 and INPP5A. Further, we found neuron or glia-specific associations with AD Braak stage progression at genes such as MCF2L, ANK1, MAP2, LRRC8B, STK32C and S100B. A comparison of our study with previous tissue-based EWAS validates multiple AD-associated DNA methylation signals and additionally specifies their origin to neuron, e.g., HOXA3 or glia (ANK1). In a meta-analysis, we reveal two novel previously unrecognized methylation changes at the key AD risk genes APP and ADAM17.

Conclusions

Our data highlight the complex interplay between disease, age and cell-type-specific methylation changes in AD risk genes thus offering new perspectives for the validation and interpretation of large EWAS results.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0211-3) contains supplementary material, which is available to authorized users.

A novel cell-type deconvolution algorithm reveals substantial contamination by immune cells in saliva, buccal and cervix

Aim: An outstanding challenge in epigenome studies is the estimation of cell-type proportions in complex epithelial tissues. Materials & methods: Here, we construct and validate a DNA methylation reference and algorithm for complex tissues that contain epithelial, immune and nonimmune stromal cells. Results: Using this reference, we show that easily accessible tissues such as saliva, buccal and cervix exhibit substantial variation in immune cell (IC) contamination. We further validate our reference in the context of oral cancer, where it correctly predicts an increased IC infiltration in cancer but suppressed in patients with highest smoking exposure. Finally, our method can improve the specificity of differentially methylated CpG calls in epithelial cancer. Conclusion: The degree and variation of IC contamination in complex epithelial tissues is substantial. We provide a valuable resource and tool for assessing the epithelial purity and IC contamination of samples and for identifying differential methylation in such complex tissues.

MethCNA: a database for integrating genomic and epigenomic data in human cancer

Background

The integration of DNA methylation and copy number alteration data promises to provide valuable insight into the underlying molecular mechanisms responsible for cancer initiation and progression. However, the generation and processing of these datasets are costly and time-consuming if carried out separately. The Illumina Infinium HumanMethylation450 BeadChip, initially designed for the evaluation of DNA methylation levels, allows copy number variant calling using bioinformatics tools.

Results

A substantial amount of Infinium HumanMethylation450 data across various cancer types has been accumulated in recent years and is a valuable resource for large-scale data analysis. Here we present MethCNA, a comprehensive database for genomic and epigenomic data integration in human cancer. In the current release, MethCNA contains about 10,000 tumor samples representing 37 cancer types. All raw array data were collected from The Cancer Genome Atlas and NCBI Gene Expression Omnibus database and analyzed using a pipeline that integrated multiple computational resources and tools. The normalized copy number aberration data and DNA methylation alterations were obtained. We provide a user-friendly web-interface for data mining and visualization.

Conclusions

The Illumina Infinium HumanMethylation450 BeadChip enables the interrogation and integration of both genomic and epigenomic data from exactly the same DNA specimen, and thus can aid in distinguishing driver from passenger mutations in cancer. We expect MethCNA will enable researchers to explore DNA methylation and copy number alteration patterns, identify key oncogenic drivers in cancer, and assist in the development of targeted therapies. MethCNA is publicly available online at http://cgma.scu.edu.cn/MethCNA.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4525-0) contains supplementary material, which is available to authorized users.