Myotonic dystrophy type 1 (DM1) clinical subtypes and CTCF site methylation status flanking the CTG expansion are mutant allele length-dependent.

Myotonic dystrophy type 1 (DM1) is a complex disease with a wide spectrum of symptoms. The exact relationship between mutant CTG repeat expansion size and clinical outcome remains unclear. DM1 congenital patients (CDM) inherit the largest expanded alleles, which are associated with abnormal and increased DNA methylation flanking the CTG repeat. However, DNA methylation at the DMPK locus remains understudied. Its relationship to DM1 clinical subtypes, expansion size and age-at-onset is not yet completely understood. Using pyrosequencing-based methylation analysis on 225 blood DNA samples from Costa Rican DM1 patients, we determined that the size of the estimated progenitor allele length (ePAL) is not only a good discriminator between CDM and non-CDM cases (with an estimated threshold at 653 CTG repeats), but also for all DM1 clinical subtypes. Secondly, increased methylation at both CTCF sites upstream and downstream of the expansion was almost exclusively present in CDM cases. Thirdly, levels of abnormal methylation were associated with clinical subtype, age and ePAL, with strong correlations between these variables. Fourthly, both ePAL and the intergenerational expansion size were significantly associated with methylation status. Finally, methylation status was associated with ePAL and maternal inheritance, with almost exclusively maternal transmission of CDM. In conclusion, increased DNA methylation at the CTCF sites flanking the DM1 expansion could be linked to ePAL, and both increased methylation and the ePAL could be considered biomarkers for the CDM phenotype.

Longitudinal increases in somatic mosaicism of the expanded CTG repeat in myotonic dystrophy type 1 are associated with variation in age-at-onset.

In myotonic dystrophy type 1 (DM1), somatic mosaicism of the (CTG)n repeat expansion is age-dependent, tissue-specific and expansion-biased. These features contribute toward variation in disease severity and confound genotype-to-phenotype analyses. To investigate how the (CTG)n repeat expansion changes over time, we collected three longitudinal blood DNA samples separated by 8-15 years and used small pool and single-molecule PCR in 43 DM1 patients. We used the lower boundary of the allele length distribution as the best estimate for the inherited progenitor allele length (ePAL), which is itself the best predictor of disease severity. Although in most patients the lower boundary of the allele length distribution was conserved over time, in many this estimate also increased with age, suggesting samples for research studies and clinical trials should be obtained as early as possible. As expected, the modal allele length increased over time, driven primarily by ePAL, age-at-sampling and the time interval. As expected, small expansions <100 repeats did not expand as rapidly as larger alleles. However, the rate of expansion of very large alleles was not obviously proportionally higher. This may, at least in part, be a result of the allele length-dependent increase in large contractions that we also observed. We also determined that individual-specific variation in the increase of modal allele length over time not accounted for by ePAL, age-at-sampling and time was inversely associated with individual-specific variation in age-at-onset not accounted for by ePAL, further highlighting somatic expansion as a therapeutic target in DM1.

DNMT3L facilitates DNA methylation partly by maintaining DNMT3A stability in mouse embryonic stem cells.

DNMT3L (DNMT3-like), a member of the DNMT3 family, has no DNA methyltransferase activity but regulates de novo DNA methylation. While biochemical studies show that DNMT3L is capable of interacting with both DNMT3A and DNMT3B and stimulating their enzymatic activities, genetic evidence suggests that DNMT3L is essential for DNMT3A-mediated de novo methylation in germ cells but is dispensable for de novo methylation during embryogenesis, which is mainly mediated by DNMT3B. How DNMT3L regulates DNA methylation and what determines its functional specificity are not well understood. Here we show that DNMT3L-deficient mouse embryonic stem cells (mESCs) exhibit downregulation of DNMT3A, especially DNMT3A2, the predominant DNMT3A isoform in mESCs. DNA methylation analysis of DNMT3L-deficient mESCs reveals hypomethylation at many DNMT3A target regions. These results confirm that DNMT3L is a positive regulator of DNA methylation, contrary to a previous report that, in mESCs, DNMT3L regulates DNA methylation positively or negatively, depending on genomic regions. Mechanistically, DNMT3L forms a complex with DNMT3A2 and prevents DNMT3A2 from being degraded. Restoring the DNMT3A protein level in DNMT3L-deficient mESCs partially recovers DNA methylation. Thus, our work uncovers a role for DNMT3L in maintaining DNMT3A stability, which contributes to the effect of DNMT3L on DNMT3A-dependent DNA methylation.

