Epigenetic polymorphism and the stochastic formation of differentially methylated regions in normal and cancerous tissues

Circulating tumor DNA to guide diagnosis and treatment of localized and locally advanced non-small cell lung cancer

Liquid biopsy is a minimally invasive method for biomarkers detection in body fluids, particularly in blood, which offers an elevated and growing number of clinical applications in oncology. As a result of the improvement in the techniques for DNA analysis, above all next-generation sequencing (NGS) assays, circulating tumor DNA (ctDNA) has become the most informing tumor-derived material for most types of cancer, including non-small cell lung cancer (NSCLC). Although ctDNA concentration is higher in patients with advanced tumors, it can be detected even in patients with early-stage disease. Therefore, numerous clinical applications of ctDNA in the management of early-stage lung cancer are emerging, such as lung cancer screening, the identification of minimal residual disease (MRD), and the prediction of relapse before radiologic progression. Moreover, a high number of clinical trials are ongoing to better define the impact of ctDNA evaluation in this setting. Aim of this review is to offer a comprehensive overview of the most relevant implementations in using ctDNA for the management of early-stage lung cancer, addressing available data, technical aspects, limitations, and future perspectives.

Single-molecule epiallelic profiling of DNA derived from routinely collected Pap specimens for noninvasive detection of ovarian cancer

Recent advances in molecular analyses of ovarian cancer have revealed a wealth of promising tumour-specific biomarkers, including protein, DNA mutations and methylation; however, reliably detecting such alterations at satisfactorily high sensitivity and specificity through low-cost methods remains challenging, especially in early-stage diseases. Here we present PapDREAM, a new approach that enables detection of rare, ovarian-cancer-specific aberrations of DNA methylation from routinely-collected cervical Pap specimens. The PapDREAM approach employs a microfluidic platform that performs highly parallelized digital high-resolution melt to analyze locus-specific DNA methylation patterns on a molecule-by-molecule basis at or near single CpG-site resolution at a fraction (< 1/10th) of the cost of next-generation sequencing techniques. We demonstrate the feasibility of the platform by assessing intermolecular heterogeneity of DNA methylation in a panel of methylation biomarker loci using DNA derived from Pap specimens obtained from a cohort of 43 women, including 18 cases with ovarian cancer and 25 cancer-free controls. PapDREAM leverages systematic multidimensional bioinformatic analyses of locus-specific methylation heterogeneity to improve upon Pap-specimen-based detection of ovarian cancer, demonstrating a clinical sensitivity of 50% at 99% specificity in detecting ovarian cancer cases with an area under the receiver operator curve of 0.90. We then establish a logistic regression model that could be used to identify high-risk patients for subsequent clinical follow-up and monitoring. The results of this study support the utility of PapDREAM as a simple, low-cost screening method with the potential to integrate with existing clinical workflows for early detection of ovarian cancer. KEY POINTS: We present a microfluidic platform for detection and analysis of rare, heterogeneously methylated DNA within Pap specimens towards detection of ovarian cancer. The platform achieves high sensitivity (fractions <0.00005%) at a suitably low cost (∼$25) for routine screening applications. Furthermore, it provides molecule-by-molecule quantitative analysis to facilitate further study on the effect of heterogeneous methylation on cancer development.

CelFiE-ISH: a probabilistic model for multi-cell type deconvolution from single-molecule DNA methylation haplotypes

Deconvolution methods infer quantitative cell type estimates from bulk measurement of mixed samples including blood and tissue. DNA methylation sequencing measures multiple CpGs per read, but few existing deconvolution methods leverage this within-read information. We develop CelFiE-ISH, which extends an existing method (CelFiE) to use within-read haplotype information. CelFiE-ISH outperforms CelFiE and other existing methods, achieving 30% better accuracy and more sensitive detection of rare cell types. We also demonstrate the importance of marker selection and of tailoring markers for haplotype-aware methods. While here we use gold-standard short-read sequencing data, haplotype-aware methods will be well-suited for long-read sequencing.

Age Prediction Using DNA Methylation Heterogeneity Metrics

Dynamic changes in genomic DNA methylation patterns govern the epigenetic developmental programs and accompany the organism's aging. Epigenetic clock (eAge) algorithms utilize DNA methylation to estimate the age and risk factors for diseases as well as analyze the impact of various interventions. High-throughput bisulfite sequencing methods, such as reduced-representation bisulfite sequencing (RRBS) or whole genome bisulfite sequencing (WGBS), provide an opportunity to identify the genomic regions of disordered or heterogeneous DNA methylation, which might be associated with cell-type heterogeneity, DNA methylation erosion, and allele-specific methylation. We systematically evaluated the applicability of five scores assessing the variability of methylation patterns by evaluating within-sample heterogeneity (WSH) to construct human blood epigenetic clock models using RRBS data. The best performance was demonstrated by the model based on a metric designed to assess DNA methylation erosion with an MAE of 3.686 years. We also trained a prediction model that uses the average methylation level over genomic regions. Although this region-based model was relatively more efficient than the WSH-based model, the latter required the analysis of just a few short genomic regions and, therefore, could be a useful tool to design a reduced epigenetic clock that is analyzed by targeted next-generation sequencing.

On epigenetic stochasticity, entropy and cancer risk

Epigenetic changes are known to accrue in normal cells as a result of ageing and cumulative exposure to cancer risk factors. Increasing evidence points towards age-related epigenetic changes being acquired in a quasi-stochastic manner, and that they may play a causal role in cancer development. Here, I describe the quasi-stochastic nature of DNA methylation (DNAm) changes in ageing cells as well as in normal cells at risk of neoplastic transformation, discussing the implications of this stochasticity for developing cancer risk prediction strategies, and in particular, how it may require a conceptual paradigm shift in how we select cancer risk markers. I also describe the mounting evidence that a significant proportion of DNAm changes in ageing and cancer development are related to cell proliferation, reflecting tissue-turnover and the opportunity this offers for predicting cancer risk via the development of epigenetic mitotic-like clocks. Finally, I describe how age-associated DNAm changes may be causally implicated in cancer development via an irreversible suppression of tissue-specific transcription factors that increases epigenetic and transcriptomic entropy, promoting a more plastic yet aberrant cancer stem-cell state. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Epigenomic Characterization of Lymphoid Neoplasms

Lymphoid neoplasms represent a heterogeneous group of disease entities and subtypes with markedly different molecular and clinical features. Beyond genetic alterations, lymphoid tumors also show widespread epigenomic changes. These severely affect the levels and distribution of DNA methylation, histone modifications, chromatin accessibility, and three-dimensional genome interactions. DNA methylation stands out as a tracer of cell identity and memory, as B cell neoplasms show epigenetic imprints of their cellular origin and proliferative history, which can be quantified by an epigenetic mitotic clock. Chromatin-associated marks are informative to uncover altered regulatory regions and transcription factor networks contributing to the development of distinct lymphoid tumors. Tumor-intrinsic epigenetic and genetic aberrations cooperate and interact with microenvironmental cells to shape the transcriptome at different phases of lymphoma evolution, and intraclonal heterogeneity can now be characterized by single-cell profiling. Finally, epigenetics offers multiple clinical applications, including powerful diagnostic and prognostic biomarkers as well as therapeutic targets.

