DNA methylation is a fundamental mechanism of epigenetic control in cells and its dysregulation is strongly implicated in cancer development. The addition of a methyl group to cytosine, particularly at CpG islands (CGIs) results in changes in gene expression and affects the transcription of genomic repeats and non-coding RNA. These result in changes within cancer cells as well as the tumor microenvironment (TME) that favor the development of disease. Patterns of tumor-specific methylation are recapitulated in cell-free DNA (cfDNA) allowing use of the latter for the study of cancer epigenetics. Liquid biopsies can be used to obtain cell-free DNA from different biofluids such as cerebrospinal fluid, pleural effusion, urine, and plasma. Of the different sources of cfDNA, plasma has seen extensive use as a source of methylated DNA for use in biomarker studies. Diagnostic, prognostic, predictive, and treatment response monitoring applications have or are being developed around methylated cfDNA obtained from plasma.
2. Biology of DNA Methylation in Malignancy
2.1. Aberrant DNA Methylation of CpG Islands Are Seen in Cancer
More than half of all human promoters overlap with 1–2 kb CG-rich regions, referred to as CpG islands (CGIs). Nearly 95% of CGI promoter regions are unmethylated in normal tissues. While CGIs are preferentially found near promoters and transcription start sites (TSS), they may also appear within gene regulatory elements, gene bodies, and intergenic regions
[1]. Cancers typically display focal regions of tumor-specific hypermethylation, especially at CGIs, against the backdrop of global hypomethylation that impacts gene expression and carcinogenesis
[2][3][4].
2.2. Gene Promoter Hypermethylation Leads to Silencing of Tumor Suppressor Genes
Hypermethylation of transcriptional regulatory regions of tumor suppressor genes (TSGs) results in their silencing and the abrogation of their anti-tumor function
[2][5]. First described in colorectal cancer (CRC), the CpG island methylator phenotype (CIMP) that is characterized by the expansive TSG promoter silencing has since been documented in other solid malignancies
[6]. DNA hypermethylation was observed in promoters of classic TSGs such as retinoblastoma transcriptional corepressor 1 (
RB1) and cell cycle regulators
CDKN2A, among many others
[1][2].
2.3. Methylation Can Lead to Epigenetic Activation of Oncogenic Pathways
Hypomethylation of oncogene promoters typically results in their expression. Oncogene products dysregulate cell proliferation, migration, and invasion. This results in elevated metastatic potential, increased stemness in cells, or immune evasion
[7]. For example, BCL6 acts as a transcriptional regulator that represses tumor suppressors thereby favoring cell proliferation and survival
[8]. Expression of BCL6 is epigenetically regulated and knockdown of UHRF1, an upstream regulator, results in loss of
BCL6 gene promoter methylation thereby favoring its expression. Notably, BCL6 was observed to be overexpressed in exhausted T-cells and may play a role in tumor immune escape
[9]. In addition to promoters, hypomethylation of other gene regulatory regions, such as enhancers, was shown to correlate with the expression of cancer-specific genes and pathways across multiple cancers
[10][11].
In contrast to promoter hypermethylation, gene body hypermethylation is posited to correlate positively with gene expression
[1][12][13]. As an example, the homeobox superfamily of transcription factors is essential to the regulation of cell growth and differentiation. These genes have high CpG density, which makes their expression susceptible to epigenetic control
[14]. In relation to this, Su et al.
[15] reported that gene-body hypermethylation results in the overexpression of approximately 43% of known homeobox genes, many of which are associated with oncogenesis.
2.4. Global DNA Hypomethylation Leads to Chromosomal Instability
DNA methylation plays a critical role in maintaining genomic integrity by regulating the transcription of constitutive heterochromatin elements, sub-telomeric regions, and genomic repeats
[16][17]. Expansive hypomethylation is common in cancer and may involve regions as large as ten megabases in size, which can render the cancer genome susceptible to genomic instability. Moreover, increased gene transcription from loss of gene promoter methylation results in the formation of DNA-RNA R-loops that can cause replication stress and DNA damage—further exacerbating genome instability
[16][18]. Reduced methylation destabilizes heterochromatin at pericentromeric regions and is associated with loss of tumor suppression and oncogenesis in gastrointestinal cancers and sarcomas
[16][17]. Telomere length is affected by methylation of sub-telomeric regions, and its dysregulation is considered a cancer hallmark
[16][19][20].
2.5. DNA Methylation Influences Gene Expression beyond Promoter Silencing
Alternative mRNA splicing is responsible for transcriptional diversity in cells
[21][22]. During splicing, non-coding introns are removed, and coding exons are ligated to form mature mRNA at specific positions called splice sites. Recognition of introns and exons by the spliceosome machinery is critical for the correct execution of this process
[21]. DNA methylation modifies the accessibility of exonic nucleosomes to RNA Polymerase II, which demonstrates the regulatory role it plays in splicing
[21][22].
Chen and Elnitski
[23] showed that methylation-correlated isoforms affect the functional protein domains of gene products. Moreover, these gene sets were enriched for oncogenes, tumor suppressors, and cancer-related pathways across multiple cancers
[23].
2.6. Methylation Can Influence the Non-Coding Genome, Leading to Cancer
DNA hypomethylation led to the expression of genomic repeat elements and transposable elements that interfered with chromosomal integrity and drove oncogene expression in a process called onco-exaptation
[16][17]. Jang et al.
[24] reported RNA-seq data from 7769 tumors across 15 cancer types from the TCGA. In this cohort, the expression of 106 oncogenes was regulated by onco-exaptation in nearly half (49.7%) of tumors. To validate this relationship, the authors performed in vitro induction of DNA methylation on the promoter of the short-interspersed nuclear element (SINE)
AluJb by CRISPR. This resulted in a 20–30% increase in methylation and consequently a 40% reduction in the expression of the oncogene
LIN28B. In addition to onco-exaptation, transposon-associated dsRNA-dependent paracrine signaling had been found to be favorable to tumor growth in breast, lung, and pancreatic cancer
[25][26].
The study of non-coding RNA is a growing area of interest owing to the emerging roles that it plays in disease
[27][28]. The interaction between non-coding RNA (ncRNA) and DNA methylation is complex. Few studies have comprehensively examined if these non-coding regions are methylated, with most investigations focusing on a specific type of ncRNA called long non-coding RNA (lncRNA).
Several comprehensive reviews discuss the association between DNA methylation and lncRNA expression in cancer
[29][30]. In contrast to lncRNAs, other types of noncoding RNA have not been examined as extensively in the context of DNA methylation. This underscores the complexity of the relationship between ncRNA and DNA methylation in cancer
[31].
2.7. Role of DNA Methylation in the Tumor Microenvironment
Methylation dynamics play a key role in contextualizing the tumor microenvironment (TME). To illustrate, TET2 demethylates promoters of genes encoding cytokines and transcription factors that influence CD4
+ T-cell fate. Expression of the cytokine IL-4 favors differentiation into Th2 cells while expression of the transcriptional regulator FOXP3 expression leads to differentiation into tissue regulatory T-cells (Treg)
[32]. In CD8
+ cells, active demethylation of gene enhancer regions by TET2 facilitates differentiation to T-effector cells. Conversely, the conditional loss of TET activity was shown to drive a shift towards a memory T-cell fate
[32][33][34]. Meanwhile, DNMT3A-directed de novo methylation of the promoter of T-cell-specific transcription factor 7 (
Tcf7) results in the suppression of memory T-cell differentiation and supports differentiation into effector subtypes
[33][34].
As with T-cells, myeloid differentiation is epigenetically influenced by DNA methylation. TET2 acts in concert with IL-4 and STAT6 to promote the expression of ITGB2, which is a cell adhesion molecule important in monocyte-to-dendritic cell differentiation and function. Dendritic cells have a central role in antigen presentation and T-cell activation within the TME
[32]. Myeloid-derived suppressor cells (MDSC) are key immune regulators in the TME that serve a pro-tumorigenic function by promoting immune tolerance to cancer. An increase in DNMT3A expression in MDSCs in the presence of Prostaglandin E2 results in a gain of DNA methylation and silencing of immunogenicity-associated genes (e.g.,
FAS, RUNX1, and
S1PR4)
[35][36]. Finally, macrophage polarization towards either the activated (M1) phenotype or the inhibitory (M2) phenotype is also subject to the influence of DNA methylation. Express of DNMT3b and DNMT1 are associated with M1-like macrophage polarization through its effect on PPAR γ1
[37]. Meanwhile, siRNA knockdown of the
DNMT3B gene increases M2 marker expression in macrophages. A similar increase in M2 markers was seen following the treatment of M1 macrophages with DNMT inhibitors (DNMTi)
[37][38][39].
The development of immune escape mechanisms that bar cytotoxic effector cells from infiltrating the TME is key to tumor progression
[20]. DNA methylation also plays a part in this process. DNMT1 activity was shown to disturb the trafficking of CD8
+ T-cells into the TME of a murine ovarian cancer cell model by repressing the tumor production of CXCL9 and CXCL10
[32]. In addition, tumors use DNA methylation to negatively regulate chemokines such as CCL5 and CCL2, which are essential in T-cell and macrophage chemotaxis. Epigenetic silencing of the
CCL2 gene by methylation was found to diminish macrophage infiltration and promote disease progression in small-cell lung cancer (SCLC)
[40][41][42]. As described previously, DNA methylation also plays a role in the activation of inhibitory immune cells such as MDSC that modulate a pro-tumor TME. Interestingly, MDSC infiltration can be reduced by treatment with DNMTi
[43].
Investigations into the mechanisms behind effector T-cell exhaustion revealed the contribution of multiple epigenetic mechanisms including DNA methylation. De novo methylation by DNMT3A established a stable exhausted state and its deletion in CAR T-cells as shown by Prinzing et al.
[44] led to a more anti-tumor state. Moreover, DNA methylation profiling of exhausted tumor-infiltrating lymphocytes (TIL) identified promoter hypomethylation of T-cell exhaustion markers PD-1 and HAVCR2. Binding sites of exhaustion-associated transcription factors such as NR4A1 were also seen to be hypomethylated, implying that activation of the T-cell inhibitory programming was governed, at least in part, by methylation
[45]. Taken together, these highlight the complex roles DNA methylation plays in the tumor TME.
Figure 1 summarizes the disparate roles DNA methylation plays in the development of cancer.
Figure 1. Influence of DNA methylation on the tumor and tumor microenvironment. Cancers are characterized by extensive genomic hypomethylation that leads to chromosomal instability and oncogene activation. Focal hypermethylation of CpG islands at promoters leads to TSG silencing. Hypermethylation of gene bodies can lead to oncogene expression and alternative splicing leading to cancer-related isoforms of genes. In the tumor microenvironment DNA methylation can lead to altered immune cell fate and trafficking. DNA methylation also underpins epigenetic mechanisms of T-cell exhaustion. Created with Biorender.
3. Cell-Free Tumor Methylome Recapitulates the Cancer Epigenome
Cell-free DNA (cfDNA) refers to DNA fragments in the noncellular component of blood. Although it has been an area of great interest over the last decade, Mandel and Métais were the first to demonstrate DNA in the blood of both healthy and diseased patients in 1948
[46]. cfDNA arises from cells through necrosis and apoptosis, but also enters the circulation via active secretion
[46][47]. Circulating DNA from normal cells is found in low concentrations in plasma (~10–15 ng/mL) but can increase due to acute stressors like disease and injury
[48]. Notably, it was shown that levels of cfDNA in individuals with cancer are higher in comparison with their counterparts without cancer
[48]. The component of cfDNA that originates from tumors is referred to as circulating tumor DNA (ctDNA) and represents a smaller fraction of overall cfDNA
[49][50]. These circulating tumor fragments are shorter than DNA from non-cancer cells and range in size from ~140 to 200 bp in length, corresponding to the size of nucleosome-associated DNA
[46][51]. The ctDNA fraction in plasma varies but increases in association with disease type and burden. Advanced disease and high-shedding tumors such as SCLC and prostate cancer may have fractions > 20% of total cfDNA
[52][53][54][55]. Conversely, reduced ctDNA fractions in the blood are detectable in early-stage disease or following effective anticancer therapy
[48][52][54][56]. In addition to plasma, ctDNA can be obtained from other biofluids such as cerebrospinal fluid, saliva, pleural effusions, urine, ascites, and stool representing multiple sources of nucleic acid for study
[48]. Depending on the source, ctDNA size may vary because of disparate conditions within the respective biological compartments from which the nucleic acids are sourced. For instance, shorter ctDNA fragments may be observed in urine compared to plasma owing to the greater nuclease activity in the former
[46].
While the term liquid biopsy is most often used to refer to the analysis of cfDNA from peripheral blood, it also pertains to the isolation and study of tumor-derived material from other bodily fluids. Liquid biopsies using standard venipuncture afford several distinct advantages over conventional tissue sampling. First, their minimally invasive nature allows for evaluation of the disease, particularly in cases where tumor access is difficult, unsafe, or ethically challenging, especially for repeated longitudinal sampling
[48]. Moreover, in addition to ctDNA, circulating tumor microRNA, extracellular vesicles, tumor-associated proteins, and circulating tumor cells (CTCs) can also be sampled at the point of care by this approach
[57]. From a testing standpoint, utilizing pooled cfDNA obtained from liquid biopsies to determine a positive signal—such as the presence of cancer in a patient—avoids dependence on specific fragments and thereby improves sensitivity. This is especially important when dealing with scenarios wherein the amount of input DNA is constrained, and adequate representation of individual DNA fragments cannot be guaranteed
[58]. Another advantage of liquid biopsies is their ability to better recapitulate the molecular heterogeneity harbored by multiple distinct clonal populations that collectively shed cfDNA into the blood. This is because the molecular diversity of disease is not as easily captured by needle aspirates of single tumor sites and multiple needle passes to sample different sites may be difficult or excessively invasive
[48]. Finally, from a therapeutic perspective, this technique allows for longitudinal sampling across successive time points providing opportunities for studying tumor evolution resulting from treatment
[48][50].
Early applications of liquid biopsy technology involve genotyping and mutation detection, but indications for its use have since expanded to include epigenetic analyses, including the study of ctDNA methylation
[50][59][60][61]. The study of ctDNA methylation in turn has several advantages. First, cancer-specific methylation changes occur earlier in tumorigenesis and represent an opportunity for early disease detection and diagnosis
[49][62]. Moreover, combining methylation with orthogonal methods such as mutation detection and fragmentomics to profile ctDNA further enhances the discriminative power of such tests by interrogating biologic features that are characteristic of tumor-derived DNA
[63]. Second, methylation patterns in ctDNA are consistent with their tissue-of-origin, allowing for the location of the source of a cancer, especially those with an unknown primary site
[60][64]. Third, tumor cells tend to exhibit more homogenous DNA methylation changes versus gene mutations that tend to exhibit greater intra-tumoral heterogeneity
[59]. The consistency of the methylation signal opens the door for evaluating tumor evolution, especially after being subjected to the selective pressure of anti-cancer therapy. In the context of treatment resistance, examining changes in differentially methylated regions (DMR) in ctDNA obtained from longitudinal blood draws can shed light on epigenetic mechanisms that underpin treatment resistance. Various bisulfite conversion (e.g., quantitative methylation-specific PCR and cfMethyl-seq) and antibody-enrichment-based (e.g., cfMeDIP-seq) platforms for the study of ctDNA methylation have been developed in the last decade. The details of these assays are beyond the scope of this research but were extensively discussed elsewhere
[65][66][67][68].
4. Applications of cfDNA Methylation as Biomarker in Cancer
4.1. Methylation as a Diagnostic Biomarker
The concordance between tumor and cell-free cancer methylomes means that it can be leveraged as a tool for cancer diagnosis. Assays that employ ctDNA methylation for this application vary in their approaches. Some use targeted detection of methylated sites in specific genes such as
SEPT9 [69][70],
SHOX2 [71],
PTGER4 [71], and
SDC2 [72] while others profile more extensive swathes of the genome.
Classifiers based on more extensive DNA methylation profiles to detect cancers via plasma-based liquid biopsies were also developed. In this regard, the most significant application of ctDNA methylome profiling is in multicancer early detection (MCED). The importance of this application lies in the fact that early diagnosis of disease allows for the administration of effective, curative treatment to optimize patient outcomes
[73].
A second platform, this time that leveraged cfMeDIP-seq, an antibody-based, non-degradative genome-wide DNA methylation enrichment assay was being developed for use in MCED. Validation of this test was performed on a cohort of 4332 consisting of untreated patients across eight different tumor types—bladder, breast, CRC, head, and neck (HNSCC), lung, ovarian, prostate, and renal cell (RCC) cancers, and a set of age and sex-matched controls
[74]. The test performed well with an AUROC of 0.94 (95% CI: 0.93, 0.96) across all tumors and was robust in each individual tumor type including low cfDNA shedding tumors—bladder, breast, prostate, and RCC—with an AUROC 0.92 (0.91, 0.94)
[74]. The high detection of low-shedding disease was particularly noteworthy since early-stage cancers with low tumor burden released small amounts of cfDNA into the circulation and were an important hurdle to overcome for viable MCED tests.
4.2. Methylation as a Prognostic Biomarker of Patient Outcomes
As a ubiquitous mode of epigenetic regulation, DNA methylation of specific genes provides value in prognosticating outcomes in cancer. Huang et al.
[75] utilized bisulfite conversion-based real-time PCR to determine the methylation status of the SEPT9 gene (
mSEPT9) in the plasma of 144 preoperative CRC patients. Patients who were
mSEPT9+ had lower disease-free survival (DFS) rates than those who were
mSEPT9− (two-year DFS: 52.1% vs. 73.9%,
p = 0.014), Moreover,
mSEPT9 was an independent predictor of prognosis (HR = 2.741,
p = 0.009) in multivariate regression analysis
[75]. Of interest, a 2021 systematic review by Hier et al.
[76] identified over 100 different genes regulated by methylation that carried prognostic significance. Apart from protein-coding genes, the loss of methylation of repeat elements was significant in early carcinogenesis and correlated with poor prognostic outcomes
[77][78].
Recognizing that the tumor methylation signal was distinct from that of the non-cancer genome, Ul Haq et al.
[79] performed a whole genome approach to the assessment of the cell-free methylome of SCLC using cfMeDIP-seq. The authors utilized an in silico approach using paired peripheral blood leukocytes to subtract non-cancer noise from cfDNA methylation data similar to previous studies
[63]. This approach identified two biologically distinct methylation-defined SCLC prognostic groups with differences in overall survival
[79]. Validation of ctDNA methylation as a prognostic biomarker in HNSCC and RCC was performed and data from these studies were presented in abstract form at the ESMO 2023 conference
[80][81].
4.3. Methylation as a Predictive Biomarker of Treatment Response
Methylation of the O
6-alkylguanine DNA alkyltransferase (
MGMT) promoter had been a known predictor of response to the alkylator temozolomide in glioma
[57]. The promoter methylation status of this gene could be assessed using ctDNA obtained from blood or cerebrospinal fluid via targeted, methylation-specific PCR and was found to be highly concordant with tumor tissue
[82]. Meanwhile, the expression of the DNA/RNA helicase Schlafen 11 (SLFN11) predicted sensitivity to DNA-damaging alkylating agents and PARP inhibitors
[83][84][85][86]. The
SLFN11 gene was observed to be under epigenetic control and methylation of its promoter in ctDNA predicted resistance to therapy in ovarian cancer with the prior mentioned agents
[87]. Hypermethylation of the promoter of the
MLH1 gene represented a principal epigenetic mechanism that led to the MSI-H/dMMR state. Tumors with deficient mismatch repair possessed a higher burden of somatic mutations, higher infiltrating T-cell counts, and showed increased PD-L1 positivity—all these predispose MSI-H/dMMR cancers to respond to immune checkpoint inhibitor (ICI) therapy
[88][89]. Wang et al.
[90] described a liquid biopsy-based approach to the detection of
MLH1 promoter methylation in CRC using methylation-sensitive restriction enzyme PCR. This assay had an AUROC value of 0.965 (95% CI: 0.94–0.99). The sensitivity and specificity of the assay were 78% and 100%, respectively (95% CI: 0.45–0.95).
Beyond individual genes, methylation-profile-based scoring also predicted anti-cancer therapy response. In an exploratory study from the phase 2 SWOG S1314 clinical trial, pre-treatment cfDNA from 72 patients was used to generate a classifier methylation-based response score (mR-score), which could predict pathologic response to neoadjuvant chemotherapy in operable, muscle-invasive bladder cancer (MIBC). The model demonstrated an AUROC of 0.636 (95% CI: 0.498–0.773) for baseline samples. Moreover, mR-scores correlated with pathologic response with complete responders demonstrating the lowest scores in the cohort and non-responders showing the highest mR-scores. Notably, this same classifier showed better performance when applied to on-treatment samples (AUROC of 0.720, 95% CI: 0.582–0.857)
[91]. Methylation-profile-based scoring was also used to predict therapy response in ICI.
4.4. Minimal Residual Disease
Effective oncologic therapy reduces tumor burden by >2–3 log
10 which corresponds to a tumor kill of >99%. Unfortunately, malignant cells may persist in patients who achieve exquisite treatment responses and represent niduses for disease relapse. These recalcitrant cancer cells, referred to as minimal residual disease (MRD), remain a key challenge in the treatment of cancer
[92]. To guarantee the durability of remission following curative-intent therapy, increasingly sensitive methods for post-treatment surveillance are necessary
[50]. Knowing that recurrent tumors detectable by conventional imaging contain >10
6 cancer cells, the likelihood of successful treatment becomes less likely once it is demonstrable by diagnostic scans
[93]. Detection of ctDNA by liquid biopsies of different biofluids, such as plasma and urine, represents an avenue for earlier detection of recurrent disease thereby facilitating cancer interception, earlier intervention, and more optimal treatment outcomes
[94][95][96][97][98][99].