Age-related epigenetic drift in the pathogenesis of MDS and AML.

The myelodysplastic syndrome (MDS) is a clonal hematologic disorder that frequently evolves to acute myeloid leukemia (AML). Its pathogenesis remains unclear, but mutations in epigenetic modifiers are common and the disease often responds to DNA methylation inhibitors. We analyzed DNA methylation in the bone marrow and spleen in two mouse models of MDS/AML, the NUP98-HOXD13 (NHD13) mouse and the RUNX1 mutant mouse model. Methylation array analysis showed an average of 512/3445 (14.9%) genes hypermethylated in NHD13 MDS, and 331 (9.6%) genes hypermethylated in RUNX1 MDS. Thirty-two percent of genes in common between the two models (2/3 NHD13 mice and 2/3 RUNX1 mice) were also hypermethylated in at least two of 19 human MDS samples. Detailed analysis of 41 genes in mice showed progressive drift in DNA methylation from young to old normal bone marrow and spleen; to MDS, where we detected accelerated age-related methylation; and finally to AML, which markedly extends DNA methylation abnormalities. Most of these genes showed similar patterns in human MDS and AML. Repeat element hypomethylation was rare in MDS but marked the transition to AML in some cases. Our data show consistency in patterns of aberrant DNA methylation in human and mouse MDS and suggest that epigenetically, MDS displays an accelerated aging phenotype.

Defects of Lipid Synthesis Are Linked to the Age-Dependent Demyelination Caused by Lamin B1 Overexpression.

Lamin B1 is a component of the nuclear lamina and plays a critical role in maintaining nuclear architecture, regulating gene expression and modulating chromatin positioning. We have previously shown that LMNB1 gene duplications cause autosomal dominant leukodystrophy (ADLD), a fatal adult onset demyelinating disease. The mechanisms by which increased LMNB1 levels cause ADLD are unclear. To address this, we used a transgenic mouse model where Lamin B1 overexpression is targeted to oligodendrocytes. These mice showed severe vacuolar degeneration of the spinal cord white matter together with marked astrogliosis, microglial infiltration, and secondary axonal damage. Oligodendrocytes in the transgenic mice revealed alterations in histone modifications favoring a transcriptionally repressed state. Chromatin changes were accompanied by reduced expression of genes involved in lipid synthesis pathways, many of which are known to play important roles in myelin regulation and are preferentially expressed in oligodendrocytes. Decreased lipogenic gene expression resulted in a significant reduction in multiple classes of lipids involved in myelin formation. Many of these gene expression changes and lipid alterations were observed even before the onset of the phenotype, suggesting a causal role. Our findings establish, for the first time, a link between LMNB1 and lipid synthesis in oligodendrocytes, and provide a mechanistic framework to explain the age dependence and white matter involvement of the disease phenotype. These results have implications for disease pathogenesis and may also shed light on the regulation of lipid synthesis pathways in myelin maintenance and turnover. SIGNIFICANCE STATEMENT: Autosomal dominant leukodystrophy (ADLD) is fatal neurological disorder caused by increased levels of the nuclear protein, Lamin B1. The disease is characterized by an age-dependent loss of myelin, the fatty sheath that covers nerve fibers. We have studied a mouse model where Lamin B1 level are increased in oligodendrocytes, the cell type that produces myelin in the CNS. We demonstrate that destruction of myelin in the spinal cord is responsible for the degenerative phenotype in our mouse model. We show that this degeneration is mediated by reduced expression of lipid synthesis genes and the subsequent reduction in myelin enriched lipids. These findings provide a mechanistic framework to explain the age dependence and tissue specificity of the ADLD disease phenotype.

CCAT2, a novel noncoding RNA mapping to 8q24, underlies metastatic progression and chromosomal instability in colon cancer.

The functional roles of SNPs within the 8q24 gene desert in the cancer phenotype are not yet well understood. Here, we report that CCAT2, a novel long noncoding RNA transcript (lncRNA) encompassing the rs6983267 SNP, is highly overexpressed in microsatellite-stable colorectal cancer and promotes tumor growth, metastasis, and chromosomal instability. We demonstrate that MYC, miR-17-5p, and miR-20a are up-regulated by CCAT2 through TCF7L2-mediated transcriptional regulation. We further identify the physical interaction between CCAT2 and TCF7L2 resulting in an enhancement of WNT signaling activity. We show that CCAT2 is itself a WNT downstream target, which suggests the existence of a feedback loop. Finally, we demonstrate that the SNP status affects CCAT2 expression and the risk allele G produces more CCAT2 transcript. Our results support a new mechanism of MYC and WNT regulation by the novel lncRNA CCAT2 in colorectal cancer pathogenesis, and provide an alternative explanation of the SNP-conferred cancer risk.