Blood DNA methylation profiling identifies cathepsin Z dysregulation in pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is characterised by pulmonary vascular remodelling causing premature death from right heart failure. Established DNA variants influence PAH risk, but susceptibility from epigenetic changes is unknown. We addressed this through epigenome-wide association study (EWAS), testing 865,848 CpG sites for association with PAH in 429 individuals with PAH and 1226 controls. Three loci, at Cathepsin Z (CTSZ, cg04917472), Conserved oligomeric Golgi complex 6 (COG6, cg27396197), and Zinc Finger Protein 678 (ZNF678, cg03144189), reached epigenome-wide significance (p < 10-7) and are hypermethylated in PAH, including in individuals with PAH at 1-year follow-up. Of 16 established PAH genes, only cg10976975 in BMP10 shows hypermethylation in PAH. Hypermethylation at CTSZ is associated with decreased blood cathepsin Z mRNA levels. Knockdown of CTSZ expression in human pulmonary artery endothelial cells increases caspase-3/7 activity (p < 10-4). DNA methylation profiles are altered in PAH, exemplified by the pulmonary endothelial function modifier CTSZ, encoding protease cathepsin Z.

Estimating genome-wide DNA methylation heterogeneity with methylation patterns

BACKGROUND

In a heterogeneous population of cells, individual cells can behave differently and respond variably to the environment. This cellular diversity can be assessed by measuring DNA methylation patterns. The loci with variable methylation patterns are informative of cellular heterogeneity and may serve as biomarkers of diseases and developmental progression. Cell-to-cell methylation heterogeneity can be evaluated through single-cell methylomes or computational techniques for pooled cells. However, the feasibility and performance of these approaches to precisely estimate methylation heterogeneity require further assessment.

RESULTS

Here, we proposed model-based methods adopted from a mathematical framework originally from biodiversity, to estimate genome-wide DNA methylation heterogeneity. We evaluated the performance of our models and the existing methods with feature comparison, and tested on both synthetic datasets and real data. Overall, our methods have demonstrated advantages over others because of their better correlation with the actual heterogeneity. We also demonstrated that methylation heterogeneity offers an additional layer of biological information distinct from the conventional methylation level. In the case studies, we showed that distinct profiles of methylation heterogeneity in CG and non-CG methylation can predict the regulatory roles between genomic elements in Arabidopsis. This opens up a new direction for plant epigenomics. Finally, we demonstrated that our score might be able to identify loci in human cancer samples as putative biomarkers for early cancer detection.

CONCLUSIONS

We adopted the mathematical framework from biodiversity into three model-based methods for analyzing genome-wide DNA methylation heterogeneity to monitor cellular heterogeneity. Our methods, namely MeH, have been implemented, evaluated with existing methods, and are open to the research community.

Genome-wide single-molecule analysis of long-read DNA methylation reveals heterogeneous patterns at heterochromatin that reflect nucleosome organisation

High-throughput sequencing technology is central to our current understanding of the human methylome. The vast majority of studies use chemical conversion to analyse bulk-level patterns of DNA methylation across the genome from a population of cells. While this technology has been used to probe single-molecule methylation patterns, such analyses are limited to short reads of a few hundred basepairs. DNA methylation can also be directly detected using Nanopore sequencing which can generate reads measuring megabases in length. However, thus far these analyses have largely focused on bulk-level assessment of DNA methylation. Here, we analyse DNA methylation in single Nanopore reads from human lymphoblastoid cells, to show that bulk-level metrics underestimate large-scale heterogeneity in the methylome. We use the correlation in methylation state between neighbouring sites to quantify single-molecule heterogeneity and find that heterogeneity varies significantly across the human genome, with some regions having heterogeneous methylation patterns at the single-molecule level and others possessing more homogeneous methylation patterns. By comparing the genomic distribution of the correlation to epigenomic annotations, we find that the greatest heterogeneity in single-molecule patterns is observed within heterochromatic partially methylated domains (PMDs). In contrast, reads originating from euchromatic regions and gene bodies have more ordered DNA methylation patterns. By analysing the patterns of single molecules in more detail, we show the existence of a nucleosome-scale periodicity in DNA methylation that accounts for some of the heterogeneity we uncover in long single-molecule DNA methylation patterns. We find that this periodic structure is partially masked in bulk data and correlates with DNA accessibility as measured by nanoNOMe-seq, suggesting that it could be generated by nucleosomes. Our findings demonstrate the power of single-molecule analysis of long-read data to understand the structure of the human methylome.

Decoding Cancer Evolution: Integrating Genetic and Non-Genetic Insights

The development of cancer begins with cells transitioning from their multicellular nature to a state akin to unicellular organisms. This shift leads to a breakdown in the crucial regulators inherent to multicellularity, resulting in the emergence of diverse cancer cell subpopulations that have enhanced adaptability. The presence of different cell subpopulations within a tumour, known as intratumoural heterogeneity (ITH), poses challenges for cancer treatment. In this review, we delve into the dynamics of the shift from multicellularity to unicellularity during cancer onset and progression. We highlight the role of genetic and non-genetic factors, as well as tumour microenvironment, in promoting ITH and cancer evolution. Additionally, we shed light on the latest advancements in omics technologies that allow for in-depth analysis of tumours at the single-cell level and their spatial organization within the tissue. Obtaining such detailed information is crucial for deepening our understanding of the diverse evolutionary paths of cancer, allowing for the development of effective therapies targeting the key drivers of cancer evolution.

Quantitative Investigation of Methylation Heterogeneity by Digital Melting Curve Analysis on a SlipChip for Atrial Fibrillation

Methylation is an essential epigenetic modification involved in regulating gene expression and maintaining genome stability. Methylation patterns can be heterogeneous, exhibiting variations in both level and density. However, current methods of methylation analysis, including sequencing, methylation-specific PCR, and high-resolution melting curve analysis (HRM), face limitations of high cost, time-consuming workflows, and the difficulty of both accurate heterogeneity analysis and precise quantification. Here, a droplet array SlipChip-based (da-SlipChip-based) digital melting curve analysis (MCA) method was developed for the accurate quantification of both methylation level (ratio of methylated molecules to total molecules) and methylation density (ratio of methylated CpG sites to total CpG sites). The SlipChip-based digital MCA system supplements an in situ thermal cycler with a fluorescence imaging module for real-time MCA. The da-SlipChip can generate 10,656 droplets of 1 nL each, which can be separated into four lanes, enabling the simultaneous analysis of four samples. This method's clinical application was demonstrated by analyzing samples from ten healthy individuals and twenty patients with atrial fibrillation (AF), the most common arrhythmia. This method can distinguish healthy individuals from those with AF of both the paroxysmal and persistent types. It also holds potential for broader application in various research and clinical settings requiring methylation analysis.