TET1 is a maintenance DNA demethylase that prevents methylation spreading in differentiated cells.

TET1 is a 5-methylcytosine dioxygenase and its DNA demethylating activity has been implicated in pluripotency and reprogramming. However, the precise role of TET1 in DNA methylation regulation outside of developmental reprogramming is still unclear. Here, we show that overexpression of the TET1 catalytic domain but not full length TET1 (TET1-FL) induces massive global DNA demethylation in differentiated cells. Genome-wide mapping reveals that 5-hydroxymethylcytosine production by TET1-FL is inhibited as DNA methylation increases, which can be explained by the preferential binding of TET1-FL to unmethylated CpG islands (CGIs) through its CXXC domain. TET1-FL specifically accumulates 5-hydroxymethylcytosine at the edges of hypomethylated CGIs, while knockdown of endogenous TET1 induces methylation spreading from methylated edges into hypomethylated CGIs. We also found that gene expression changes after TET1-FL overexpression are relatively small and independent of its dioxygenase function. Thus, our results identify TET1 as a maintenance DNA demethylase that does not purposely decrease methylation levels, but specifically prevents aberrant methylation spreading into CGIs in differentiated cells.

RUNX3 promoter hypermethylation is frequent in leukaemia cell lines and associated with acute myeloid leukaemia inv(16) subtype.

Correlative and functional studies support the involvement of the RUNX gene family in haematological malignancies. To elucidate the role of epigenetics in RUNX inactivation, we evaluated promoter DNA methylation of RUNX1, 2, and 3 in 23 leukaemia cell lines and samples from acute myeloid leukaemia (AML), acute lymphocytic leukaemia (ALL) and myelodysplatic syndromes (MDS) patients. RUNX1 and RUNX2 gene promoters were mostly unmethylated in cell lines and clinical samples. Hypermethylation of RUNX3 was frequent among cell lines (74%) and highly variable among patient samples, with clear association to cytogenetic status. High frequency of RUNX3 hypermethylation (85% of the 20 studied cases) was found in AML patients with inv(16)(p13.1q22) compared to other AML subtypes (31% of the other 49 cases). RUNX3 hypermethylation was also frequent in ALL (100% of the six cases) but low in MDS (21%). In support of a functional role, hypermethylation of RUNX3 was correlated with low levels of protein, and treatment of cell lines with the DNA demethylating agent, decitabine, resulted in mRNA re-expression. Furthermore, relapse-free survival of non-inv(16)(p13.1q22) AML patients without RUNX3 methylation was significantly better (P = 0·016) than that of methylated cases. These results suggest that RUNX3 silencing is an important event in inv(16)(p13.1q22) leukaemias.

Colorectal carcinomas with CpG island methylator phenotype 1 frequently contain mutations in chromatin regulators.

BACKGROUND & AIMS: Subgroups of colorectal carcinomas (CRCs) characterized by DNA methylation anomalies are termed CpG island methylator phenotype (CIMP)1, CIMP2, or CIMP-negative. The pathogenesis of CIMP1 colorectal carcinomas, and their effects on patients' prognoses and responses to treatment, differ from those of other CRCs. We sought to identify genetic somatic alterations associated with CIMP1 CRCs. METHODS: We examined genomic DNA samples from 100 primary CRCs, 10 adenomas, and adjacent normal-appearing mucosae from patients undergoing surgery or colonoscopy at 3 tertiary medical centers. We performed exome sequencing of 16 colorectal tumors and their adjacent normal tissues. Extensive comparison with known somatic alterations in CRCs allowed segregation of CIMP1-exclusive alterations. The prevalence of mutations in selected genes was determined from an independent cohort. RESULTS: We found that genes that regulate chromatin were mutated in CIMP1 CRCs; the highest rates of mutation were observed in CHD7 and CHD8, which encode members of the chromodomain helicase/adenosine triphosphate-dependent chromatin remodeling family. Somatic mutations in these 2 genes were detected in 5 of 9 CIMP1 CRCs. A prevalence screen showed that nonsilencing mutations in CHD7 and CHD8 occurred significantly more frequently in CIMP1 tumors (18 of 42 [43%]) than in CIMP2 (3 of 34 [9%]; P < .01) or CIMP-negative tumors (2 of 34 [6%]; P < .001). CIMP1 markers had increased binding by CHD7, compared with all genes. Genes altered in patients with CHARGE syndrome (congenital malformations involving the central nervous system, eye, ear, nose, and mediastinal organs) who had CHD7 mutations were also altered in CRCs with mutations in CHD7. CONCLUSIONS: Aberrations in chromatin remodeling could contribute to the development of CIMP1 CRCs. A better understanding of the biological determinants of CRCs can be achieved when these tumors are categorized according to their epigenetic status.

Regulation of AURKC expression by CpG island methylation in human cancer cells.

AURKC, a member of the Aurora kinase gene family, is highly expressed in testis but is either moderately expressed or repressed in most somatic cells. Varying expression of AURKC has been observed in human cancers, but the underlying mechanisms of differential expression have been investigated only to a limited extent. We investigated the role of promoter CpG methylation in the regulation of AURKC gene expression in human cancer cells, in relation to a recently reported AURKC transcription repressor PLZF/ZBTB16, implicated in transformation and tumorigenesis. AURKC and PLZF/ZBTB16 expression profiles were investigated in reference to CpG methylation status on the AURKC promoter experimentally, and also in The Cancer Genome Atlas (TCGA) dataset involving multiple cancer types. AURKC promoter showed dense to moderate hypermethylation correlating with low to moderate expression of the gene in normal somatic cells and cancer cell lines, while testis with high expression revealed marked hypo-methylation. Treatment with the demethylating agent, 5-aza-dC, but not the histone deacetylase (HDAC) inhibitor, TSA, led to elevated expression in cancer cell lines, indicating that promoter DNA methylation negatively regulates AURKC expression. High expression of PLZF in PLZF-transfected cells treated with 5-aza-dC only partially repressed expression of AURKC despite 5-aza-dC also inducing elevated PLZF expression. Analyses of the TCGA data showed differential expression of AURKC in multiple cancer types and stronger correlation of AURKC expression with CpG methylation compared to PLZF levels. These findings demonstrate that differential promoter CpG methylation is an important mechanism regulating AURKC expression in cancer cells.

Epigenetic mechanisms of induced pluripotency.

Reprogramming of somatic cells to induced pluripotent stem cells (iPSCs) requires profound alterations in the epigenetic landscape. During reprogramming, a change in chromatin structure resets the gene expression and stabilises self-renewal. Reprogramming is a highly inefficient process, in part due to multiple epigenetic barriers. Although many epigenetic factors have already been shown to affect self-renewal and pluripotency in embryonic stem cells (ESCs), only a few of them have been examined in the context of dedifferentiation. In order to improve current protocols of iPSCs generation, it is essential to identify epigenetic drivers and blockages of somatic cell reprogramming.

Repetitive elements and enforced transcriptional repression co-operate to enhance DNA methylation spreading into a promoter CpG-island.

Repression of many tumor suppressor genes in cancer is concurrent with aberrantly increased DNA methylation levels at promoter CpG islands (CGIs). About one-fourth of empirically defined human promoters are surrounded by or contain clustered repetitive elements. It was previously observed that a sharp transition of methylation exists between highly methylated repetitive elements and unmethylated promoter-CGIs in normal tissues. The factors that lead to aberrant CGI hypermethylation in cancer remain poorly understood. Here, we established a site-specific integration system with enforced local transcriptional repression in colorectal cancer cells and monitored the occurrence of initial de novo methylation at specific CG sites adjacent to the CGI of the INSL6 promoter, which could be accelerated by binding a KRAB-containing transcriptional factor. Additional repetitive elements from P16 and RIL (PDLIM4), if situated adjacent to the promoter of INSL6, could confer DNA methylation spreading into the CGI particularly in the setting of KRAB-factor binding. However, a repressive chromatin alone was not sufficient to initiate DNA methylation, which required specific DNA sequences and was integration-site (and/or cell-line) specific. Overall, these results demonstrate a requirement for specific DNA sequences to trigger de novo DNA methylation, and repetitive elements as cis-regulatory factors to cooperate with advanced transcriptional repression in promoting methylation spreading.