There's more to it: uncovering genomewide DNA methylation heterogeneity

Tweetable abstract Monitoring changes in methylation heterogeneity can be powerful in detecting disease progression early. This editorial highlights the importance of profiling methylation heterogeneity and identifies existing measures and research gaps.

DRSN4mCPred: accurately predicting sites of DNA N4-methylcytosine using deep residual shrinkage network for diagnosis and treatment of gastrointestinal cancer in the precision medicine era

Introduction

The DNA N4-methylcytosine (4mC) site levels of those suffering from digestive system cancers were higher, and the pathogenesis of digestive system cancers may also be related to the changes in DNA 4mC levels. Identifying DNA 4mC sites is a very important step in studying the analysis of biological function and cancer prediction. Extracting accurate features from DNA sequences is the key to establishing a prediction model of effective DNA 4mC sites. This study sought to develop a new predictive model, DRSN4mCPred, which aimed to improve the performance of the predicting DNA 4mC sites.

Methods

The model adopted multi-scale channel attention to extract features and used attention feature fusion (AFF) to fuse features. In order to capture features information more accurately and effectively, this model utilized Deep Residual Shrinkage Network with Channel-Wise thresholds (DRSN-CW) to eliminate noise-related features and achieve a more precise feature representation, thereby, distinguishing the sites in DNA with 4mC and non-4mC. Additionally, the predictive model incorporated an inverted residual block, a Multi-scale Channel Attention Module (MS-CAM), a Bi-directional Long Short Term Memory Network (Bi-LSTM), AFF, and DRSN-CW.

Results and Discussion

The results indicated the predictive model DRSN4mCPred had extremely good performance in predicting the DNA 4mC sites across different species. This paper will potentially provide support for the diagnosis and treatment of gastrointestinal cancer based on artificial intelligence in the precise medical era.

DNA methylation entropy is associated with DNA sequence features and developmental epigenetic divergence

Epigenetic information defines tissue identity and is largely inherited in development through DNA methylation. While studied mostly for mean differences, methylation also encodes stochastic change, defined as entropy in information theory. Analyzing allele-specific methylation in 49 human tissue sample datasets, we find that methylation entropy is associated with specific DNA binding motifs, regulatory DNA, and CpG density. Then applying information theory to 42 mouse embryo methylation datasets, we find that the contribution of methylation entropy to time- and tissue-specific patterns of development is comparable to the contribution of methylation mean, and methylation entropy is associated with sequence and chromatin features conserved with human. Moreover, methylation entropy is directly related to gene expression variability in development, suggesting a role for epigenetic entropy in developmental plasticity.

Metheor: Ultrafast DNA methylation heterogeneity calculation from bisulfite read alignments

Phased DNA methylation states within bisulfite sequencing reads are valuable source of information that can be used to estimate epigenetic diversity across cells as well as epigenomic instability in individual cells. Various measures capturing the heterogeneity of DNA methylation states have been proposed for a decade. However, in routine analyses on DNA methylation, this heterogeneity is often ignored by computing average methylation levels at CpG sites, even though such information exists in bisulfite sequencing data in the form of phased methylation states, or methylation patterns. In this study, to facilitate the application of the DNA methylation heterogeneity measures in downstream epigenomic analyses, we present a Rust-based, extremely fast and lightweight bioinformatics toolkit called Metheor. As the analysis of DNA methylation heterogeneity requires the examination of pairs or groups of CpGs throughout the genome, existing softwares suffer from high computational burden, which almost make a large-scale DNA methylation heterogeneity studies intractable for researchers with limited resources. In this study, we benchmark the performance of Metheor against existing code implementations for DNA methylation heterogeneity measures in three different scenarios of simulated bisulfite sequencing datasets. Metheor was shown to dramatically reduce the execution time up to 300-fold and memory footprint up to 60-fold, while producing identical results with the original implementation, thereby facilitating a large-scale study of DNA methylation heterogeneity profiles. To demonstrate the utility of the low computational burden of Metheor, we show that the methylation heterogeneity profiles of 928 cancer cell lines can be computed with standard computing resources. With those profiles, we reveal the association between DNA methylation heterogeneity and various omics features. Source code for Metheor is at https://github.com/dohlee/metheor and is freely available under the GPL-3.0 license.

Gene Regulatory Network Inference Using Convolutional Neural Networks from scRNA-seq Data

In recent years, with the rapid development of single-cell sequencing technology, this brings new opportunities and challenges to reconstruct gene regulatory networks. On the one hand, scRNA-seq data reveal statistical information of gene expression at single-cell resolution, which is beneficial to construct gene expression regulatory networks. On the other hand, the noise and dropout of single-cell data bring great difficulties to the analysis of scRNA-seq data, resulting in lower accuracy of gene regulatory networks reconstructed by traditional methods. In this article, we propose a novel supervised convolutional neural network (CNNSE), which can extract gene expression information from 2D co-expression matrices of gene doublets and identify interactions between genes. Our method can avoid the loss of extreme point interference by constructing a 2D co-expression matrix of gene pairs and significantly improve the regulation precision between gene pairs. And the CNNSE model is able to obtain detailed and high-level semantic information from the 2D co-expression matrix. Our method achieves satisfactory results on simulated data [accuracy (ACC): 0.712, F1: 0.724]. On two real scRNA-seq datasets, our method exhibits higher stability and accuracy in inference tasks compared with other existing gene regulatory network inference algorithms.

Single-molecule footprinting identifies context-dependent regulation of enhancers by DNA methylation

Enhancers are cis-regulatory elements that control the establishment of cell identities during development. In mammals, enhancer activation is tightly coupled with DNA demethylation. However, whether this epigenetic remodeling is necessary for enhancer activation is unknown. Here, we adapted single-molecule footprinting to measure chromatin accessibility and transcription factor binding as a function of the presence of methylation on the same DNA molecules. We leveraged natural epigenetic heterogeneity at active enhancers to test the impact of DNA methylation on their chromatin accessibility in multiple cell lineages. Although reduction of DNA methylation appears dispensable for the activity of most enhancers, we identify a class of cell-type-specific enhancers where DNA methylation antagonizes the binding of transcription factors. Genetic perturbations reveal that chromatin accessibility and transcription factor binding require active demethylation at these loci. Thus, in addition to safeguarding the genome from spurious activation, DNA methylation directly controls transcription factor occupancy at active enhancers.

DNA methylation entropy as a measure of stem cell replication and aging

BACKGROUND

Epigenetic marks are encoded by DNA methylation and accumulate errors as organisms age. This drift correlates with lifespan, but the biology of how this occurs is still unexplained. We analyze DNA methylation with age in mouse intestinal stem cells and compare them to nonstem cells.