Tumor- and circulating-free DNA methylation identifies clinically relevant small cell lung cancer subtypes.

Small cell lung cancer (SCLC) is an aggressive malignancy composed of distinct transcriptional subtypes, but implementing subtyping in the clinic has remained challenging, particularly due to limited tissue availability. Given the known epigenetic regulation of critical SCLC transcriptional programs, we hypothesized that subtype-specific patterns of DNA methylation could be detected in tumor or blood from SCLC patients. Using genomic-wide reduced-representation bisulfite sequencing (RRBS) in two cohorts totaling 179 SCLC patients and using machine learning approaches, we report a highly accurate DNA methylation-based classifier (SCLC-DMC) that can distinguish SCLC subtypes. We further adjust the classifier for circulating-free DNA (cfDNA) to subtype SCLC from plasma. Using the cfDNA classifier (cfDMC), we demonstrate that SCLC phenotypes can evolve during disease progression, highlighting the need for longitudinal tracking of SCLC during clinical treatment. These data establish that tumor and cfDNA methylation can be used to identify SCLC subtypes and might guide precision SCLC therapy.

Genome architecture marked by retrotransposons modulates predisposition to DNA methylation in cancer.

Epigenetic silencing plays an important role in cancer development. An attractive hypothesis is that local DNA features may participate in differential predisposition to gene hypermethylation. We found that, compared with methylation-resistant genes, methylation-prone genes have a lower frequency of SINE and LINE retrotransposons near their transcription start site. In several large testing sets, this distribution was highly predictive of promoter methylation. Genome-wide analysis showed that 22% of human genes were predicted to be methylation-prone in cancer; these tended to be genes that are down-regulated in cancer and that function in developmental processes. Moreover, retrotransposon distribution marks a larger fraction of methylation-prone genes compared to Polycomb group protein (PcG) marking in embryonic stem cells; indeed, PcG marking and our predictive model based on retrotransposon frequency appear to be correlated but also complementary. In summary, our data indicate that retrotransposon elements, which are widespread in our genome, are strongly associated with gene promoter DNA methylation in cancer and may in fact play a role in influencing epigenetic regulation in normal and abnormal physiological states.

GSK-3484862 targets DNMT1 for degradation in cells.

Maintenance of genomic methylation patterns at DNA replication forks by DNMT1 is the key to faithful mitotic inheritance. DNMT1 is often overexpressed in cancer cells and the DNA hypomethylating agents azacytidine and decitabine are currently used in the treatment of hematologic malignancies. However, the toxicity of these cytidine analogs and their ineffectiveness in treating solid tumors have limited wider clinical use. GSK-3484862 is a newly-developed, dicyanopyridine containing, non-nucleoside DNMT1-selective inhibitor with low cellular toxicity. Here, we show that GSK-3484862 targets DNMT1 for protein degradation in both cancer cell lines and murine embryonic stem cells (mESCs). DNMT1 depletion was rapid, taking effect within hours following GSK-3484862 treatment, leading to global hypomethylation. Inhibitor-induced DNMT1 degradation was proteasome-dependent, with no discernible loss of DNMT1 mRNA. In mESCs, GSK-3484862-induced Dnmt1 degradation requires the Dnmt1 accessory factor Uhrf1 and its E3 ubiquitin ligase activity. We also show that Dnmt1 depletion and DNA hypomethylation induced by the compound are reversible after its removal. Together, these results indicate that this DNMT1-selective degrader/inhibitor will be a valuable tool for dissecting coordinated events linking DNA methylation to gene expression and identifying downstream effectors that ultimately regulate cellular response to altered DNA methylation patterns in a tissue/cell-specific manner.

Stem Cell Theory of Cancer: Rude Awakening or Bad Dream from Cancer Dormancy?

To be dormant or not depends on the origin and nature of both the cell and its niche. Similar to other cancer hallmarks, dormancy is ingrained with stemness, and stemness is embedded within dormancy. After all, cancer dormancy is dependent on multiple factors such as cell cycle arrest, metabolic inactivity, and the microenvironment. It is the net results and sum effects of a myriad of cellular interactions, interconnections, and interplays. When we unite all cancer networks and integrate all cancer hallmarks, we practice and preach a unified theory of cancer. From this perspective, we review cancer dormancy in the context of a stem cell theory of cancer. We revisit the seed and soil hypothesis of cancer. We reexamine its implications in both primary tumors and metastatic lesions. We reassess its roles in cell cycle arrest, metabolic inactivity, and stemness property. Cancer dormancy is particularly revealing when it informs us about the mysteries of late relapse, prolonged remission, and second malignancy. It is paradoxically rewarding when it delivers us the promises and power of cancer prevention and maintenance therapy in patient care.

Genotype-to-Phenotype Associations in the Aggressive Variant Prostate Cancer Molecular Profile (AVPC-m) Components.

The aggressive variant prostate cancer molecular profile (AVPC-m), composed of combined defects in TP53, RB1 and PTEN, characterizes a subset of prostate cancers linked to androgen indifference and platinum sensitivity. To contribute to the optimization of the AVPC-m assessment for inclusion in prospective clinical trials, we investigated the status of the AVPC-m components in 28 patient tumor-derived xenografts (PDXs) developed at MDACC. We subjected single formalin-fixed, paraffin-embedded (FFPE) blocks from each PDX to immunohistochemistry (IHC), targeted next-generation genomic sequencing (NGS) and Clariom-S Affymetrix human microarray expression profiling. Standard validated IHC assays and a 10% labeling index cutoff resulted in high reproducibility across three separate laboratories and three independent readers for all tumor suppressors, as well as strong correlations with loss-of-function transcriptional scores (LOF-TS). Adding intensity assessment to labeling indices strengthened the association between IHC results and LOF-TS for TP53 and RB1, but not for PTEN. For TP53, genomic alterations determined by NGS had slightly higher agreement scores with LOF-TS than aberrant IHC, while for RB1 and PTEN, NGS and IHC determinations resulted in similar agreement scores with LOF-TS. Nonetheless, our results indicate that the AVPC-m components can be assessed reproducibly by IHC using various widely available standardized assays.

An Sp1/Sp3 binding polymorphism confers methylation protection.

Hundreds of genes show aberrant DNA hypermethylation in cancer, yet little is known about the causes of this hypermethylation. We identified RIL as a frequent methylation target in cancer. In search for factors that influence RIL hypermethylation, we found a 12-bp polymorphic sequence around its transcription start site that creates a long allele. Pyrosequencing of homozygous tumors revealed a 2.1-fold higher methylation for the short alleles (P<0.001). Bisulfite sequencing of cancers heterozygous for RIL showed that the short alleles are 3.1-fold more methylated than the long (P<0.001). The comparison of expression levels between unmethylated long and short EBV-transformed cell lines showed no difference in expression in vivo. Electrophorectic mobility shift assay showed that the inserted region of the long allele binds Sp1 and Sp3 transcription factors, a binding that is absent in the short allele. Transient transfection of RIL allele-specific transgenes showed no effects of the additional Sp1 site on transcription early on. However, stable transfection of methylation-seeded constructs showed gradually decreasing transcription levels from the short allele with eventual spreading of de novo methylation. In contrast, the long allele showed stable levels of expression over time as measured by luciferase and approximately 2-3-fold lower levels of methylation by bisulfite sequencing (P<0.001), suggesting that the polymorphic Sp1 site protects against time-dependent silencing. Our finding demonstrates that, in some genes, hypermethylation in cancer is dictated by protein-DNA interactions at the promoters and provides a novel mechanism by which genetic polymorphisms can influence an epigenetic state.

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.

Evolution of DNA methylome from precancerous lesions to invasive lung adenocarcinomas.

The evolution of DNA methylome and methylation intra-tumor heterogeneity (ITH) during early carcinogenesis of lung adenocarcinoma has not been systematically studied. We perform reduced representation bisulfite sequencing of invasive lung adenocarcinoma and its precursors, atypical adenomatous hyperplasia, adenocarcinoma in situ and minimally invasive adenocarcinoma. We observe gradual increase of methylation aberrations and significantly higher level of methylation ITH in later-stage lesions. The phylogenetic patterns inferred from methylation aberrations resemble those based on somatic mutations suggesting parallel methylation and genetic evolution. De-convolution reveal higher ratio of T regulatory cells (Tregs) versus CD8 + T cells in later-stage diseases, implying progressive immunosuppression with neoplastic progression. Furthermore, increased global hypomethylation is associated with higher mutation burden, copy number variation burden and AI burden as well as higher Treg/CD8 ratio, highlighting the potential impact of methylation on chromosomal instability, mutagenesis and tumor immune microenvironment during early carcinogenesis of lung adenocarcinomas.