RESULTS

Age-related changes in DNA methylation are identical in stem and nonstem cells, affect most prominently CpG islands and correlate weakly with gene expression. Age-related DNA methylation entropy, measured by the Jensen-Shannon Distribution, affects up to 25% of the detectable CpG sites and is a better measure of aging than individual CpG methylation. We analyze this entropy as a function of age in seven other tissues (heart, kidney, skeletal muscle, lung, liver, spleen, and blood) and it correlates strikingly with tissue-specific stem cell division rates. Thus, DNA methylation drift and increased entropy with age are primarily caused by and are sensors for, stem cell replication in adult tissues.

CONCLUSIONS

These data have implications for the mechanisms of tissue-specific functional declines with aging and for the development of DNA-methylation-based biological clocks.

Epigenetics as a mediator of plasticity in cancer

The concept of an epigenetic landscape describing potential cellular fates arising from pluripotent cells, first advanced by Conrad Waddington, has evolved in light of experiments showing nondeterministic outcomes of regulatory processes and mathematical methods for quantifying stochasticity. In this Review, we discuss modern approaches to epigenetic and gene regulation landscapes and the associated ideas of entropy and attractor states, illustrating how their definitions are both more precise and relevant to understanding cancer etiology and the plasticity of cancerous states. We address the interplay between different types of regulatory landscapes and how their changes underlie cancer progression. We also consider the roles of cellular aging and intrinsic and extrinsic stimuli in modulating cellular states and how landscape alterations can be quantitatively mapped onto phenotypic outcomes and thereby used in therapy development.

Computational Methods for Single-cell DNA Methylome Analysis

Dissecting intercellular epigenetic differences is key to understanding tissue heterogeneity. Recent advances in single-cell DNA methylome profiling have presented opportunities to resolve this heterogeneity at the maximum resolution. While these advances enable us to explore frontiers of chromatin biology and better understand cell lineage relationships, they pose new challenges in data processing and interpretation. This review surveys the current state of computational tools developed for single-cell DNA methylome data analysis. We discuss critical components of single-cell DNA methylome data analysis, including data preprocessing, quality control, imputation, dimensionality reduction, cell clustering, supervised cell annotation, cell lineage reconstruction, gene activity scoring, and integration with transcriptome data. We also highlight unique aspects of single-cell DNA methylome data analysis and discuss how techniques common to other single-cell omics data analyses can be adapted to analyze DNA methylomes. Finally, we discuss existing challenges and opportunities for future development.

Non-Association of Driver Alterations in PTEN with Differential Gene Expression and Gene Methylation in IDH1 Wildtype Glioblastomas

During oncogenesis, alterations in driver genes called driver alterations (DAs) modulate the transcriptome, methylome and proteome through oncogenic signaling pathways. These modulatory effects of any DA may be analyzed by examining differentially expressed mRNAs (DEMs), differentially methylated genes (DMGs) and differentially expressed proteins (DEPs) between tumor samples with and without that DA. We aimed to analyze these modulations with 12 common driver genes in Isocitrate Dehydrogenase 1 wildtype glioblastomas (IDH1-W-GBs). Using Cbioportal, groups of tumor samples with and without DAs in these 12 genes were generated from the IDH1-W-GBs available from "The Cancer Genomics Atlas Firehose Legacy Study Group" (TCGA-FL-SG) on Glioblastomas (GBs). For all 12 genes, samples with and without DAs were compared for DEMs, DMGs and DEPs. We found that DAs in PTEN were unassociated with any DEM or DMG in contrast to DAs in all other drivers, which were associated with several DEMs and DMGs. This contrasting PTEN-related property of being unassociated with differential gene expression or methylation in IDH1-W-GBs was unaffected by concurrent DAs in other common drivers or by the types of DAs affecting PTEN. From the lists of DEMs and DMGs associated with some common drivers other than PTEN, enriched gene ontology terms and insights into the co-regulatory effects of these drivers on the transcriptome were obtained. The findings from this study can improve our understanding of the molecular mechanisms underlying gliomagenesis with potential therapeutic benefits.

Inferring the Cancer Cellular Epigenome Heterogeneity via DNA Methylation Patterns

Tumor cells evolve through space and time, generating genetically and phenotypically diverse cancer cell populations that are continually subjected to the selection pressures of their microenvironment and cancer treatment.

cvlr: finding heterogeneously methylated genomic regions using ONT reads

Summary

Nanopore reads encode information on the methylation status of cytosines in CpG dinucleotides. The length of the reads makes it comparatively easy to look at patterns consisting of multiple loci; here, we exploit this property to search for regions where one can define subpopulations of molecules based on methylation patterns. As an example, we run our clustering algorithm on known imprinted genes; we also scan chromosome 15 looking for windows corresponding to heterogeneous methylation. Our software can also compute the covariance of methylation across these regions while keeping into account the mixture of different types of reads.

Availability and implementation

https://github.com/EmanueleRaineri/cvlr.

Contact

simon.heath@cnag.crg.eu.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

MC profiling: a novel approach to analyze DNA methylation heterogeneity in genome-wide bisulfite sequencing data

DNA methylation is an epigenetic mark implicated in crucial biological processes. Most of the knowledge about DNA methylation is based on bulk experiments, in which DNA methylation of genomic regions is reported as average methylation. However, average methylation does not inform on how methylated cytosines are distributed in each single DNA molecule. Here, we propose Methylation Class (MC) profiling as a genome-wide approach to the study of DNA methylation heterogeneity from bulk bisulfite sequencing experiments. The proposed approach is built on the concept of MCs, groups of DNA molecules sharing the same number of methylated cytosines. The relative abundances of MCs from sequencing reads incorporates the information on the average methylation, and directly informs on the methylation level of each molecule. By applying our approach to publicly available bisulfite-sequencing datasets, we individuated cell-to-cell differences as the prevalent contributor to methylation heterogeneity. Moreover, we individuated signatures of loci undergoing imprinting and X-inactivation, and highlighted differences between the two processes. When applying MC profiling to compare different conditions, we identified methylation changes occurring in regions with almost constant average methylation. Altogether, our results indicate that MC profiling can provide useful insights on the epigenetic status and its evolution at multiple genomic regions.

iDNA-ABF: multi-scale deep biological language learning model for the interpretable prediction of DNA methylations

In this study, we propose iDNA-ABF, a multi-scale deep biological language learning model that enables the interpretable prediction of DNA methylations based on genomic sequences only. Benchmarking comparisons show that our iDNA-ABF outperforms state-of-the-art methods for different methylation predictions. Importantly, we show the power of deep language learning in capturing both sequential and functional semantics information from background genomes. Moreover, by integrating the interpretable analysis mechanism, we well explain what the model learns, helping us build the mapping from the discovery of important sequential determinants to the in-depth analysis of their biological functions.

Locally correlated kinetics of post-replication DNA methylation reveals processivity and region specificity in DNA methylation maintenance

DNA methylation occurs predominantly on cytosine-phosphate-guanine (CpG) dinucleotides in the mammalian genome, and the methylation landscape is maintained over mitotic cell division. It has been posited that coupling of maintenance methylation activity among neighbouring CpGs is critical to stability over cellular generations; however, the mechanism is unclear. We used mathematical models and stochastic simulation to analyse data from experiments that probe genome-wide methylation of nascent DNA post-replication in cells. We find that DNA methylation maintenance rates on individual CpGs are locally correlated, and the degree of this correlation varies by genomic regional context. By using theory of protein diffusion along DNA, we show that exponential decay of methylation rate correlation with genomic distance is consistent with enzyme processivity. Our results provide quantitative evidence of genome-wide methyltransferase processivity in vivo. We further developed a method to disentangle different mechanistic sources of kinetic correlations. From the experimental data, we estimate that an individual methyltransferase methylates neighbour CpGs processively if they are 36 basepairs apart, on average. But other mechanisms of coupling dominate for longer inter-CpG distances. Our study demonstrates that quantitative insights into enzymatic mechanisms can be obtained from replication-associated, cell-based genome-wide measurements, by combining data-driven statistical analyses with hypothesis-driven mathematical modelling.

White Paper: Mimetics of Class 2 Tumor Suppressor Proteins as Novel Drug Candidates for Personalized Cancer Therapy

Simple Summary

A concept is presented for a new therapeutic approach, still in its early stages, which focuses on the phenotypic mimicry (“mimesis”) of proteins encoded by highly disease-relevant class 2 tumor suppressor genes that are silenced by DNA promoter methylation. Proteins derived from tumor suppressor genes are usually considered control systems of cells against oncogenic properties. Thus they represent the brakes in the “car-of-life.” Restoring this “brake function” in tumors by administering mimetic drugs may have a significant therapeutic effect. The proposed approach could thus open up a new, hitherto unexploited area of research for the development of anticancer drugs for difficult-to-treat cancers.

Abstract

The aim of our proposed concept is to find new target structures for combating cancers with unmet medical needs. This, unfortunately, still applies to the majority of the clinically most relevant tumor entities such as, for example, liver cancer, pancreatic cancer, and many others. Current target structures almost all belong to the class of oncogenic proteins caused by tumor-specific genetic alterations, such as activating mutations, gene fusions, or gene amplifications, often referred to as cancer “driver alterations” or just “drivers.” However, restoring the lost function of tumor suppressor genes (TSGs) could also be a valid approach to treating cancer. TSG-derived proteins are usually considered as control systems of cells against oncogenic properties; thus, they represent the brakes in the “car-of-life.” Restoring these tumor-defective brakes by gene therapy has not been successful so far, with a few exceptions. It can be assumed that most TSGs are not being inactivated by genetic alteration (class 1 TSGs) but rather by epigenetic silencing (class 2 TSGs or short “C2TSGs”). Reactivation of C2TSGs in cancer therapy is being addressed by the use of DNA demethylating agents and histone deacetylase inhibitors which act on the whole cancer cell genome. These epigenetic therapies have neither been particularly successful, probably because they are “shotgun” approaches that, although acting on C2TSGs, may also reactivate epigenetically silenced oncogenic sequences in the genome. Thus, new strategies are needed to exploit the therapeutic potential of C2TSGs, which have also been named DNA methylation cancer driver genes or “DNAme drivers” recently. Here we present a concept for a new translational and therapeutic approach that focuses on the phenotypic imitation (“mimesis”) of proteins encoded by highly disease-relevant C2TSGs/DNAme drivers. Molecular knowledge on C2TSGs is used in two complementary approaches having the translational concept of defining mimetic drugs in common: First, a concept is presented how truncated and/or genetically engineered C2TSG proteins, consisting solely of domains with defined tumor suppressive function can be developed as biologicals. Second, a method is described for identifying small molecules that can mimic the effect of the C2TSG protein lost in the cancer cell. Both approaches should open up a new, previously untapped discovery space for anticancer drugs.

Epialleles and epiallelic heterogeneity in hematological malignancies

DNA methylation has a well-established role in the pathogenesis, prognosis, and response to treatment in all the spectra of hematological malignancies. However, most of the data reported involve average DNA methylation observed in a sample. The emergence of bisulfite sequencing methods such as enhanced reduced representation that permit analyze adjacent CpGs led to exciting findings. Among these are the epialleles shift and the resulting epigenetic heterogeneity observed in leukemias and lymphomas. Epialleles seem to have an influential role as the cause of mutations that characterize leukemias, may stratify groups with different prognosis and response to treatment, and may be redistributed in the genome at different time points of the disease promoting activation of alternate transcriptional networks. Epiallelic shift may be responsible for the intratumor heterogeneity observed within the cells of the same tumor which increases with disease aggressiveness. It may also responsible for the interpatient heterogeneity explaining why blood cancers exhibit different behavior among different patients. Understanding better epiallelic conformation and the consequent chromatin conformational changes and the pathways that may be affected will permit deeper understanding of hematological malignancies pathogenesis and treatment.

Interplay between chromatin marks in development and disease

DNA methylation (DNAme) and histone post-translational modifications (PTMs) have important roles in transcriptional regulation. Although many reports have characterized the functions of such chromatin marks in isolation, recent genome-wide studies reveal surprisingly complex interactions between them. Here, we focus on the interplay between DNAme and methylation of specific lysine residues on the histone H3 tail. We describe the impact of genetic perturbation of the relevant methyltransferases in the mouse on the landscape of chromatin marks as well as the transcriptome. In addition, we discuss the specific neurodevelopmental growth syndromes and cancers resulting from pathogenic mutations in the human orthologues of these genes. Integrating these observations underscores the fundamental importance of crosstalk between DNA and histone H3 methylation in development and disease.

Epigenetic Inheritance From Normal Origin Cells Can Determine the Aggressive Biology of Tumor-Initiating Cells and Tumor Heterogeneity

The acquisition of genetic- and epigenetic-abnormalities during transformation has been recognized as the two fundamental factors that lead to tumorigenesis and determine the aggressive biology of tumor cells. However, there is a regularity that tumors derived from less-differentiated normal origin cells (NOCs) usually have a higher risk of vascular involvement, lymphatic and distant metastasis, which can be observed in both lymphohematopoietic malignancies and somatic cancers. Obviously, the hypothesis of genetic- and epigenetic-abnormalities is not sufficient to explain how the linear relationship between the cellular origin and the biological behavior of tumors is formed, because the cell origin of tumor is an independent factor related to tumor biology. In a given system, tumors can originate from multiple cell types, and tumor-initiating cells (TICs) can be mapped to different differentiation hierarchies of normal stem cells, suggesting that the heterogeneity of the origin of TICs is not completely chaotic. TIC’s epigenome includes not only genetic- and epigenetic-abnormalities, but also established epigenetic status of genes inherited from NOCs. In reviewing previous studies, we found much evidence supporting that the status of many tumor-related “epigenetic abnormalities” in TICs is consistent with that of the corresponding NOC of the same differentiation hierarchy, suggesting that they may not be true epigenetic abnormalities. So, we speculate that the established statuses of genes that control NOC’s migration, adhesion and colonization capabilities, cell-cycle quiescence, expression of drug transporters, induction of mesenchymal formation, overexpression of telomerase, and preference for glycolysis can be inherited to TICs through epigenetic memory and be manifested as their aggressive biology. TICs of different origins can maintain different degrees of innate stemness from NOC, which may explain why malignancies with stem cell phenotypes are usually more aggressive.

Genetic variation at mouse and human ribosomal DNA influences associated epigenetic states

BACKGROUND

Ribosomal DNA (rDNA) displays substantial inter-individual genetic variation in human and mouse. A systematic analysis of how this variation impacts epigenetic states and expression of the rDNA has thus far not been performed.

RESULTS

Using a combination of long- and short-read sequencing, we establish that 45S rDNA units in the C57BL/6J mouse strain exist as distinct genetic haplotypes that influence the epigenetic state and transcriptional output of any given unit. DNA methylation dynamics at these haplotypes are dichotomous and life-stage specific: at one haplotype, the DNA methylation state is sensitive to the in utero environment, but refractory to post-weaning influences, whereas other haplotypes entropically gain DNA methylation during aging only. On the other hand, individual rDNA units in human show limited evidence of genetic haplotypes, and hence little discernible correlation between genetic and epigenetic states. However, in both species, adjacent units show similar epigenetic profiles, and the overall epigenetic state at rDNA is strongly positively correlated with the total rDNA copy number. Analysis of different mouse inbred strains reveals that in some strains, such as 129S1/SvImJ, the rDNA copy number is only approximately 150 copies per diploid genome and DNA methylation levels are < 5%.

CONCLUSIONS

Our work demonstrates that rDNA-associated genetic variation has a considerable influence on rDNA epigenetic state and consequently rRNA expression outcomes. In the future, it will be important to consider the impact of inter-individual rDNA (epi)genetic variation on mammalian phenotypes and diseases.

Optical Nanopore Sensors for Quantitative Analysis

Nanopore sensors have received significant interest for the detection of clinically important biomarkers with single-molecule resolution. These sensors typically operate by detecting changes in the ionic current through a nanopore due to the translocation of an analyte. Recently, there has been interest in developing optical readout strategies for nanopore sensors for quantitative analysis. This is because they can utilize wide-field microscopy to independently monitor many nanopores within a high-density array. This significantly increases the amount of statistics that can be obtained, thus enabling the analysis of analytes present at ultralow concentrations. Here, we review the use of optical nanopore sensing strategies for quantitative analysis. We discuss optical nanopore sensing assays that have been developed to detect clinically relevant biomarkers, the potential for multiplexing such measurements, and techniques to fabricate high density arrays of nanopores with a view toward the use of these devices for clinical applications.

A novel workflow for the qualitative analysis of DNA methylation data

•

A novel R package (EpiStatProfiler) for the qualitative analysis of DNA methylation data.

•

A novel workflow for the analysis of CG and non-CG epialleles starting from any type of bisulfite sequencing data.

•

EpiStatProfiler can perform strand-specific characterization of epialleles composition.

•

Important loci can be annotated along with their biological role and potential functions.

•

EpiStatProfiler has the ability to identify loci whose epiallelic profile is associated with disease pathogenesis.

DNA methylation is an epigenetic modification that plays a pivotal role in major biological mechanisms, such as gene regulation, genomic imprinting, and genome stability. Different combinations of methylated cytosines for a given DNA locus generate different epialleles and alterations of these latter have been associated with several pathological conditions. Existing computational methods and statistical tests relevant to DNA methylation analysis are mostly based on the comparison of average CpG sites methylation levels and they often neglect non-CG methylation. Here, we present EpiStatProfiler, an R package that allows the analysis of CpG and non-CpG based epialleles starting from bisulfite sequencing data through a collection of dedicated extraction functions and statistical tests. EpiStatProfiler is provided with a set of useful auxiliary features, such as customizable genomic ranges, strand-specific epialleles analysis, locus annotation and gene set enrichment analysis. We showcase the package functionalities on two public datasets by identifying putative relevant loci in mice harboring the Huntington’s disease-causing Htt gene mutation and in Ctcf +/− mice compared to their wild-type counterparts. To our knowledge, EpiStatProfiler is the first package providing functionalities dedicated to the analysis of epialleles composition derived from any kind of bisulfite sequencing experiment.

OUP accepted manuscript

Motivation

Intermediately methylated regions occupy a significant fraction of the human genome and are closely associated with epigenetic regulations or cell-type deconvolution of bulk data. However, these regions show distinct methylation patterns, corresponding to different biological mechanisms. Although there have been some metrics developed for investigating these regions, the high noise sensitivity limits the utility for distinguishing distinct methylation patterns.

Results

We proposed a method named MeConcord to measure local methylation concordance across reads and CpG sites, respectively. MeConcord showed the most stable performance in distinguishing distinct methylation patterns (‘identical’, ‘uniform’ and ‘disordered’) compared with other metrics. Applying MeConcord to the whole genome data across 25 cell lines or primary cells or tissues, we found that distinct methylation patterns were associated with different genomic characteristics, such as CTCF binding or imprinted genes. Further, we showed the differences of CpG island hypermethylation patterns between senescence and tumorigenesis by using MeConcord. MeConcord is a powerful method to study local read-level methylation patterns for both the whole genome and specific regions of interest.

Availability and implementation

MeConcord is available at https://github.com/WangLabTHU/MeConcord.

Supplementary information

Supplementary data are available at Bioinformatics online.

DNA methylation in genes associated with the evolution of ageing and disease: A critical review

Ageing is characterised by a physical decline in biological functioning which results in a progressive risk of mortality with time. As a biological phenomenon, it is underpinned by the dysregulation of a myriad of complex processes. Recently, however, ever-increasing evidence has associated epigenetic mechanisms, such as DNA methylation (DNAm) with age-onset pathologies, including cancer, cardiovascular disease, and Alzheimer's disease. These diseases compromise healthspan. Consequently, there is a medical imperative to understand the link between epigenetic ageing, and healthspan. Evolutionary theory provides a unique way to gain new insights into epigenetic ageing and health. This review will: (1) provide a brief overview of the main evolutionary theories of ageing; (2) discuss recent genetic evidence which has revealed alleles that have pleiotropic effects on fitness at different ages in humans; (3) consider the effects of DNAm on pleiotropic alleles, which are associated with age related disease; (4) discuss how age related DNAm changes resonate with the mutation accumulation, disposable soma and programmed theories of ageing; (5) discuss how DNAm changes associated with caloric restriction intersect with the evolution of ageing; and (6) conclude by discussing how evolutionary theory can be used to inform investigations which quantify age-related DNAm changes which are linked to age onset pathology.

A mathematical model which examines age-related stochastic fluctuations in DNA maintenance methylation

Due to its complexity and its ubiquitous nature the ageing process remains an enduring biological puzzle. Many molecular mechanisms and biochemical process have become synonymous with ageing. However, recent findings have pinpointed epigenetics as having a key role in ageing and healthspan. In particular age related changes to DNA methylation offer the possibility of monitoring the trajectory of biological ageing and could even be used to predict the onset of diseases such as cancer, Alzheimer's disease and cardiovascular disease. At the molecular level emerging evidence strongly suggests the regulatory processes which govern DNA methylation are subject to intracellular stochasticity. It is challenging to fully understand the impact of stochasticity on DNA methylation levels at the molecular level experimentally. An ideal solution is to use mathematical models to capture the essence of the stochasticity and its outcomes. In this paper we present a novel stochastic model which accounts for specific methylation levels within a gene promoter. Uncertainty of the eventual site-specific methylation levels for different values of methylation age, depending on the initial methylation levels were analysed. Our model predicts the observed bistable levels in CpG islands. In addition, simulations with various levels of noise indicate that uncertainty predominantly spreads through the hypermethylated region of stability, especially for large values of input noise. A key outcome of the model is that CpG islands with high to intermediate methylation levels tend to be more susceptible to dramatic DNA methylation changes due to increasing methylation age.

Genome-wide detection of cytosine methylations in plant from Nanopore data using deep learning

In plants, cytosine DNA methylations (5mCs) can happen in three sequence contexts as CpG, CHG, and CHH (where H = A, C, or T), which play different roles in the regulation of biological processes. Although long Nanopore reads are advantageous in the detection of 5mCs comparing to short-read bisulfite sequencing, existing methods can only detect 5mCs in the CpG context, which limits their application in plants. Here, we develop DeepSignal-plant, a deep learning tool to detect genome-wide 5mCs of all three contexts in plants from Nanopore reads. We sequence Arabidopsis thaliana and Oryza sativa using both Nanopore and bisulfite sequencing. We develop a denoising process for training models, which enables DeepSignal-plant to achieve high correlations with bisulfite sequencing for 5mC detection in all three contexts. Furthermore, DeepSignal-plant can profile more 5mC sites, which will help to provide a more complete understanding of epigenetic mechanisms of different biological processes.

RLM: Fast and simplified extraction of Read-Level Methylation metrics from bisulfite sequencing data

 

Bisulfite sequencing data provide value beyond the straightforward methylation assessment by analyzing single-read patterns. Over the past years, various informative metrics have been established to explore this information. However, limited compatibility with alignment tools, reference genomes or the measurements they provide present a bottleneck for most groups to include this information as standard analysis. To address this, we developed RLM, a fast and scalable tool for the computation of frequently used Read-Level Methylation statistics. RLM supports several common alignment tools, works independently of the reference genome and handles all frequently used sequencing experiment designs. RLM can process large input files with a billion reads in just a few hours on common workstations.

AVAILABILITY

https://github.com/sarahet/RLM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Epigenetic encoding, heritability and plasticity of glioma transcriptional cell states

Single-cell RNA sequencing has revealed extensive transcriptional cell state diversity in cancer, often observed independently of genetic heterogeneity, raising the central question of how malignant cell states are encoded epigenetically. To address this, here we performed multiomics single-cell profiling-integrating DNA methylation, transcriptome and genotype within the same cells-of diffuse gliomas, tumors characterized by defined transcriptional cell state diversity. Direct comparison of the epigenetic profiles of distinct cell states revealed key switches for state transitions recapitulating neurodevelopmental trajectories and highlighted dysregulated epigenetic mechanisms underlying gliomagenesis. We further developed a quantitative framework to directly measure cell state heritability and transition dynamics based on high-resolution lineage trees in human samples. We demonstrated heritability of malignant cell states, with key differences in hierarchal and plastic cell state architectures in IDH-mutant glioma versus IDH-wild-type glioblastoma, respectively. This work provides a framework anchoring transcriptional cancer cell states in their epigenetic encoding, inheritance and transition dynamics.

Single-cell multimodal glioma analyses identify epigenetic regulators of cellular plasticity and environmental stress response

Glioma intratumoral heterogeneity enables adaptation to challenging microenvironments and contributes to therapeutic resistance. We integrated 914 single-cell DNA methylomes, 55,284 single-cell transcriptomes and bulk multi-omic profiles across 11 adult IDH mutant or IDH wild-type gliomas to delineate sources of intratumoral heterogeneity. We showed that local DNA methylation disorder is associated with cell-cell DNA methylation differences, is elevated in more aggressive tumors, links with transcriptional disruption and is altered during the environmental stress response. Glioma cells under in vitro hypoxic and irradiation stress increased local DNA methylation disorder and shifted cell states. We identified a positive association between genetic and epigenetic instability that was supported in bulk longitudinally collected DNA methylation data. Increased DNA methylation disorder associated with accelerated disease progression and recurrently selected DNA methylation changes were enriched for environmental stress response pathways. Our work identified an epigenetically facilitated adaptive stress response process and highlights the importance of epigenetic heterogeneity in shaping therapeutic outcomes.

Age-related demethylation of the TDP-43 autoregulatory region in the human motor cortex

In amyotrophic lateral sclerosis (ALS), TAR DNA-binding protein 43 (TDP-43), which is encoded by TARDBP, forms aggregates in the motor cortex. This aggregate formation may be triggered by an increase in the TDP-43 level with aging. However, the amount of TDP-43 is autoregulated by alternative splicing of the TARDBP 3′UTR, and how this autoregulation is affected by aging remains to be elucidated. We found that DNA demethylation in the autoregulatory region in the TARDBP 3′UTR reduced alternative splicing and increased TARDBP mRNA expression. Furthermore, in the human motor cortex, we found that this region was demethylated with aging, resulting in increased expression of TARDBP mRNA. The acceleration of DNA demethylation in the motor cortex was associated with the age of ALS onset. In summary, the dysregulation of TDP-43 autoregulation by age-related DNA demethylation in the motor cortex may explain the contribution of aging and motor system selectivity in ALS.

In order to assess the effects of aging on the autoregulation of TAR DNA-binding protein 43 (TDP-43) and the potential effects of this on the role of TDP-43 in Amyotrophic Lateral Sclerosis (ALS), Koike et al examined post-mortem motor cortex tissue from ALS patients. They found that DNA demethylation in the autoregulatory region of the TARDBP 3′UTR, which encodes TDP-43, increased with age and was associated with the onset age of ALS and thus could be indicative of a role for dysregulation of TDP-43 autoregulation in ALS pathology.

DNA methylation landscapes of 1538 breast cancers reveal a replication-linked clock, epigenomic instability and cis-regulation

DNA methylation is aberrant in cancer, but the dynamics, regulatory role and clinical implications of such epigenetic changes are still poorly understood. Here, reduced representation bisulfite sequencing (RRBS) profiles of 1538 breast tumors and 244 normal breast tissues from the METABRIC cohort are reported, facilitating detailed analysis of DNA methylation within a rich context of genomic, transcriptional, and clinical data. Tumor methylation from immune and stromal signatures are deconvoluted leading to the discovery of a tumor replication-linked clock with genome-wide methylation loss in non-CpG island sites. Unexpectedly, methylation in most tumor CpG islands follows two replication-independent processes of gain (MG) or loss (ML) that we term epigenomic instability. Epigenomic instability is correlated with tumor grade and stage, TP53 mutations and poorer prognosis. After controlling for these global trans-acting trends, as well as for X-linked dosage compensation effects, cis-specific methylation and expression correlations are uncovered at hundreds of promoters and over a thousand distal elements. Some of these targeted known tumor suppressors and oncogenes. In conclusion, this study demonstrates that global epigenetic instability can erode cancer methylomes and expose them to localized methylation aberrations in-cis resulting in transcriptional changes seen in tumors.

The concurrence of DNA methylation and demethylation is associated with transcription regulation

The mammalian DNA methylome is formed by two antagonizing processes, methylation by DNA methyltransferases (DNMT) and demethylation by ten-eleven translocation (TET) dioxygenases. Although the dynamics of either methylation or demethylation have been intensively studied in the past decade, the direct effects of their interaction on gene expression remain elusive. Here, we quantify the concurrence of DNA methylation and demethylation by the percentage of unmethylated CpGs within a partially methylated read from bisulfite sequencing. After verifying ‘methylation concurrence’ by its strong association with the co-localization of DNMT and TET enzymes, we observe that methylation concurrence is strongly correlated with gene expression. Notably, elevated methylation concurrence in tumors is associated with the repression of 40~60% of tumor suppressor genes, which cannot be explained by promoter hypermethylation alone. Furthermore, methylation concurrence can be used to stratify large undermethylated regions with negligible differences in average methylation into two subgroups with distinct chromatin accessibility and gene regulation patterns. Together, methylation concurrence represents a unique methylation metric important for transcription regulation and is distinct from conventional metrics, such as average methylation and methylation variation.

The global pattern of the mammalian methylome is formed by changes in methylation and demethylation. Here the authors describe a metric methylation concurrence that measures the ratio of unmethylated CpGs inside the partially methylated reads and show that methylation concurrence is associated with epigenetically regulated tumour suppressor genes.

The DNA methylation landscape of multiple myeloma shows extensive inter- and intrapatient heterogeneity that fuels transcriptomic variability

BACKGROUND

Cancer evolution depends on epigenetic and genetic diversity. Historically, in multiple myeloma (MM), subclonal diversity and tumor evolution have been investigated mostly from a genetic perspective.

METHODS

Here, we performed an analysis of 42 MM samples from 21 patients by using enhanced reduced representation bisulfite sequencing (eRRBS). We combined several metrics of epigenetic heterogeneity to analyze DNA methylation heterogeneity in MM patients.

RESULTS

We show that MM is characterized by the continuous accumulation of stochastic methylation at the promoters of development-related genes. High combinatorial entropy change is associated with poor outcomes in our pilot study and depends predominantly on partially methylated domains (PMDs). These PMDs, which represent the major source of inter- and intrapatient DNA methylation heterogeneity in MM, are linked to other key epigenetic aberrations, such as CpG island (CGI)/transcription start site (TSS) hypermethylation and H3K27me3 redistribution as well as 3D organization alterations. In addition, transcriptome analysis revealed that intratumor methylation heterogeneity was associated with low-level expression and high variability.

CONCLUSIONS

We propose that disrupted DNA methylation in MM is responsible for high epigenetic and transcriptomic instability allowing tumor cells to adapt to environmental changes by tapping into a pool of evolutionary trajectories.

Evaluation of whole-genome DNA methylation sequencing library preparation protocols

Background

With rapidly dropping sequencing cost, the popularity of whole-genome DNA methylation sequencing has been on the rise. Multiple library preparation protocols currently exist. We have performed 22 whole-genome DNA methylation sequencing experiments on snap frozen human samples, and extensively benchmarked common library preparation protocols for whole-genome DNA methylation sequencing, including three traditional bisulfite-based protocols and a new enzyme-based protocol. In addition, different input DNA quantities were compared for two kits compatible with a reduced starting quantity. In addition, we also present bioinformatic analysis pipelines for sequencing data from each of these library types.

Results

An assortment of metrics were collected for each kit, including raw read statistics, library quality and uniformity metrics, cytosine retention, and CpG beta value consistency between technical replicates. Overall, the NEBNext Enzymatic Methyl-seq and Swift Accel-NGS Methyl-Seq kits performed quantitatively better than the other two protocols. In addition, the NEB and Swift kits performed well at low-input amounts, validating their utility in applications where DNA is the limiting factor.

Results

The NEBNext Enzymatic Methyl-seq kit appeared to be the best option for whole-genome DNA methylation sequencing of high-quality DNA, closely followed by the Swift kit, which potentially works better for degraded samples. Further, a general bioinformatic pipeline is applicable across the four protocols, with the exception of extra trimming needed for the Swift Biosciences’s Accel-NGS Methyl-Seq protocol to remove the Adaptase sequence.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00401-y.

Discovery of Candidate DNA Methylation Cancer Driver Genes

Epigenetic alterations, such as promoter hypermethylation, may drive cancer through tumor suppressor gene inactivation. However, we have limited ability to differentiate driver DNA methylation (DNAme) changes from passenger events. We developed DNAme driver inference-MethSig-accounting for the varying stochastic hypermethylation rate across the genome and between samples. We applied MethSig to bisulfite sequencing data of chronic lymphocytic leukemia (CLL), multiple myeloma, ductal carcinoma in situ, glioblastoma, and to methylation array data across 18 tumor types in TCGA. MethSig resulted in well-calibrated quantile-quantile plots and reproducible inference of likely DNAme drivers with increased sensitivity/specificity compared with benchmarked methods. CRISPR/Cas9 knockout of selected candidate CLL DNAme drivers provided a fitness advantage with and without therapeutic intervention. Notably, DNAme driver risk score was closely associated with adverse outcome in independent CLL cohorts. Collectively, MethSig represents a novel inference framework for DNAme driver discovery to chart the role of aberrant DNAme in cancer. SIGNIFICANCE: MethSig provides a novel statistical framework for the analysis of DNA methylation changes in cancer, to specifically identify candidate DNA methylation driver genes of cancer progression and relapse, empowering the discovery of epigenetic mechanisms that enhance cancer cell fitness.This article is highlighted in the In This Issue feature, p. 2113.